Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 813,204 T→C 7.2% intergenic (+108/+23) 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,2040TC 7.2% 512.5 / 19.1 167intergenic (+108/+23)PP_0699/PP_t70LysE family transporter/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base T (70/85);  new base C (8/4);  total (78/89)
Fisher's exact test for biased strand distribution p-value = 2.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

Only 100 of 171 total aligned reads displayed.

TCGGGCTTGGCCTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTCAACACTCATT  >  NC_002947/813056‑813350
                                                                                                                                                    |                                                                                                                                                  
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)
   ggCTTGGCCTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCgg                                                                                                                                                >  2:3318700/1‑150 (MQ=255)
      ttGGCCTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAATGAAAAAGGGCTTACCTTGCGGTaa                                                                                                                                             >  2:3025262/1‑150 (MQ=255)
      ttGGCCTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCAGTGAGTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTaa                                                                                                                                             <  2:4050684/150‑1 (MQ=255)
          ccTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGccc                                                                                                                                         <  2:888125/150‑1 (MQ=255)
          ccTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGccc                                                                                                                                         >  1:2525838/1‑150 (MQ=255)
          ccTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCGCAAGGTAAGccc                                                                                                                                         <  2:29947/150‑1 (MQ=255)
           cTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCt                                                                                                                                        >  2:2411542/1‑150 (MQ=255)
             ggCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCttt                                                                                                                                      <  1:4196640/150‑1 (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)
                                    gTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAgg                                                                                                               >  1:3003355/1‑150 (MQ=255)
                                     ttGAGGCGACGATCGCAGCGCCGACACCTGTGATTCGGCGCCCGGCTGCGGCGGGATGGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAACAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAggg                                                                                                              <  2:4347712/150‑1 (MQ=255)
                                      tGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGa                                                                                                             <  1:1448357/150‑1 (MQ=255)
                                       gAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGAc                                                                                                            >  1:5235770/1‑150 (MQ=255)
                                         ggCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACtt                                                                                                          <  2:2987076/150‑1 (MQ=255)
                                            gacgaTAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGaa                                                                                                       <  2:1595421/150‑1 (MQ=255)
                                               gaTAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAAccc                                                                                                    >  1:3265091/1‑150 (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)
                                                   gcagcGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGa                                                                                                <  1:3214167/150‑1 (MQ=255)
                                                   gcagcGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGAATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGa                                                                                                <  1:1779527/150‑1 (MQ=255)
                                                     agcGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCGCAAGGTAAGcc                                                                                                                                          >  2:203859/1‑106 (MQ=255)
                                                     agcGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCGCAAGGTAAGcc                                                                                                                                          <  1:203859/106‑1 (MQ=255)
                                                        gcCGAAACCTGTGCTTCGGCGCCCGGCTGCTGCGGGATGGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCa                                                                                           <  2:3940751/150‑1 (MQ=255)
                                                        gcCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCa                                                                                           <  1:3910793/150‑1 (MQ=255)
                                                         ccGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAg                                                                                          >  1:1078014/1‑150 (MQ=255)
                                                           gAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCa                                                                                        <  1:2170883/150‑1 (MQ=255)
                                                             aaCCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAAccc                                                                                                    >  2:4866490/1‑136 (MQ=255)
                                                              aCCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTa                                                                                     <  2:2629995/150‑1 (MQ=255)
                                                              aCCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTa                                                                                     <  2:2629957/150‑1 (MQ=255)
                                                                cTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATg                                                                                   >  1:2618029/1‑150 (MQ=255)
                                                                 tgtgATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGa                                                                                  >  1:3842835/1‑150 (MQ=255)
                                                                     aTTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAatgc                                                                              >  1:1256880/1‑150 (MQ=255)
                                                                     aTTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAatgc                                                                              >  1:1257087/1‑150 (MQ=255)
                                                                       tCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAatgcta                                                                            >  1:2913723/1‑150 (MQ=255)
                                                                        cGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAatgctat                                                                           >  1:2621968/1‑150 (MQ=255)
                                                                        atccgcCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAatgctat                                                                           <  1:228974/147‑1 (MQ=255)
                                                                         ggCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTTCCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAatgctatg                                                                          <  1:429491/150‑1 (MQ=255)
                                                                         ggCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAatgctggg                                                                          <  1:915591/150‑4 (MQ=255)
                                                                              ccGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAATGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTa                                                                     >  2:2743623/1‑150 (MQ=255)
                                                                              ccGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTa                                                                     >  1:4388567/1‑150 (MQ=255)
                                                                                   tgctgcGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCaa                                                                <  1:3735500/150‑1 (MQ=255)
                                                                                     ctgcGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACTTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACt                                                              <  2:47743/150‑1 (MQ=255)
                                                                                       gcGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGa                                                            <  1:3219023/150‑1 (MQ=255)
                                                                                       gcGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGa                                                            <  1:3218963/150‑1 (MQ=255)
                                                                                           gATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTtggt                                                                                                                          <  1:703077/84‑1 (MQ=255)
                                                                                           gATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTtggt                                                                                                                          >  2:703077/1‑84 (MQ=255)
                                                                                             ttggttggCGCGGGCATCATCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCCACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATg                                                      <  1:5130097/150‑1 (MQ=255)
                                                                                                gttggCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAg                                                   >  2:1996300/1‑150 (MQ=255)
                                                                                                  tggCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGcc                                                 <  2:4381655/150‑1 (MQ=255)
                                                                                                    gcgcGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGa                                               <  1:3334632/150‑1 (MQ=255)
                                                                                                      gcgGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATg                                             <  2:3040860/150‑1 (MQ=255)
                                                                                                       cgGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATgg                                            <  2:4588956/150‑1 (MQ=255)
                                                                                                            aTCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGgcgcgc                                       >  2:1606448/1‑150 (MQ=255)
                                                                                                             tCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCa                                      >  1:1286901/1‑150 (MQ=255)
                                                                                                               aGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCatt                                    >  1:3304659/1‑150 (MQ=255)
                                                                                                                gCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCatta                                   <  2:366456/150‑1 (MQ=255)
                                                                                                                  aaCTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCattatt                                 <  2:5431179/150‑1 (MQ=255)
                                                                                                                  aaCTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCattatt                                 >  2:3307707/1‑150 (MQ=255)
                                                                                                                   aCTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATtc                                <  2:2913723/150‑1 (MQ=255)
                                                                                                                   aCTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATtc                                >  2:4302418/1‑150 (MQ=255)
                                                                                                                     tttttCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCCTTATtctc                              >  1:376670/1‑150 (MQ=255)
                                                                                                                     tttttCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGGCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATtctc                              >  2:3851103/1‑150 (MQ=255)
                                                                                                                       tttCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCtt                            <  1:4106777/150‑1 (MQ=255)
                                                                                                                         tCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGa                          <  1:2846983/150‑1 (MQ=255)
                                                                                                                             caaacaaGACAAAGGGCTTACCGCAAGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGc                      >  2:2169429/6‑150 (MQ=38)
                                                                                                                              aaacaaGACAAAGGGCTTACCGCAAGGTAAGCCCTTTTTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTa                                                    >  1:2967782/5‑119 (MQ=37)
                                                                                                                               gCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCGTTGTCTTGTTTGGTGGATACACAGGGACTTGAACACCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACTccc                    >  2:3634243/1‑150 (MQ=255)
                                                                                                                                 aaaGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGccccc                  >  2:4666971/1‑150 (MQ=255)
                                                                                                                                 aaaGAAAAAGGGCTTACCTTGCGGTAAGACCTTTGTCTTGTTTGGTGGCTACAAAGGGACTTGAACCACGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATACTCTTGAACGccccc                  >  1:4023432/1‑150 (MQ=255)
                                                                                                                                  aaGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCct                 >  2:3182130/1‑150 (MQ=255)
                                                                                                                                    gAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCctct               <  1:2061393/150‑1 (MQ=255)
                                                                                                                                       aaaGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCt                                                         <  2:3424696/105‑1 (MQ=255)
                                                                                                                                       aaaGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCt                                                         >  1:3424696/1‑105 (MQ=255)
                                                                                                                                       aaaGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTc            >  1:4661146/1‑150 (MQ=255)
                                                                                                                                       aaaGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTc            <  1:5367866/150‑1 (MQ=255)
                                                                                                                                       aaaGGGCTTACCGCAAGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTc            >  2:604012/1‑150 (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)
                                                                                                                                                   ttGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTg                           <  2:2374627/123‑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: