Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 813,205 G→A 7.3% intergenic (+109/+22) PP_0699 → / ← PP_t70 LysE family transporter/tRNA‑Met

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947813,2050GA 7.3% 492.8 / 19.8 164intergenic (+109/+22)PP_0699/PP_t70LysE family transporter/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base G (66/85);  new base A (8/4);  total (75/89)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 169 total aligned reads displayed.

TCGGGCTTGGCCTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTCAACACTCATT  >  NC_002947/813056‑813350
                                                                                                                                                     |                                                                                                                                                 
tCGGGCTTGGCCTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTg                                                                                                                                                   >  1:3235114/1‑150 (MQ=255)
tCGGGCTTGGCCTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTg                                                                                                                                                   >  2:4769897/1‑150 (MQ=255)
tCGGGCTTGGCCTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTg                                                                                                                                                   >  1:2609861/1‑150 (MQ=255)
tCGGGCTTGGCCTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTg                                                                                                                                                   >  2:1684061/1‑150 (MQ=255)
 cGGGCTTGGCCTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGc                                                                                                                                                  >  2:2061393/1‑150 (MQ=255)
 cGGGCTTGGCCTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGc                                                                                                                                                  <  2:4453098/150‑1 (MQ=255)
 cGGGCTTGGCCTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGAGAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGc                                                                                                                                                  <  1:992917/150‑1 (MQ=255)
  gggCTTGGCCTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCg                                                                                                                                                 <  2:2117178/150‑1 (MQ=255)
      ttGGCCTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAATGAAAAAGGGCTTACCTTGCGGTaa                                                                                                                                             >  2:3025262/1‑150 (MQ=255)
      ttGGCCTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTaa                                                                                                                                             >  1:4038221/1‑150 (MQ=255)
      ttGGCCTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCAGTGAGTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTaa                                                                                                                                             <  2:4050684/150‑1 (MQ=255)
          ccTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGccc                                                                                                                                         >  1:2525838/1‑150 (MQ=255)
          ccTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGccc                                                                                                                                         <  2:888125/150‑1 (MQ=255)
          ccTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCGCAAGGTAAGccc                                                                                                                                         <  2:29947/150‑1 (MQ=255)
           cTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCt                                                                                                                                        >  2:2411542/1‑150 (MQ=255)
            tGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGcccgt                                                                                                                                       >  2:3214167/1‑148 (MQ=255)
            tGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCtt                                                                                                                                       >  1:5270417/1‑150 (MQ=255)
              gCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTg                                                                                                                                     <  1:4856827/150‑1 (MQ=255)
                  tgGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCtt                                                                                                                                 <  1:4201250/150‑1 (MQ=255)
                     cGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGtt                                                                                                                              <  2:5194232/150‑1 (MQ=255)
                       gTGACCGGAATATGTTGAGGCGACGATAGCAGCGCCGAAGCCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCGTTGGCTTGTTtg                                                                                                                            >  2:2681258/1‑150 (MQ=255)
                       gTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTTTCtttgcct                                                                                                                            <  2:2126486/150‑6 (MQ=255)
                         gACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTtggt                                                                                                                          <  2:209602/150‑1 (MQ=255)
                          aCCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTtggtg                                                                                                                         >  2:1014195/1‑150 (MQ=255)
                           ccGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTtggtgg                                                                                                                        >  2:3568409/1‑150 (MQ=255)
                             ggAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCt                                                                                                                      <  1:4906142/150‑1 (MQ=255)
                              gAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTa                                                                                                                     >  1:41490/1‑150 (MQ=255)
                                   cGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTCTGGTGGCTACACAg                                                                                                                <  1:1066695/150‑1 (MQ=255)
                                     ttGAGGCGACGATCGCAGCGCCGACACCTGTGATTCGGCGCCCGGCTGCGGCGGGATGGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAACAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAggg                                                                                                              <  2:4347712/150‑1 (MQ=255)
                                      tGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGa                                                                                                             <  1:1448357/150‑1 (MQ=255)
                                         ggCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACtt                                                                                                          <  2:2987091/150‑1 (MQ=255)
                                            gacgaTAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGaa                                                                                                       <  2:1595421/150‑1 (MQ=255)
                                               gaTAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAAccc                                                                                                    <  2:3832171/150‑1 (MQ=255)
                                               gaTAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAAccc                                                                                                    >  1:3666638/1‑150 (MQ=255)
                                               gaTAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAAccc                                                                                                    <  2:2266661/150‑1 (MQ=255)
                                                aTAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAAcccc                                                                                                   <  2:2418859/150‑1 (MQ=255)
                                                   gcagcGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGAATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGa                                                                                                <  1:1779527/150‑1 (MQ=255)
                                                     agcGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCGCAAGGTAAGcc                                                                                                                                          <  1:203859/106‑1 (MQ=255)
                                                        gcCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCa                                                                                           <  1:3910793/150‑1 (MQ=255)
                                                        gcCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCa                                                                                           <  2:3664645/150‑1 (MQ=255)
                                                         ccGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAg                                                                                          >  1:1078014/1‑150 (MQ=255)
                                                             aaCCTGTGATTCGGCGCCCGTCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCATAGAAAAAGGGCTTACCTTGCTGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCAtt                                                                                      <  1:2961817/150‑1 (MQ=255)
                                                             aaCCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAAccc                                                                                                    >  2:4866490/1‑136 (MQ=255)
                                                              aCCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTa                                                                                     <  2:2629995/150‑1 (MQ=255)
                                                                  gtgATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGaa                                                                                 <  1:3079641/150‑1 (MQ=255)
                                                                   tgATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGGCTTGAACCCCGGACCCCAGCATTATGAat                                                                                <  1:3852036/150‑1 (MQ=255)
                                                                   tgATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAat                                                                                <  1:4200971/150‑1 (MQ=255)
                                                                     aTTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAatgc                                                                              >  1:1256880/1‑150 (MQ=255)
                                                                     aTTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAatgc                                                                              >  1:1257087/1‑150 (MQ=255)
                                                                      ttCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAatgct                                                                             >  1:3527025/1‑150 (MQ=255)
                                                                        cGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAatgctat                                                                           >  1:2621968/1‑150 (MQ=255)
                                                                        atccgcCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAatgctat                                                                           <  1:228974/147‑1 (MQ=255)
                                                                         ggCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAatgctggg                                                                          <  1:915591/150‑4 (MQ=255)
                                                                         ggCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAatgctatg                                                                          <  2:2034066/150‑1 (MQ=255)
                                                                           cgcCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCAGCATTATGAATGCTATGctc                                                                       <  2:3132602/150‑1 (MQ=255)
                                                                             cccGGCTGCTGCGGGATTTGTTTGCGCGGGCATAAGCAAATTTTTAGCAGGCAAATAAAAAGGGCTTAACTTTCGGTAATCCCATTGTCTTTTTTGGTGGCTACACAGGGACATTAAACCCGGAACCCACCATTATTAaatctattagat                                                                      >  2:228974/1‑145 (MQ=255)
                                                                              ccGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAATGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTa                                                                     >  2:2743623/1‑150 (MQ=255)
                                                                              ccGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTa                                                                     >  1:4388567/1‑150 (MQ=255)
                                                                                ggCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAAc                                                                   <  2:4448927/150‑1 (MQ=255)
                                                                                   tgctgcGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCaa                                                                >  1:4322113/1‑150 (MQ=255)
                                                                                   tgctgcGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCaa                                                                <  1:3735500/150‑1 (MQ=255)
                                                                                     ctgcGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACTTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACt                                                              <  2:47743/150‑1 (MQ=255)
                                                                                       gcGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGa                                                            <  1:3218963/150‑1 (MQ=255)
                                                                                           gATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTtggt                                                                                                                          >  2:703077/1‑84 (MQ=255)
                                                                                           gATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTtggt                                                                                                                          <  1:703077/84‑1 (MQ=255)
                                                                                             ttggttggCGCGGGCATCATCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCCACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATg                                                      <  1:5130097/150‑1 (MQ=255)
                                                                                                gttggCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAg                                                   >  2:1996300/1‑150 (MQ=255)
                                                                                                    gcgcGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGa                                               >  1:2166025/1‑150 (MQ=255)
                                                                                                    gcgcGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGa                                               <  1:3334632/150‑1 (MQ=255)
                                                                                                       cgGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATgg                                            <  2:4588956/150‑1 (MQ=255)
                                                                                                            aTCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGgcgcgc                                       >  2:1606448/1‑150 (MQ=255)
                                                                                                             tCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCa                                      >  2:336752/1‑150 (MQ=255)
                                                                                                             tCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCa                                      >  1:1286901/1‑150 (MQ=255)
                                                                                                               aGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCatt                                    >  1:3304659/1‑150 (MQ=255)
                                                                                                                gCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGTGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCatta                                   >  2:3484939/1‑150 (MQ=255)
                                                                                                                  aaCTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCattatt                                 <  2:5431179/150‑1 (MQ=255)
                                                                                                                  aaCTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCattatt                                 >  2:3307707/1‑150 (MQ=255)
                                                                                                                   aCTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATtc                                >  2:4302418/1‑150 (MQ=255)
                                                                                                                     tttttCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGGCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATtctc                              >  2:3851103/1‑150 (MQ=255)
                                                                                                                       tttCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCtt                            <  1:4106777/150‑1 (MQ=255)
                                                                                                                         tCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGa                          <  1:2847129/150‑1 (MQ=255)
                                                                                                                             caaacaaGACAAAGGGCTTACCGCAAGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGc                      >  2:2169429/6‑150 (MQ=38)
                                                                                                                              aaacaaGACAAAGGGCTTACCGCAAGGTAAGCCCTTTTTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTa                                                    <  2:2967782/115‑1 (MQ=37)
                                                                                                                              aaacaaGACAAAGGGCTTACCGCAAGGTAAGCCCTTTTTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTa                                                    >  1:2967782/5‑119 (MQ=37)
                                                                                                                               gCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCGTTGTCTTGTTTGGTGGATACACAGGGACTTGAACACCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACTccc                    >  2:3634243/1‑150 (MQ=255)
                                                                                                                                  aaGACAAAGGGCTTACCGCAAGGTAAGCCCTTTTTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCct                 >  1:3216618/1‑150 (MQ=39)
                                                                                                                                    gAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCctct               <  2:4038221/150‑1 (MQ=255)
                                                                                                                                    gAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCctct               <  1:2061393/150‑1 (MQ=255)
                                                                                                                                      aaaaGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGt             <  2:2609861/150‑1 (MQ=255)
                                                                                                                                      aaaaGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGt             <  1:3296220/150‑1 (MQ=255)
                                                                                                                                       aaaGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTc            >  1:4661146/1‑150 (MQ=255)
                                                                                                                                       aaaGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTc            <  1:5367866/150‑1 (MQ=255)
                                                                                                                                       aaaGGGCTTACCGCAAGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTc            >  2:604012/1‑150 (MQ=255)
                                                                                                                                          gggCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTCAac         <  1:3773721/150‑1 (MQ=255)
                                                                                                                                            gCTTACCTTGCGGTAAGAACTTTGTCTTGTTTGGTGGATACACAGGGACTTGAACCCCGCACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTCAacac       >  2:911285/1‑150 (MQ=255)
                                                                                                                                             cTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTCAACACt      >  1:171271/1‑150 (MQ=255)
                                                                                                                                             cTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTCAACACt      >  1:3969990/1‑150 (MQ=255)
                                                                                                                                              ttACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTCAACACTc     >  1:4250887/1‑150 (MQ=255)
                                                                                                                                              ttACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTCAACACTc     >  2:752278/1‑150 (MQ=255)
                                                                                                                                                 ccTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTCAACACTCAtt  <  1:4769897/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
TCGGGCTTGGCCTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTCAACACTCATT  >  NC_002947/813056‑813350

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