Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,955,966 T→G 18.6% intergenic (‑122/+35) fliL ← / ← fliK flagellar protein FliL/flagellar hook‑length control protein FliK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,955,9660TG18.6% 452.3 / 92.1 234intergenic (‑122/+35)fliL/fliKflagellar protein FliL/flagellar hook‑length control protein FliK
Reads supporting (aligned to +/- strand):  ref base T (99/87);  new base G (23/22);  total (123/111)
Fisher's exact test for biased strand distribution p-value = 8.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.64e-01

Only 100 of 234 total aligned reads displayed.

GTCTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCACCCCGGA  >  NC_002947/4955818‑4956114
                                                                                                                                                    |                                                                                                                                                     
gTCTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTc                                                                                                                                                      <  1:2772127/150‑1 (MQ=255)
 tCTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCg                                                                                                                                                     >  1:3697397/1‑150 (MQ=255)
  cTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGc                                                                                                                                                    <  2:142284/150‑1 (MQ=255)
  cTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGc                                                                                                                                                    <  2:650291/150‑1 (MQ=255)
   tttGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGcc                                                                                                                                                   >  1:5009121/1‑150 (MQ=255)
   tttGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGcc                                                                                                                                                   >  1:2619403/1‑150 (MQ=255)
     tGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCcg                                                                                                                                                 >  1:2330261/1‑150 (MQ=255)
      gACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGTAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTACTGGCTTCTTCGCGGGCTCGCCcgc                                                                                                                                                <  2:5158578/150‑1 (MQ=255)
      gACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCcgc                                                                                                                                                >  2:5815768/1‑150 (MQ=255)
       aCTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCcgcg                                                                                                                                               >  1:6807373/1‑150 (MQ=255)
       aCTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCcgcg                                                                                                                                               <  2:4568538/150‑1 (MQ=255)
        cTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGGGCTGGCTTCTTCGCGGGCGAGCCCGCga                                                                                                                                              <  1:4837355/150‑1 (MQ=255)
         tGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGACGGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGAAGGTTTTGGGGTGCTGGCTTCTTCGCGGGCGAGCCCGCgaa                                                                                                                                             <  1:6027705/150‑1 (MQ=255)
         tGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCgaa                                                                                                                                             <  2:2970365/150‑1 (MQ=255)
         tGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCCGCgaa                                                                                                                                             <  2:1548688/150‑1 (MQ=255)
         tGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCCGCgaa                                                                                                                                             <  2:1197957/150‑1 (MQ=255)
         tGCTTCGCTCTTCGCCAGGCCAATAATCCGTCGTCCATCGGGGGTTCAACGGAGAATGGCTGGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCCGCgaa                                                                                                                                             <  2:484910/150‑1 (MQ=255)
          gCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCgaag                                                                                                                                            <  2:6420704/150‑1 (MQ=255)
           cgtCGCTCTTCGCCATGCCAATAATCCGGCGTCCATCGGGGGTTCACCGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCgaaga                                                                                                                                           <  2:204441/148‑1 (MQ=255)
            ttCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTTCAGGGGATGGATTCCAGCAGGTTTTGCGGGGCTGGCTTCTTCGCGGGCGAGCCCGCgaagaa                                                                                                                                          <  1:4291706/150‑1 (MQ=255)
            ttCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCgaagaa                                                                                                                                          <  1:3940462/150‑1 (MQ=255)
            ttCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCCGCgaagaa                                                                                                                                          <  1:3231761/150‑1 (MQ=255)
            ttCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCCGCgaagaa                                                                                                                                          <  2:4991602/150‑1 (MQ=255)
            ttCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTCCGCGGGCTCGCCCGCgaagaa                                                                                                                                          <  1:193467/150‑1 (MQ=255)
                ctctTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAagcca                                                                                                                                      <  2:483596/150‑1 (MQ=255)
                ctctTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAagcca                                                                                                                                      >  1:4471282/1‑150 (MQ=255)
                   ttCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGCAGAagccagcc                                                                                                                                   >  1:1232234/1‑150 (MQ=255)
                   ttCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAagccagcc                                                                                                                                   >  1:6565863/1‑150 (MQ=255)
                    tCGCCATGCCAATAATCCGGCGTCCATCGGGGGTTCAAAGGAGACTGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGCAGAAGCCAGccc                                                                                                                                  <  2:3369955/150‑1 (MQ=255)
                        cATGCAAATACTCCGTCGTCCATCGGGGGTTCAACGGAGAATGGGTTGGCTTGAGCATGTGTTATGCCAGATTTTGCTGGGGATGGATTCCAGCTGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGc                                                                                                                              <  1:1603453/150‑1 (MQ=255)
                             cAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACa                                                                                                                         >  2:389191/1‑150 (MQ=255)
                              aatactCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCCGCgaagaa                                                                                                                                          <  1:2385124/132‑1 (MQ=255)
                              aataatCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGGTCGCCCGCGAAGAAGCCAGCCCACGCCAACAg                                                                                                                        >  1:2853499/1‑150 (MQ=255)
                              aataatCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAg                                                                                                                        >  1:2717356/1‑150 (MQ=255)
                                    ccGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGATATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCa                                                                                                                  <  2:6741741/150‑1 (MQ=255)
                                       tcgtcCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGc                                                                                                               <  1:3402734/150‑1 (MQ=255)
                                       tcgtcCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGc                                                                                                               <  2:5590251/150‑1 (MQ=255)
                                         gtcCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCgt                                                                                                             >  2:5249374/1‑150 (MQ=255)
                                          tcCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCgta                                                                                                            <  2:2369946/150‑1 (MQ=255)
                                          tcCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCgta                                                                                                            <  1:5565377/150‑1 (MQ=255)
                                                gggggTTCAAATGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGGTATCAGGCGTAGTAATc                                                                                                      >  2:5355560/1‑150 (MQ=255)
                                                gggggTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCCCCAAGAAGCCAGCCCACGCCAACACTTATCAGGCGTAGTaaaa                                                                                                      >  1:176172/1‑148 (MQ=255)
                                                  gggTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGCGCCCCCGAAGAAGCCACCCCAGGCCAACAGTTAAAAGGCGTAGTAATCGa                                                                                                    >  1:3166563/1‑150 (MQ=255)
                                                    gTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGAcc                                                                                                  <  2:4375271/150‑1 (MQ=255)
                                                    gTTCAAAGGACAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGAcc                                                                                                  <  2:3260161/150‑1 (MQ=255)
                                                       cAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCa                                                                                               <  2:3953196/150‑1 (MQ=255)
                                                                 tggcttggctTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGa                                                                                     <  2:1163368/150‑1 (MQ=255)
                                                                 tggcttggctTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGa                                                                                     <  2:4387829/150‑1 (MQ=255)
                                                                        gcttgAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGCGCCCGCCAAGAAGCCACCCCAGGCCAACACTTAAAATGCGTGGTAATAGACCCCACAGTAGCCCACAaccac                                                                              >  1:6399107/1‑150 (MQ=255)
                                                                         cttgAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACGGTCGCCGACAaccacc                                                                             >  2:2402706/1‑150 (MQ=255)
                                                                            tcgCATGTGTTATGCAAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACctgc                                                                          <  1:1092269/148‑1 (MQ=255)
                                                                              gCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATAAGGCGTAGTAATCGACCCCACTGTCGCCCACAACCACctgctg                                                                        >  2:2571030/1‑150 (MQ=255)
                                                                               cAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCtgctgc                                                                       <  2:3635522/150‑1 (MQ=255)
                                                                               cAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCtgctgc                                                                       >  2:6679948/1‑150 (MQ=255)
                                                                               cAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCACGCCAACAGTTATAAGGCGTAGTAATAGACCGCACAGTCGCCGACAACCACCagctgc                                                                       >  1:5349270/1‑150 (MQ=255)
                                                                                 aGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCACCCCACGCCCACAGTTATCAGGCGTAGTAATCGACCGCACAGTCGCCGACAACCACCAGCAGCTc                                                                     >  2:5753/1‑150 (MQ=255)
                                                                                     ttATGCCATATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCaaa                                                                 <  1:4336966/150‑1 (MQ=255)
                                                                                     ttATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCACGCCCCCGAAGAAGCCAGCCCACGCCAACACTTATCAGGCGTAGTAATAGACCCCACAGTCGCCGACAAACACCAGCGGCACCaaa                                                                 >  1:6532373/1‑150 (MQ=255)
                                                                                        tGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGc                                                              >  1:3367803/1‑150 (MQ=255)
                                                                                          ccAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCCGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCGCCCAAGGCCg                                                            >  1:2050034/1‑150 (MQ=255)
                                                                                          ccAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCg                                                            >  2:2066359/1‑150 (MQ=255)
                                                                                          ccAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCg                                                            >  1:5612302/1‑150 (MQ=255)
                                                                                          ccAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCAGCTGCTCCAAAGGCCg                                                            >  2:2201333/1‑150 (MQ=255)
                                                                                           cataTTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGATCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCgg                                                           <  2:811489/147‑1 (MQ=255)
                                                                                           cAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCgg                                                           <  2:6152470/150‑1 (MQ=255)
                                                                                           cAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCgg                                                           >  2:5138890/1‑150 (MQ=255)
                                                                                           cAGATTTGGCAGGGGAGGGATTCCAGCTGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCgg                                                           <  2:672231/150‑1 (MQ=255)
                                                                                           aataTTTTGCAGGGGATGGATTCCAGAAGGTTTTGCGTTGCTGGCTTCTTCGCGGGCTCGCCAGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCgg                                                           <  2:4234762/147‑1 (MQ=255)
                                                                                            aGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCCGGCCAACAGTTATAAGGCGTAGTAATCGACCGCACAGTCGCCGACAACCACCAGCAGCTCCAAAGGCCggg                                                          >  1:917694/1‑150 (MQ=255)
                                                                                            aGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGCGCCCCCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAAACGACCGCACAGTCGCCCACAACCCCCAGCTGCTCCAAAGGCCggg                                                          >  2:1693101/1‑150 (MQ=255)
                                                                                                   gCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATc                                                  <  2:741475/150‑1 (MQ=255)
                                                                                                   gCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATc                                                  <  2:1636543/150‑1 (MQ=255)
                                                                                                     aGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGa                                                >  1:1498580/1‑150 (MQ=255)
                                                                                                         gatggatTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGGGCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATAAGGCGTTGTAATCGACCGCACAGTCGCCGACAAACACCGGCTGCTCCAAAGGCCGGGCGCGATCGCCCGc                                            >  1:1427752/1‑150 (MQ=255)
                                                                                                              atTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCCGCCCACGCCAACAGTTATCAGGCGTAGTAATAGACCGCACGGTCGCCGACAACCACCAGCGTATCCAAAGGCCGGGCGCGATCGGCCGCTTCgg                                       >  2:3148328/1‑150 (MQ=255)
                                                                                                              atTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCAGCTGCTCCAAAGGCCGGGCGCTGGTCGACCgacacg                                        >  1:2101433/1‑145 (MQ=255)
                                                                                                                  cagcagGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCCGCCCAGGCCAACAGTTATCAGGCGTAGTAAAAGACCCCACTGTCGCCCACAACCACCTGCAGCTCCAAAGGCCCGGCGCTATCGACCCCATCGGcccc                                   >  2:1917759/1‑150 (MQ=255)
                                                                                                                    gcagGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCAGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCAc                                 >  2:3558355/1‑150 (MQ=255)
                                                                                                                     ctggTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCTAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCAcc                                <  1:713108/148‑1 (MQ=255)
                                                                                                                               ggTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGcc                      <  2:3776881/150‑1 (MQ=255)
                                                                                                                                 tGCTGGCTTCTTCGCGGGCGCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCg                    <  2:4257345/150‑1 (MQ=255)
                                                                                                                                    tGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACACCCCCCTGCTGCTCCAAAGGCCGGGCGCTAGCGACCGCTTCGGCCCCACCCGGCGCGGCCCGgcg                 >  2:3690284/1‑150 (MQ=255)
                                                                                                                                    tGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCAGCTCGGCCCGgcg                 >  2:5386210/1‑150 (MQ=255)
                                                                                                                                     ggCTTCTGCGCGGGCGCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTACTCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGgcgc                <  1:5447344/150‑1 (MQ=255)
                                                                                                                                       cttcttCGCGGGCGAGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGgcgcgc              >  2:278898/1‑150 (MQ=255)
                                                                                                                                        ttcttcGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGcgcc             >  1:5277322/1‑150 (MQ=255)
                                                                                                                                        ttcttcGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGcgcc             >  2:5831223/1‑150 (MQ=255)
                                                                                                                                        ttcttcGCGGGCGAGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGcgcc             >  2:276664/1‑150 (MQ=255)
                                                                                                                                        ttcttcGCGGGCGAGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGcgcc             >  1:1380785/1‑150 (MQ=255)
                                                                                                                                         tcttcGGGGGCCCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGcgccg            <  1:173622/150‑1 (MQ=255)
                                                                                                                                           ttcGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCgccgcc          >  2:5647195/1‑150 (MQ=255)
                                                                                                                                           ttcGCGGGCGAGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTTCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCgccgcc          >  1:4198596/1‑150 (MQ=255)
                                                                                                                                            tcGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCa         >  1:1817825/1‑150 (MQ=255)
                                                                                                                                             cgcgGGCTCGCCCGCCAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTACTAATCGACCGCACTGTCGCCGACAACCACCTGCAGCTCCAAAGGCCGGGCGCTAACGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCAc        >  2:1749740/1‑150 (MQ=255)
                                                                                                                                             cgcgGGCGAGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCAc        >  1:1108929/1‑150 (MQ=255)
                                                                                                                                              gcgGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCAcc       >  2:6216448/1‑150 (MQ=255)
                                                                                                                                                gggCTCGCCCGCGAAGACGCCAGCCCAGGCCACCAGTTATCCGGCGTAGTACTCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCGATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCAcccc     >  1:1462806/1‑150 (MQ=255)
                                                                                                                                                gggCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCAcccc     >  1:517166/1‑150 (MQ=255)
                                                                                                                                                  gCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAACGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCActaaag   <  2:2816238/150‑6 (MQ=255)
                                                                                                                                                   cTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCACCCCGGa  >  2:6143889/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTCTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCACCCCGGA  >  NC_002947/4955818‑4956114

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: