Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,955,966 T→G 15.2% 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,9660TG15.2% 522.0 / 75.2 236intergenic (‑122/+35)fliL/fliKflagellar protein FliL/flagellar hook‑length control protein FliK
Reads supporting (aligned to +/- strand):  ref base T (110/87);  new base G (18/19);  total (128/108)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

Only 100 of 238 total aligned reads displayed.

GTCTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCACCCCGGAA  >  NC_002947/4955818‑4956115
                                                                                                                                                    |                                                                                                                                                     
gTCTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTc                                                                                                                                                      <  1:7284946/150‑1 (MQ=255)
gTCTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTc                                                                                                                                                      <  2:7093535/150‑1 (MQ=255)
  cTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGc                                                                                                                                                    >  1:6174897/1‑150 (MQ=255)
   tttGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGGTCGcc                                                                                                                                                   >  1:7023114/1‑150 (MQ=255)
   tttGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGcc                                                                                                                                                   >  2:6758130/1‑150 (MQ=255)
     tGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCcg                                                                                                                                                 >  2:4291197/1‑150 (MQ=255)
     tGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCcg                                                                                                                                                 <  1:2860286/150‑1 (MQ=255)
      gACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCcgc                                                                                                                                                <  2:2465332/150‑1 (MQ=255)
      gACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCcgc                                                                                                                                                <  2:6184048/150‑1 (MQ=255)
         tGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCgaa                                                                                                                                             <  2:93220/150‑1 (MQ=255)
         tGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCCGCgaa                                                                                                                                             <  1:3917648/150‑1 (MQ=255)
          gCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGTTTCTTCGCGGGCGAGCCCGCgaag                                                                                                                                            <  2:5826783/150‑1 (MQ=255)
          gCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCCGCgaag                                                                                                                                            <  2:6585319/150‑1 (MQ=255)
            ttCGCTCTTCGCCATGCCAATAATCCGTCGTCCTTCGGGGGTTCAAAGGAGATTGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCgaagaa                                                                                                                                          <  1:3372010/150‑1 (MQ=255)
            ttCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCgaagaa                                                                                                                                          <  1:6758130/150‑1 (MQ=255)
            ttCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCgaagaa                                                                                                                                          <  1:967984/150‑1 (MQ=255)
            ttCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCgaagaa                                                                                                                                          <  2:5361834/150‑1 (MQ=255)
            ttCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCgaagaa                                                                                                                                          >  1:6671538/1‑150 (MQ=255)
            ttCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGAGGGCTCGCCCGCgaagaa                                                                                                                                          <  1:6378610/150‑1 (MQ=255)
                ctctTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGCGCCCGCGAAGAagccc                                                                                                                                      >  1:6921233/1‑149 (MQ=255)
                ctctTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCCGCGAAGAagcca                                                                                                                                      <  1:1385324/150‑1 (MQ=255)
                ctctTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATACCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAagcca                                                                                                                                      <  2:2784083/150‑1 (MQ=255)
                 tctTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAagccag                                                                                                                                     >  1:316933/1‑150 (MQ=255)
                 tctTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAagccag                                                                                                                                     >  2:2941608/1‑150 (MQ=255)
                  ctTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAagccagc                                                                                                                                    <  2:7247673/150‑1 (MQ=255)
                          tGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCa                                                                                                                            >  1:5533068/1‑150 (MQ=255)
                          tGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCa                                                                                                                            >  1:5984321/1‑150 (MQ=255)
                             cAATAATCCGTGGTCCATCGGGGTTTCAAAGGAGAATGGCTTGGCTTGAGCATGTGTTATGCCAGATTTTGCTGGGGATGGATCCCAGCAGGTTTTGCGGTGCTGGCTTCTGCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACa                                                                                                                         <  1:3579655/150‑1 (MQ=255)
                              agtaatCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAg                                                                                                                        >  1:4005943/3‑150 (MQ=255)
                              aataatCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAg                                                                                                                        >  1:2648597/1‑150 (MQ=255)
                                  atCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTAt                                                                                                                    <  2:4847752/150‑1 (MQ=255)
                                   tCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCCGCgaag                                                                                                                                            <  1:5676077/125‑1 (MQ=255)
                                   tCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCCGCgaag                                                                                                                                            >  2:5676077/1‑125 (MQ=255)
                                       tcgtcCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCCACAGTTATCAGGc                                                                                                               >  2:3204304/1‑150 (MQ=255)
                                       tcgtcCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTAGGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGc                                                                                                               <  2:6389527/150‑1 (MQ=255)
                                         gtcCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCgt                                                                                                             >  2:4683543/1‑150 (MQ=255)
                                          tcCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCgta                                                                                                            <  2:4472375/150‑1 (MQ=255)
                                          tcCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCgta                                                                                                            <  1:4283728/150‑1 (MQ=255)
                                           ccATCGGGGTTTCAAAGGAGAAGGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCGGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCgtag                                                                                                           <  1:5357866/150‑1 (MQ=255)
                                             aTCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCCGCgaagaa                                                                                                                                          >  2:947937/1‑117 (MQ=255)
                                                 ggggTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGGAGTAATCg                                                                                                     >  2:1663874/1‑150 (MQ=255)
                                                        aaaGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCACCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCAc                                                                                              >  1:5156693/1‑150 (MQ=255)
                                                         aaGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACt                                                                                             >  2:2819612/1‑150 (MQ=255)
                                                         aaGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGGAGTAATCGACCGCACt                                                                                             >  2:1768364/1‑150 (MQ=255)
                                                                 tggcttggctTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGCACAAGCCAGCCCCCGCCCACACATATCAGGCGGAGTAAACGACCGCACAGTCGCCgc                                                                                     >  1:6482857/1‑149 (MQ=255)
                                                                     ttggcttgAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAACAAGCCAGCCCAGGCCAACAGTTATCATGCGGAGTAATCGACCGCACAGTCGCCGACAac                                                                                 >  1:2160952/1‑150 (MQ=255)
                                                                        gcttgAGCAAGTGTTATGCCAGATTTGGCTGGGTATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAaccac                                                                              <  1:1687992/150‑1 (MQ=255)
                                                                          ttgAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCCACAGTTATAAGGCGTAGTAATCGAACGCAAAGTCGCCGACAaccaccg                                                                            >  2:1237576/1‑149 (MQ=255)
                                                                            gAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACctgc                                                                          >  1:4956095/1‑150 (MQ=255)
                                                                             aGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACctgct                                                                         <  2:1779378/150‑1 (MQ=255)
                                                                              gCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAACCCAGCCCAGGCCCACAGTTATCAGGCGTAGTAATCGACCGCACAGTCGCCGACAACCACctgctg                                                                        >  2:1907801/1‑150 (MQ=255)
                                                                                 aGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCCGGCGCGCCCCCGGAGAAGCCAGCCCACGCCAAAACATATAAGGCGGAGTAAAAGACCCCACGGTCGCCCCCAACCACCagcggcgc                                                                     >  1:2566836/1‑148 (MQ=255)
                                                                                     ttATGCCAGATTTTGCAGGGGTTGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCaaa                                                                 <  1:3647584/150‑1 (MQ=255)
                                                                                     ttATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCaaa                                                                 <  2:4201555/150‑1 (MQ=255)
                                                                                     ttATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGGTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCaaa                                                                 >  2:6867582/1‑150 (MQ=255)
                                                                                      tttGCCAGATTTTGCTGGGTTTGGATACCAGCAGGTTTTGCGGTGCTGGATTCTTCGCGGGCTCCCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAg                                                                <  1:5089469/148‑1 (MQ=255)
                                                                                          ccAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACAGTCGCCGACAACCACCAGCTGCTCCAAAGGCCg                                                            >  2:5012669/1‑150 (MQ=255)
                                                                                          ccAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGGAGTAATCGACCGCACTGTCGCCGACAACCACCAGCTGCTCCAAAGGCCg                                                            >  2:5098850/1‑150 (MQ=255)
                                                                                          ccAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATAAGGCGTAGTAATAGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCg                                                            >  2:1544668/1‑150 (MQ=255)
                                                                                          ccAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCACCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACAGTCGCCGACAACCACCAGCTGCTCCAAAGGCCg                                                            >  1:2272714/1‑150 (MQ=255)
                                                                                           cattttttGCTGGGGTTGGTTTCCAGCTGGTTTTGCGGTGGTGGATTCTTCGGGGGCTCGCCCGCTAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCgg                                                           <  1:5825351/146‑1 (MQ=255)
                                                                                           cAGATTTTGCAGGGGTTGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCgg                                                           <  2:2182795/150‑1 (MQ=255)
                                                                                           cAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCgg                                                           <  2:4013392/150‑1 (MQ=255)
                                                                                           cAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCCCCAACAACCCAGCCCACCCCCACACTTATAAGGCGGAGTAAACGCCCCCCCAGTCGCCCCCAACCACCAGCAGCGCCAAAGGcccg                                                           >  1:1557791/1‑148 (MQ=255)
                                                                                           cAGATTTTGCAGGGGAGGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGTCCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCgg                                                           <  1:5472100/150‑1 (MQ=255)
                                                                                                 ttGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTa                                                     <  2:4594756/150‑1 (MQ=255)
                                                                                                   gCAGGGGATGTATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATc                                                   <  2:2001398/150‑1 (MQ=255)
                                                                                                        ggatggatTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAACAAGCCAGCCCAGGCCAACAGGTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTCTCTACCg                                              >  2:5238452/1‑150 (MQ=255)
                                                                                                         gatggatTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGcc                                                             >  2:7411927/1‑134 (MQ=255)
                                                                                                           tggaaTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTAGCCACCGAAGAACCCAGCCCAGGCCAACAGATATCAGGCGTAGTAAACTACCGCACTGTAGCCGACAACCACCAGGTGGACCAAAGGCCCGGCGCTATCGACCGCtt                                           >  2:6008953/1‑150 (MQ=255)
                                                                                                                   agcagGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCCAAGAAGCCAGCCCAGGCCAACATGTATCAGGCGTAGTAATCGACCGCACTGTCGCCGCCAACCACCGGCTGCGCCAAAGGCCGGGCGCTATCGGCCGCTGCGGccccc                                   >  1:7266835/1‑149 (MQ=255)
                                                                                                                    gctggTTGTGCGGTGGTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCAc                                  <  1:4540204/147‑1 (MQ=255)
                                                                                                                     cagGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCAcc                                 <  2:4149885/150‑1 (MQ=255)
                                                                                                                      agGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCAccc                                <  2:6708542/150‑1 (MQ=255)
                                                                                                                        gTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTg                              >  2:4710407/1‑150 (MQ=255)
                                                                                                                         ttttGCGGTGCTGGCTTCTTCGGGGGCTCGCCCGCGTAGATGCCAGCCCAGGCCAAGAGTTTTCAGGGGTAGTTATGGTCGGGACTGTGGGCGACTACCACCTGCTGGTCCTAAGGGGGGGCGGTATCGACCGCTTCGGCCCCACCCTgg                             >  2:7308856/1‑149 (MQ=255)
                                                                                                                            tggggTGCTGGCTTCTTCGGGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCg                          <  1:80935/147‑1 (MQ=255)
                                                                                                                               ggTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGcc                       <  2:6671538/150‑1 (MQ=255)
                                                                                                                               ggTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGGTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCCCCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCAGCTCGGcc                       >  2:462376/1‑150 (MQ=255)
                                                                                                                                     ggCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGgcgc                 <  2:2216963/150‑1 (MQ=255)
                                                                                                                                        ttcttcGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCgta                                                                                                            >  1:1079593/1‑56 (MQ=255)
                                                                                                                                        ttcttcGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCgta                                                                                                            <  2:1079593/56‑1 (MQ=255)
                                                                                                                                        ttcttcGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGc                                                                                        <  2:3391208/76‑1 (MQ=255)
                                                                                                                                        ttcttcGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGcgcc              >  2:5549746/1‑150 (MQ=255)
                                                                                                                                        ttcttcGCGGGCGAGCCCGCGAAGAAGCCAGCACAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGa                                                 >  1:1007691/1‑115 (MQ=255)
                                                                                                                                         tcttcGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGcgccg             >  1:171434/1‑150 (MQ=255)
                                                                                                                                         tcttcGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGGCGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGcgccg             >  2:2646926/1‑150 (MQ=255)
                                                                                                                                           ttcGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACctgc                                                                          <  1:3869510/87‑1 (MQ=255)
                                                                                                                                           ttcGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACctgc                                                                          >  2:3869510/1‑87 (MQ=255)
                                                                                                                                           ttcGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCgccgcc           >  2:3615705/1‑150 (MQ=255)
                                                                                                                                           ttcGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCgccgcc           >  2:7142064/1‑150 (MQ=255)
                                                                                                                                           ttcGCGGGCGAGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCgccgcc           >  1:5731881/1‑150 (MQ=255)
                                                                                                                                           ttcGCGGGCGAGCCCGCGAAGAAGCCAGCACAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCgccgcc           >  1:441025/1‑150 (MQ=255)
                                                                                                                                            tcGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCa          >  2:208536/1‑150 (MQ=255)
                                                                                                                                            tcGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCa          >  2:6843582/1‑150 (MQ=255)
                                                                                                                                            tcGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCa          >  2:6029484/1‑150 (MQ=255)
                                                                                                                                               cgGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCAccc       >  2:3607537/1‑150 (MQ=255)
                                                                                                                                               cgGGCGAGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCAccc       <  1:1136332/150‑1 (MQ=255)
                                                                                                                                                gggCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCAcccc      <  2:6080744/150‑1 (MQ=255)
                                                                                                                                                    tCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCACCCCGGaa  <  1:4893811/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTCTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCACCCCGGAA  >  NC_002947/4955818‑4956115

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: