Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,171,852 C→T 6.3% intergenic (+19/+27) qseC → / ← ygiZ quorum sensing sensory histidine kinase in two‑component regulatory system with QseB/inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,171,8520CT 6.3% 156.3 / 2.7 64intergenic (+19/+27)qseC/ygiZquorum sensing sensory histidine kinase in two‑component regulatory system with QseB/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (33/26);  new base T (2/2);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.85e-01

CCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAA  >  NC_000913/3171706‑3171995
                                                                                                                                                  |                                                                                                                                               
ccGGACAACCCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCg                                                                                                                                               <  2:876722/149‑1 (MQ=255)
  ggACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGtt                                                                                                                                             >  1:76735/1‑149 (MQ=255)
      aaaCCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaaa                                                                                                                                         >  2:423705/1‑149 (MQ=255)
       aaCCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAc                                                                                                                                        <  2:843069/149‑1 (MQ=255)
       aaCCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAc                                                                                                                                        >  1:559132/1‑149 (MQ=255)
        aCCGCTACCGGCAGCGGGCTTGGGCTATAGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGTAATGCGGGACAAGCTGGATTTGAGGAGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGTGCAGTTAAAACg                                                                                                                                       >  2:771258/1‑149 (MQ=255)
          cGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACgcc                                                                                                                                     >  2:926444/1‑147 (MQ=255)
            cTACCTGCAGCGGGCTTGGGATATCTATTGGCCAGCTAATCGCAAAATTGCATGGCATGAATGTTGAATTTGGGAATCCGGAACAAGGTGGATTTGCGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActc                                                                                                                                   <  2:873228/149‑1 (MQ=255)
            cTACCGGCTGCTGGCTTGGGCTTTCGATTGTCCCGCTAATCGCCAAATTGCATGGCATGAATGTTCAATTTGGGAATGCGGAACAAGGTGTATTTGTGGCGAAGGTAACCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActc                                                                                                                                   <  1:1020570/149‑1 (MQ=255)
                 ggCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATa                                                                                                                              <  1:21309/149‑1 (MQ=255)
                   cAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGATGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGa                                                                                                                            <  1:1150291/149‑1 (MQ=255)
                        ggCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTa                                                                                                                       <  2:57797/149‑1 (MQ=255)
                        ggCTTGGGCTAGCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCGGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTa                                                                                                                       <  1:62141/149‑1 (MQ=255)
                         gCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGTCATGAATGTTGAAATTGGGAATGCGGAACAAGGTTGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                                                                      <  1:541408/149‑1 (MQ=255)
                           ttGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTa                                                                                                                    >  1:604341/1‑149 (MQ=255)
                               gCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTTGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAg                                                                                                                >  1:204460/1‑149 (MQ=255)
                               gCTATCCTTTGTCCAGCGAATACCCAAATTGCATGGCATGAATTTTGAATTTGGGAATGCGCAACAAGGGCGAATTGACGCGAAGGTAAGCTGGAAACTTCTCTGTCTATGCGCGCCGTTAAAACGACTCTCATAGATTTTACTAATAg                                                                                                                <  1:324208/149‑1 (MQ=255)
                                     gATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAat                                                                                                          >  1:647938/1‑149 (MQ=255)
                                     gATTGTCCAGCGAATCGCCAAATTGCATGGCAGGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTATCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAat                                                                                                          >  1:259132/1‑149 (MQ=255)
                                      aTTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTATAAGCAAAata                                                                                                         >  1:985728/1‑149 (MQ=255)
                                         gTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAACGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaa                                                                                                      <  2:654230/149‑1 (MQ=255)
                                              gCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGt                                                                                                 <  2:1010376/149‑1 (MQ=255)
                                                 aaTCGCCAAATTGCATGTCTGGAAGGTTGAATATGGTAATGCGGAACAATGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGAAAAATATAAACCGTccc                                                                                              <  2:1070671/149‑1 (MQ=255)
                                                     gCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTTACTTCTCTGTCTATGAGAGCCGTGAAAACGACTCTCATAGATTTTACGAATAGCAAAATATAAACCGTACCCaaa                                                                                          >  2:713087/1‑149 (MQ=255)
                                                      ccAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaa                                                                                         <  1:948062/149‑1 (MQ=255)
                                                        aaaTTGCATGGCATGAATGTGGAATTTGGGAATGCGCAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaaa                                                                                       <  2:264766/149‑1 (MQ=255)
                                                           ttGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcc                                                                                    >  1:624358/1‑149 (MQ=255)
                                                               atggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcctcca                                                                                >  1:765181/1‑149 (MQ=255)
                                                                    atgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcac                                                                           <  1:423705/149‑1 (MQ=255)
                                                                    atgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcac                                                                           >  2:900790/1‑149 (MQ=255)
                                                                    atgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATACATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcac                                                                           >  1:924841/1‑149 (MQ=255)
                                                                           ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAg                                                                    >  1:1100054/1‑149 (MQ=255)
                                                                           ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATGTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAACCACCAACCACAAACCAg                                                                    >  2:514193/1‑149 (MQ=255)
                                                                            tGAATTTGGGAATGCGGAACAAGGTTGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGACGTTAAAACGACTATCATAGATTTTACTAATTGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGa                                                                   <  1:1076641/149‑1 (MQ=255)
                                                                             gcaTTTGGGAATGCGGAACAAGGTGGATTTGCGGCGAAGGTAAGCTGGTAACTTATCTGTCTATGAGAGCCGTTAAAACGCCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAACGCCACCAACCACAAACCAGAc                                                                  >  1:887786/3‑149 (MQ=255)
                                                                                tttGGTCATGCGGAACAAGGTGTTTTTGAGTCGAAGGAAAGCTGGTAACTTCTCTGTCTAGGAGAGCCGTTAAAACGACACTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaaa                                                               <  2:765181/149‑1 (MQ=255)
                                                                                        tGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                                                       >  1:617232/1‑149 (MQ=255)
                                                                                         tcGGAAGAAGGTGGATTTGATGCGATGGTAAGCTGTTAACTTCTCTGTATATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg                                                      <  1:1145526/148‑1 (MQ=255)
                                                                                         gCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg                                                      >  2:222336/1‑149 (MQ=255)
                                                                                         gCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg                                                      >  2:279130/1‑149 (MQ=255)
                                                                                          ctgAACATGGTGGATTTGAGGCTAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga                                                     <  2:106334/147‑1 (MQ=255)
                                                                                          cGGACCAAGGTGGATTTGATGCGAAGGTAAGCTGGTAACTGCTCTGTCTCTGAGAGCCGTTAAAACGACTCTCATAGATTTTACTTATAGCAAAATATAAACCGTCCACAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga                                                     <  2:523090/149‑1 (MQ=255)
                                                                                             aaCATGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAATAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGaa                                                  <  2:553146/149‑1 (MQ=255)
                                                                                                aagggtgATTTTCGGCGAAGGTAAGCTGTTAACTTATCTGTCTATGAGAGCCGTTAAAACTACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATAc                                               <  2:1075019/143‑1 (MQ=255)
                                                                                                 aGGTGGATTTGAGGCGAAGGTAAGCTGGTATCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAATCCGTCCCCAAAAACGCCACCAACCACAAACCAGACAAATACCCAGAGAGAATACa                                              >  1:1031435/1‑149 (MQ=255)
                                                                                                   gTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACata                                            >  2:85781/1‑149 (MQ=255)
                                                                                                        tttGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAccc                                       >  1:292923/1‑149 (MQ=255)
                                                                                                             ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAAc                                  <  1:411688/149‑1 (MQ=255)
                                                                                                              gCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACt                                 >  1:480721/1‑149 (MQ=255)
                                                                                                                gAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCg                               >  1:336029/1‑149 (MQ=255)
                                                                                                                 aaGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGa                              >  2:23492/1‑149 (MQ=255)
                                                                                                                 aaGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCAAAAACCAGACAAATACCGAGCGAGAATAAATATACCCAAAACTCGa                              >  1:851139/1‑149 (MQ=255)
                                                                                                                    gTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACAGAGCGAGAATACATATACCCAAAACTCGaaaa                           >  2:1106201/1‑147 (MQ=255)
                                                                                                                     tAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAAAAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACAt                          <  2:43176/149‑1 (MQ=255)
                                                                                                                       aGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACAttt                        <  2:608985/149‑1 (MQ=255)
                                                                                                                                   ctctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACGAATAGCAAAATATAAACCGTCCCCAAAAAATCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAATCTCGAACATTTCCCGCATAAaga            >  1:904561/1‑149 (MQ=255)
                                                                                                                                    tctGTCTATGAGAGCCGATAAAACGACTCTCATAGATATTACTAATAGCAAAATATAAACCGACCCCAAAAAATCCACCAACCACAAACCAGACAAATAACGAGCGAGAATACAAATACACAAAACACGAACATTACCCTCATAAagat           >  2:28359/1‑148 (MQ=255)
                                                                                                                                     aaatCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGt          >  2:1097633/4‑149 (MQ=39)
                                                                                                                                      tGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAACCCACCAACCACAAACCAGACAAATACCCAGCGAGAATACATATAACCAAAACTCGAACATTTCCCGCAAAAAGAGtt         >  2:41392/1‑149 (MQ=255)
                                                                                                                                          tATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACACACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCt     >  1:27118/1‑149 (MQ=255)
                                                                                                                                           aTGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCtt    <  2:76735/149‑1 (MQ=255)
                                                                                                                                            tGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTta   <  2:878778/149‑1 (MQ=255)
                                                                                                                                             tagagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAACGTCCCCAAAAAATCAACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAATCCTCGAACATTTCCCGCATAAAGAGTTTTCTtaa  <  1:713087/148‑1 (MQ=255)
                                                                                                                                             gagaGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAAGACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaa  >  2:571819/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
CCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAA  >  NC_000913/3171706‑3171995

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑