Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,639,617 C→T 13.6% intergenic (‑165/‑146) czcR‑III ← / → PP_1439 response regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,639,6170CT13.6% 437.4 / 48.4 170intergenic (‑165/‑146)czcR‑III/PP_1439response regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (86/60);  new base T (14/9);  total (100/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 173 total aligned reads displayed.

TTGAAAGGGGACGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGA  >  NC_002947/1639468‑1639763
                                                                                                                                                     |                                                                                                                                                  
 tGAAAGGGGACGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCcgcg                                                                                                                                                   >  1:7203724/1‑150 (MQ=255)
   aaaGGGGACGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGcc                                                                                                                                                 <  1:8107474/150‑1 (MQ=255)
   aaaGGGGACGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGcc                                                                                                                                                 >  1:4739866/1‑150 (MQ=255)
   aaaGGGGACGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTAGGGCCTCTTCGCGAGTAAACCCGCGcc                                                                                                                                                 >  2:3015486/1‑150 (MQ=255)
       gggACGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTa                                                                                                                                             <  2:2225930/150‑1 (MQ=255)
        ggACGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAg                                                                                                                                            <  2:7222180/150‑1 (MQ=255)
        ggACGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAg                                                                                                                                            >  1:5151823/1‑150 (MQ=255)
         gACGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAgg                                                                                                                                           >  2:1735676/1‑150 (MQ=255)
         gACGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAgg                                                                                                                                           >  2:5758503/1‑150 (MQ=255)
           cGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGaa                                                                                                                                         <  1:402198/150‑1 (MQ=255)
            gATTCTTGACCCTGTAGTGACATCCGGGTCACGCGCCCCAATGTGACGGCCGGAAACAAGGCCGGCGGTGCAGATACCGGATTGGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCGCTTCGCGAGTACACCCGCGCCTGTAGGAAc                                                                                                                                        <  2:5399859/150‑1 (MQ=255)
                cTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGt                                                                                                                                    <  1:5369527/150‑1 (MQ=255)
                   gACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGTTCCTACAGGCGCGGGTTTa                                                                                                                                 <  1:8155474/150‑1 (MQ=38)
                   gACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTa                                                                                                                                 >  2:3163828/1‑150 (MQ=255)
                      ccTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTActc                                                                                                                              <  2:1969292/150‑3 (MQ=37)
                      ccTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGTTCCTACAGGCGCGGGTTTActc                                                                                                                              <  2:2704458/150‑3 (MQ=38)
                      ccTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGTTCCTACAGGCGCGGGTTTActc                                                                                                                              <  1:6341027/150‑3 (MQ=38)
                      ccTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTAccc                                                                                                                              <  2:4291585/150‑1 (MQ=255)
                      ccTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCAGTAGGAACGGGTTTAccc                                                                                                                              >  1:7760618/1‑150 (MQ=255)
                        tGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAGGGCCTGCGGTGCAGATACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCcgcg                                                                                                                           >  1:6987676/1‑150 (MQ=255)
                           aGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGaa                                                                                                                         <  2:4262417/150‑1 (MQ=37)
                           aGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGaa                                                                                                                         <  1:1206055/150‑1 (MQ=37)
                             tGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCTCCtac                                                                                                                                              >  1:4904222/1‑125 (MQ=255)
                             tGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAaga                                                                                                                       <  1:986422/150‑1 (MQ=255)
                                 aTCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGcc                                                                                                                   <  2:8345142/150‑1 (MQ=255)
                                      ggTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAAc                                                                                                              >  1:5826386/1‑150 (MQ=255)
                                         cAAGCGCCCCATTGTTAGGGCATGAAAAATGGCCTGGGGTGCTGATAACGGATTTGTAATTCTGCTGTTATCGGGCTGTGTGTTCGGGCCTCTTCGCGCGTAAACCCGCGCGTGTAGGACCGGGTTTACCCGCGAAGAGGCCCGAACAGt                                                                                                           <  2:7760618/150‑1 (MQ=255)
                                         cAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGt                                                                                                           >  2:4441961/1‑150 (MQ=255)
                                           aGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGTTTTACCCGCGATGAGGCCCGAACAGTcc                                                                                                         <  2:7116751/150‑1 (MQ=255)
                                            gccccccATTGTTCCGGCATGAAAAAAGGCCTGCGGTGCAGATACCGATTTTTTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCGTTTCGTAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCca                                                                                                        <  2:3842278/147‑1 (MQ=38)
                                                cccAATGTTACGGCATGAAAAAAGGCCTGCGGTGTAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGTCCCGAACAGTCcacaca                                                                                                    >  2:4108767/1‑150 (MQ=255)
                                                cccAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCcacaca                                                                                                    >  2:1868343/1‑150 (MQ=255)
                                                cccAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCcacaca                                                                                                    >  1:2613177/1‑150 (MQ=255)
                                                cccAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCcacaca                                                                                                    >  1:7683562/1‑150 (MQ=255)
                                                 ccAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACaa                                                                                                   >  2:8536195/1‑150 (MQ=255)
                                                           ggCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTcg                                                                                                                             <  2:8334211/114‑1 (MQ=21)
                                                           ggCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTcg                                                                                                                             >  1:8334211/1‑114 (MQ=21)
                                                             cATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGATTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCTAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCa                                                                                       <  2:982233/150‑1 (MQ=255)
                                                             cATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGATGAGGCCCGAACAGTCCACACAAGGCCACATTGCa                                                                                       <  2:8624394/150‑1 (MQ=255)
                                                             cATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCa                                                                                       <  2:8839473/150‑1 (MQ=255)
                                                              aTGAAAAATGGCCTGCGGTGAAGATACCGTTTTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGACCGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAg                                                                                      <  2:3927464/150‑1 (MQ=255)
                                                              aTGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGACGAGGCCCGAACAGTCCACACAAGGCCACATTGCAg                                                                                      >  2:4917245/1‑150 (MQ=255)
                                                                 aaaaaaGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGAc                                                                                   >  2:9121140/1‑150 (MQ=255)
                                                                  aaaatGGCCTGCGGTGATTATATCGGTTTTGTATTTCTGGTTTTTTCTGGGTGTGTTTTGGGGACTCTTCGCGAGTAAACCGGCGCTTGTTGGTCCGGGTTTACCCGCGAAGAGGCCGGAACAGTCCACACAAGGCCACATTGCAGGACt                                                                                  <  1:3015486/150‑1 (MQ=255)
                                                                  aaaaaGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACt                                                                                  >  2:4914987/1‑150 (MQ=255)
                                                                                   atctAAGGGATTTGTATTTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGCGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCAcc                                                                 <  1:5124945/147‑1 (MQ=255)
                                                                                    gATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCt                                                                >  1:5727780/1‑150 (MQ=255)
                                                                                      tAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTg                                                              >  2:2483658/1‑150 (MQ=255)
                                                                                      tAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTg                                                              >  2:3219515/1‑150 (MQ=255)
                                                                                          ggATTTGTAATTCTGCTGTTATCTGGCTTTCTGTTCGGGCCTCTTCGCGCGTAAACCCGCGCCTTTTGGAAGGGGTTTACCCGCTAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGt                                                          <  1:615789/150‑1 (MQ=255)
                                                                                          ggATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGt                                                          >  1:6762652/1‑150 (MQ=255)
                                                                                              ttGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTagc                                                      <  2:517012/150‑1 (MQ=255)
                                                                                                gTATTTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTagcca                                                    <  2:1365981/150‑1 (MQ=255)
                                                                                                gTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTagcca                                                    >  1:8699006/1‑150 (MQ=255)
                                                                                                gTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCCCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACACTCCACACAATCCCAAATAGCAGGACTGAAATCCCCGCCTCACCATGCCGTTagcca                                                    >  2:5505958/1‑150 (MQ=38)
                                                                                                    ttCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGccc                                                <  2:3562616/150‑1 (MQ=255)
                                                                                                    ttCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGccc                                                >  2:2409111/1‑150 (MQ=255)
                                                                                                    ttCTGCTGTTATCTGGCTATCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGccc                                                <  2:7586475/150‑1 (MQ=255)
                                                                                                       tgctgTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCg                                             >  2:8032931/1‑150 (MQ=255)
                                                                                                        gctgTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGTCCCGAACAGTCCACACAAGGCCACATTGCATGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGc                                            <  1:4108767/150‑1 (MQ=255)
                                                                                                           gTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCGTGCCGTTAGCCAGCCCTCGCTAg                                         >  1:5233545/1‑150 (MQ=255)
                                                                                                            ttATGTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCGTGTCGTACGGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAgg                                        <  1:1734760/150‑1 (MQ=255)
                                                                                                             tATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGACGGGGTTTACCCGCTAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGa                                       <  1:5758503/150‑1 (MQ=255)
                                                                                                                 tGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGc                                   <  2:6623955/150‑1 (MQ=255)
                                                                                                                    gtgtttgtGCGGGCCTCTTCGCGATAAACCCCGCGCGTGTTTCAAGGGTTTAACCCGCTAAGTGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCgt                                <  1:4977460/149‑1 (MQ=255)
                                                                                                                       tctgtTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTc                             <  2:2181592/150‑1 (MQ=255)
                                                                                                                         tgtTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAAAAGTCCACACAAGGCCACAATGCAGGACTGTAAACCCCCCCTCACCATGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTg                           >  1:5502805/1‑150 (MQ=255)
                                                                                                                          gtTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGGGGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTCATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGa                          >  2:7523505/1‑150 (MQ=255)
                                                                                                                              gggCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCgg                      >  1:6755207/1‑150 (MQ=255)
                                                                                                                               ggCCTATTCGCGAGTAAACCCGCGCATGTAGGAACGGTTTTACCCGCTAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGt                     <  1:199863/150‑1 (MQ=255)
                                                                                                                                gCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGtt                    >  1:323926/1‑150 (MQ=255)
                                                                                                                                  ctctTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGttat                  >  1:8363248/1‑150 (MQ=255)
                                                                                                                                  ctctTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGttat                  >  1:3752551/1‑150 (MQ=255)
                                                                                                                                   tctTCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGttatt                 >  2:2344204/1‑150 (MQ=25)
                                                                                                                                   tctTCGCGATTAAACCCGCGCCTTTTGGAACGGGTTTACCCGCTATGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGttatt                 <  1:585468/150‑1 (MQ=255)
                                                                                                                                     ttCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTAc               >  2:4885608/1‑150 (MQ=37)
                                                                                                                                     ttCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTAc               >  1:912509/1‑150 (MQ=255)
                                                                                                                                     ttCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTAc               >  2:8220543/1‑150 (MQ=255)
                                                                                                                                       cgcgGGTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGAAGAGGCCCGAACAGACCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGcc                                  <  2:7748376/129‑1 (MQ=21)
                                                                                                                                        gcgGGTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAg            >  2:8736456/1‑150 (MQ=37)
                                                                                                                                        gcgGGTAAACCCGTTCCTACAGGCGCGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAg            >  2:5363892/1‑150 (MQ=37)
                                                                                                                                        gcgAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAg            >  1:1567052/1‑150 (MQ=255)
                                                                                                                                        gcgAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAg            >  1:1050028/1‑150 (MQ=255)
                                                                                                                                        gcgAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAg            >  1:8885005/1‑150 (MQ=255)
                                                                                                                                        gcgAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAg            >  2:3488062/1‑150 (MQ=255)
                                                                                                                                         cgAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAgg           <  2:5655770/150‑1 (MQ=255)
                                                                                                                                          gggTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAAGc          >  2:44853/3‑150 (MQ=37)
                                                                                                                                          gggTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGAAGAGGCCCGAACAGACCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGc          >  2:2845572/3‑150 (MQ=37)
                                                                                                                                           aGTAAACCCGCGCCTTTAGTAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCa         <  2:4232536/150‑1 (MQ=255)
                                                                                                                                           aGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCa         <  2:4233341/150‑1 (MQ=255)
                                                                                                                                           aGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCa         >  1:231992/1‑150 (MQ=255)
                                                                                                                                            gTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCaa        >  2:3318923/1‑150 (MQ=37)
                                                                                                                                            gTAAACCCGCGCCTGTAGGAACGGGTTTACTCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCaa        >  2:3678109/1‑150 (MQ=255)
                                                                                                                                              aaaCCCGTTCCTACAGGCGCGGGTTTACTCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAac      >  1:7416245/1‑150 (MQ=38)
                                                                                                                                              aaaCCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAac      >  2:4636728/1‑150 (MQ=255)
                                                                                                                                                 cccGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACg   >  2:6444962/1‑150 (MQ=255)
                                                                                                                                                 cccGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACg   >  1:2166541/1‑150 (MQ=255)
                                                                                                                                                  ccGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGa  >  1:4970299/1‑150 (MQ=255)
                                                                                                                                                  ccGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGa  <  2:3880349/150‑1 (MQ=255)
                                                                                                                                                  ccGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGa  >  1:4970123/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TTGAAAGGGGACGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGA  >  NC_002947/1639468‑1639763

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: