Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,955,967 C→A 13.8% intergenic (‑123/+34) 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,9670CA13.8% 562.7 / 76.6 231intergenic (‑123/+34)fliL/fliKflagellar protein FliL/flagellar hook‑length control protein FliK
Reads supporting (aligned to +/- strand):  ref base C (107/91);  new base A (15/17);  total (122/109)
Fisher's exact test for biased strand distribution p-value = 5.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.74e-01

Only 100 of 237 total aligned reads displayed.

GTCTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCACCCCGGAAA  >  NC_002947/4955818‑4956116
                                                                                                                                                     |                                                                                                                                                     
gTCTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTc                                                                                                                                                       <  1:1400008/150‑1 (MQ=255)
  cTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGc                                                                                                                                                     <  2:650291/150‑1 (MQ=255)
   tttGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGcc                                                                                                                                                    >  1:5009121/1‑150 (MQ=255)
     tGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGTGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCcg                                                                                                                                                  <  2:3584643/150‑1 (MQ=255)
     tGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCcg                                                                                                                                                  >  1:2330261/1‑150 (MQ=255)
      gACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGTAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTACTGGCTTCTTCGCGGGCTCGCCcgc                                                                                                                                                 <  2:5158578/150‑1 (MQ=255)
      gACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCcgc                                                                                                                                                 >  2:5645912/1‑150 (MQ=255)
      gACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCcgc                                                                                                                                                 >  2:5815768/1‑150 (MQ=255)
       aCTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCcgcg                                                                                                                                                >  1:6807373/1‑150 (MQ=255)
        cTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGGGCTGGCTTCTTCGCGGGCGAGCCCGCga                                                                                                                                               <  1:4837355/150‑1 (MQ=255)
         tGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTTGCTTCTTCGCGGGCTCGCCCGCgaa                                                                                                                                              >  1:1436708/1‑150 (MQ=255)
         tGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCgaa                                                                                                                                              <  1:3529542/150‑1 (MQ=255)
         tGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCgaa                                                                                                                                              >  1:2090994/1‑150 (MQ=255)
         tGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCgaa                                                                                                                                              <  2:2970365/150‑1 (MQ=255)
         tGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCgaa                                                                                                                                              <  2:2333178/150‑1 (MQ=255)
         tGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCCGCgaa                                                                                                                                              <  2:1548688/150‑1 (MQ=255)
         tGCTTCGCTCTTCGCCAGGCCAATAATCCGTCGTCCATCGGGGGTTCAACGGAGAATGGCTGGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCCGCgaa                                                                                                                                              <  2:484910/150‑1 (MQ=255)
         ggCTTCGCTCGTCGCCATGCCAATCATCCGGCGTCCATCGGGGGTTCACAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATGTTGCAGGGGATGGAGTCCAGCAGGTTGTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCgaa                                                                                                                                              <  2:4244322/149‑1 (MQ=255)
          gCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCgaag                                                                                                                                             <  2:6420704/150‑1 (MQ=255)
           cgtCGCTCTTCGCCATGCCAATAATCCGGCGTCCATCGGGGGTTCACCGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCgaaga                                                                                                                                            <  2:204441/148‑1 (MQ=255)
            ttCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCgaagaa                                                                                                                                           <  1:3940462/150‑1 (MQ=255)
            ttCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTCCGCGGGCTCGCCCGCgaagaa                                                                                                                                           <  1:193467/150‑1 (MQ=255)
             tCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCCGCGaagaag                                                                                                                                          <  1:1477426/150‑1 (MQ=255)
                ctctTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAagcca                                                                                                                                       >  1:4471282/1‑150 (MQ=255)
                ctctTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAagcca                                                                                                                                       <  2:483596/150‑1 (MQ=255)
                   ttCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGCAGAagccagcc                                                                                                                                    >  1:1232234/1‑150 (MQ=255)
                       ccATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTGGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAgg                                                                                                                                <  2:2653919/150‑1 (MQ=255)
                           gCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCTGCGAAGAAGCCAGCCCAGGCCaa                                                                                                                            >  1:1679149/1‑150 (MQ=255)
                           gCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGCGCCCGCGAAGAAGCCAGCCCAGGCCaa                                                                                                                            >  1:1790314/1‑150 (MQ=255)
                             cAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACa                                                                                                                          >  2:3809155/1‑150 (MQ=255)
                             cAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACa                                                                                                                          <  2:3897629/150‑1 (MQ=255)
                              aataatCCGTCTTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGACAGATTTTGAAGGGGATGGATTCCAGCAGGTTTTGCGGCGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGACAGCCCAGGCCAACAg                                                                                                                         >  2:3270984/1‑150 (MQ=255)
                              aataatCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGGTCGCCCGCGAAGAAGCCAGCCCACGCCAACAg                                                                                                                         >  1:2853499/1‑150 (MQ=255)
                              aataatCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAg                                                                                                                         >  1:2717356/1‑150 (MQ=255)
                                    cggtcgtcCTGGGGGGGTACAACGCAGATTGGTTTGGCTTGCGCATTTTTTATCCATTATTTGGCTGGGTTTGGATCCCCGCTGGTTTTGGGGTGCGGGCTTCTTCGGGGGCCCGCCGGCCATGAAGCCAGCCCAGGCCAACAGTTATCa                                                                                                                   <  1:98967/148‑1 (MQ=255)
                                    ccGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGATATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCa                                                                                                                   <  2:6741741/150‑1 (MQ=255)
                                     cgtcgtCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAg                                                                                                                  >  1:3185251/1‑150 (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:1432423/150‑1 (MQ=255)
                                          tcCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCgta                                                                                                             <  1:5565377/150‑1 (MQ=255)
                                          tcCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCgta                                                                                                             <  2:2369946/150‑1 (MQ=255)
                                             aTCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCgtagta                                                                                                          >  2:5039640/1‑150 (MQ=255)
                                             aTCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCgtagta                                                                                                          >  2:4329593/1‑150 (MQ=255)
                                               cGGGGGTTCAAAGGAGATTGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAAt                                                                                                        <  2:3640543/150‑1 (MQ=255)
                                               cGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAAt                                                                                                        <  2:3697397/150‑1 (MQ=255)
                                                gggggTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAAAGTTATCAGGCGTAGTaaac                                                                                                       >  1:175908/1‑148 (MQ=255)
                                                  gggTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGCGCCCCCGAAGAAGCCACCCCAGGCCAACAGTTAAAAGGCGTAGTAATCGa                                                                                                     >  1:3166563/1‑150 (MQ=255)
                                                     ttCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCg                                                                                                  >  2:525387/1‑150 (MQ=255)
                                                     ttCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCg                                                                                                  >  2:1964028/1‑150 (MQ=255)
                                                          aGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTg                                                                                             >  1:2767029/1‑150 (MQ=255)
                                                              gactggcttggctTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGc                                                                                         <  2:5829295/147‑1 (MQ=255)
                                                                 tggcttggctTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATTGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGa                                                                                      <  2:2881266/150‑1 (MQ=255)
                                                                 tggcttggctTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGa                                                                                      <  2:1163368/150‑1 (MQ=255)
                                                                        gcttgAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGCGCCCGCCAAGAAGCCACCCCAGGCCAACACTTAAAATGCGTGGTAATAGACCCCACAGTAGCCCACAaccac                                                                               >  1:6399107/1‑150 (MQ=255)
                                                                         cttgAGCAAGTGTTATGCCAGCTTTGGCCGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAaccacc                                                                              <  2:2470070/150‑1 (MQ=255)
                                                                            tcgCATGTGTTATGCAAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACctgc                                                                           <  1:1092269/148‑1 (MQ=255)
                                                                            gAGCAAGTGTTATGCCATATTTTGCTGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGATTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACctgc                                                                           <  1:1980145/150‑1 (MQ=255)
                                                                               cAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCtgctgc                                                                        >  2:6679948/1‑150 (MQ=255)
                                                                               cAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCACGCCAACAGTTATAAGGCGTAGTAATAGACCGCACAGTCGCCGACAACCACCagctgc                                                                        >  1:5349270/1‑150 (MQ=255)
                                                                                 aGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTc                                                                      <  2:1647052/150‑1 (MQ=255)
                                                                                 aGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCACCCCACGCCCACAGTTATCAGGCGTAGTAATCGACCGCACAGTCGCCGACAACCACCAGCAGCTc                                                                      >  2:5753/1‑150 (MQ=255)
                                                                                   tgtTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGtt                                                                                                                       >  1:6680810/1‑99 (MQ=255)
                                                                                     ttATGCCATATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCaaa                                                                  <  1:4336966/150‑1 (MQ=255)
                                                                                     ttATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCACGCCCCCGAAGAAGCCAGCCCACGCCAACACTTATCAGGCGTAGTAATAGACCCCACAGTCGCCGACAAACACCAGCGGCACCaaa                                                                  >  1:6532373/1‑150 (MQ=255)
                                                                                        tGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGc                                                               >  1:3367803/1‑150 (MQ=255)
                                                                                         gCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGGAGTAATCGACCGCACTGTCGCCGACACCCACCTGCTGCTCCAAAGGcc                                                              >  2:4385023/1‑150 (MQ=255)
                                                                                          ccAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCg                                                             >  2:2066359/1‑150 (MQ=255)
                                                                                          ccAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCg                                                             >  1:5612302/1‑150 (MQ=255)
                                                                                          ccAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGGCGCCGACAACCACCTGCTGCTCCAAAGGCCg                                                             >  1:2552236/1‑150 (MQ=255)
                                                                                          ccAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACAGTCGCCGACAACCACCTGCTGCGCCAAAGGCCg                                                             >  1:1533316/1‑150 (MQ=255)
                                                                                           cataTTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGATCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCgg                                                            <  2:811489/147‑1 (MQ=255)
                                                                                           cAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCgg                                                            >  2:5138890/1‑150 (MQ=255)
                                                                                           cAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCgg                                                            <  2:6152470/150‑1 (MQ=255)
                                                                                           cAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGATTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCgg                                                            <  2:3258343/150‑1 (MQ=255)
                                                                                           cAGATTTGGCAGGGGAGGGATTCCAGCTGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCgg                                                            <  2:672231/150‑1 (MQ=255)
                                                                                            aGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGCGCCCCCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAAACGACCGCACAGTCGCCCACAACCCCCAGCTGCTCCAAAGGCCggg                                                           >  2:1693101/1‑150 (MQ=255)
                                                                                                   gCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATc                                                    <  2:741475/150‑1 (MQ=255)
                                                                                                            gttttCCAGCTGTTTTTGCGGTGCTGGCTTCTCCGCGGGCGCGCCCGCTAACAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTc                                           <  1:2381246/147‑1 (MQ=255)
                                                                                                            ggatTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTc                                           <  2:4396432/150‑1 (MQ=255)
                                                                                                              atTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCgg                                         >  2:2792310/1‑150 (MQ=255)
                                                                                                                    gcagGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCGATCGACCGCGTCGGCCCCAc                                   >  2:1323264/1‑150 (MQ=255)
                                                                                                                     ctggTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCTAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCAcc                                  <  1:713108/148‑1 (MQ=255)
                                                                                                                                    tGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCAGCTCGGCCCGgcg                   >  2:5386210/1‑150 (MQ=255)
                                                                                                                                     ggCTTCTGCGCGGGCGCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTACTCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGgcgc                  <  1:5447344/150‑1 (MQ=255)
                                                                                                                                        ttcttcGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGcgcc               >  1:5277322/1‑150 (MQ=255)
                                                                                                                                        ttcttcGCGGGCGAGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGcgcc               >  2:276664/1‑150 (MQ=255)
                                                                                                                                        ttcttcGCGGGCGAGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGcgcc               >  1:1380785/1‑150 (MQ=255)
                                                                                                                                         tcttcGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACCGTTATCAGGCGGAGTAATCGACCGCACGGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGcgccg              >  2:4568638/1‑150 (MQ=255)
                                                                                                                                           ttcGCGGTCGAGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCgccgcc            >  2:3512319/1‑150 (MQ=255)
                                                                                                                                           ttcGCGGGCGAGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTTCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCgccgcc            >  1:4198596/1‑150 (MQ=255)
                                                                                                                                           ttcGCGGGCGAGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCgccgcc            >  2:4018378/1‑150 (MQ=255)
                                                                                                                                            tcGCGGGCGAGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCa           >  1:6135321/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)
                                                                                                                                                gggCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCAcccc       >  1:517166/1‑150 (MQ=255)
                                                                                                                                                  gCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCACCCCgg     <  1:4273703/150‑1 (MQ=255)
                                                                                                                                                   cTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCACCCCGGa    >  2:6143889/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTCTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCACCCCGGAAA  >  NC_002947/4955818‑4956116

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: