Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,087,878 C→T 5.5% intergenic (+648/+270) PP_1865 → / ← aroH transposase/3‑deoxy‑7‑phosphoheptulonate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,087,8780CT 5.5% 309.4 / 3.2 107intergenic (+648/+270)PP_1865/aroHtransposase/3‑deoxy‑7‑phosphoheptulonate synthase
Reads supporting (aligned to +/- strand):  ref base C (49/52);  new base T (4/2);  total (53/54)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 108 total aligned reads displayed.

GGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGTT  >  NC_002947/2087739‑2088027
                                                                                                                                            |                                                                                                                                                     
ggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCaaa                                                                                                                                             <  2:2302538/150‑1 (MQ=35)
ggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCaaa                                                                                                                                             <  2:1350663/150‑1 (MQ=35)
 gCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACTAAAg                                                                                                                                            >  1:1120367/1‑150 (MQ=21)
  cATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGt                                                                                                                                           >  1:1043022/1‑150 (MQ=34)
   aTGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTa                                                                                                                                          >  1:2088673/1‑150 (MQ=34)
    tGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTaa                                                                                                                                         <  2:2156156/150‑1 (MQ=35)
    tGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTaa                                                                                                                                         <  2:1079628/150‑1 (MQ=35)
         gCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGATGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACccgc                                                                                                                                    >  2:1535226/1‑150 (MQ=25)
           ccGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACcctccg                                                                                                                                  >  2:1421584/1‑150 (MQ=21)
           ccGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACccgccg                                                                                                                                  <  2:2042531/150‑1 (MQ=35)
           ccGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACccgccg                                                                                                                                  >  2:1026149/1‑150 (MQ=35)
            cGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACccgccct                                                                                                                                 >  1:2711554/1‑148 (MQ=25)
              tCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCaaa                                                                                                                                             <  2:1106324/136‑1 (MQ=35)
              tCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGGTTTGGTTACTTTGGCCACGACCaaa                                                                                                                                             >  1:1106324/1‑136 (MQ=17)
                   cgcTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCg                                                                                                                          <  1:172072/150‑1 (MQ=35)
                   cgcTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCg                                                                                                                          >  1:958154/1‑150 (MQ=34)
                     tgaggCCCCATGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCgg                                                                                                                         >  2:1598427/4‑150 (MQ=255)
                          ccccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGt                                                                                                                   >  2:1904526/1‑150 (MQ=35)
                          ccccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGt                                                                                                                   <  2:202772/150‑1 (MQ=35)
                            ccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGa                                                                                                                 >  1:820836/1‑150 (MQ=35)
                            ccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGa                                                                                                                 >  1:942784/1‑150 (MQ=35)
                             cAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTGGCCCGCCGGAAGGGCGGAAAGGTGaa                                                                                                                <  1:1237054/150‑1 (MQ=25)
                             cAGGATGGCGCCTACAGAGCTGACCCTCCAGGAACATGTCCTGAGTGAGGGGACCCCGGAGCGTAGCGTAGCTGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTACCCCGCCGGAAGGGCGGAAATGTGaa                                                                                                                <  1:1300626/150‑1 (MQ=11)
                              aGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGTGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGAGAAt                                                                                                               >  2:608145/1‑150 (MQ=255)
                                tcctGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTGTTGGTTACTGTGGCCACGACCATAGGAACCCGCCGGTAGGGCGGAAAGGTGAATa                                                                                                              >  2:2464334/4‑150 (MQ=11)
                                    ggggCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGgcg                                                                                                          >  2:830022/1‑150 (MQ=35)
                                     gggCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGgcgc                                                                                                         >  1:899819/1‑150 (MQ=35)
                                      ggCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAg                                                                                                                     >  2:1579984/1‑137 (MQ=37)
                                      ggCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAg                                                                                                                     <  1:1579984/137‑1 (MQ=37)
                                           gCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACAt                                                                                                   <  1:1904526/150‑1 (MQ=38)
                                              gcgcGGTCCTCCCGGGAGCAAGGCCGGAGTGAGGGGCGTCGTGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAGCCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGt                                                                                                <  1:2464334/150‑1 (MQ=11)
                                                 cgGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCCCGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGtgaa                                                                                             <  2:1982631/150‑1 (MQ=255)
                                                 cgGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGtgaa                                                                                             <  2:1393487/150‑1 (MQ=38)
                                                   gTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGGCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGtgaatg                                                                                           >  1:757106/1‑150 (MQ=255)
                                                    tcctccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAggg                                                                                                                            >  1:951885/1‑116 (MQ=38)
                                                    tcctccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAggg                                                                                                                            <  2:951885/116‑1 (MQ=38)
                                                    tcctccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGtgaatga                                                                                          >  2:2310415/1‑150 (MQ=39)
                                                     cctccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTACCCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGtgaatgaa                                                                                         <  2:2088673/150‑1 (MQ=255)
                                                      ctccCGGGAGCAAGCCCGGAGTTAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGAACACCGTAGTGTGGCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTgaatgaat                                                                                        >  2:2385530/1‑150 (MQ=14)
                                                      ctccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTgaatgaat                                                                                        <  1:760977/150‑1 (MQ=39)
                                                      ctccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTgaatgaat                                                                                        >  2:1627730/1‑150 (MQ=39)
                                                       tccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTgaatgaat                                                                                        <  1:1321907/149‑1 (MQ=39)
                                                       tccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTgaatgaat                                                                                        >  2:1321907/1‑149 (MQ=39)
                                                         ccGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTc                                                                                     >  2:554259/1‑150 (MQ=37)
                                                         ccGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTc                                                                                     >  2:3019434/1‑150 (MQ=39)
                                                         ccGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTc                                                                                     >  1:2756395/1‑150 (MQ=39)
                                                         ccGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTc                                                                                     >  2:1194304/1‑150 (MQ=39)
                                                         ccGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTc                                                                                     >  2:880521/1‑150 (MQ=38)
                                                           gggAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGGCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTa                                                                                   <  2:757106/150‑1 (MQ=255)
                                                                 aaGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGAt                                                                             <  2:354999/150‑1 (MQ=39)
                                                                   gCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGt                                                                                                                   >  1:801563/1‑110 (MQ=38)
                                                                   gCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGt                                                                                                                   <  2:801563/110‑1 (MQ=38)
                                                                   gCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTgaatgaat                                                                                        <  2:309092/137‑1 (MQ=39)
                                                                   gCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTgaatgaat                                                                                        >  1:309092/1‑137 (MQ=40)
                                                                   gCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGc                                                                           >  2:2802459/1‑150 (MQ=38)
                                                                     ccGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTAc                                                                                  <  2:516602/141‑1 (MQ=38)
                                                                     ccGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTAc                                                                                  >  1:516602/1‑141 (MQ=38)
                                                                                ggACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCTCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAc                                                              <  1:1421584/150‑1 (MQ=255)
                                                                                ggACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGc                                                                                                                           >  2:1699707/1‑89 (MQ=40)
                                                                                ggACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGc                                                                                                                           <  1:1699707/89‑1 (MQ=40)
                                                                                   ccccGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTCCTTTGGCCACGACAAAAGTTACCCGCCGGAAGGGCGGAAAGTTTAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCTCACAAccca                                                           <  2:1120367/150‑1 (MQ=14)
                                                                                    cccGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaa                                                          >  2:2085934/1‑150 (MQ=255)
                                                                                      cGGAGCGTAGTGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaa                                                        <  1:608145/150‑1 (MQ=255)
                                                                                      cGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaa                                                        <  1:830022/150‑1 (MQ=255)
                                                                                      cGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaa                                                        <  1:1194304/150‑1 (MQ=255)
                                                                                         agcgtagcgtAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcc                                                     <  1:2310415/150‑1 (MQ=255)
                                                                                             tagcgtagGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGTCGTGGCCAAAGTaa                                                                                                                                         <  2:1822707/62‑1 (MQ=14)
                                                                                             tagcgtagGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGTCGTGGCCAAAGTaa                                                                                                                                         >  1:1822707/1‑62 (MQ=14)
                                                                                             tagcgtagGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAACGCACAAcccaaaagcccaaa                                                 >  2:1354049/1‑150 (MQ=255)
                                                                                                cgtagGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTACCCCGCCGGAAGGGCGGAAAGGTGAATCGGCGCCTACATCGTGAATGAACGTCTACGTGATGCGAAAGCGCACAAcccaaaagccccaaagc                                              <  1:1598427/150‑1 (MQ=14)
                                                                                                cgtagGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagc                                              <  1:1627730/150‑1 (MQ=255)
                                                                                                 gtagGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcc                                             <  1:880521/150‑1 (MQ=255)
                                                                                                     gggCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaa                                         >  1:1879770/1‑150 (MQ=255)
                                                                                                       gCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaag                                       <  1:2802459/150‑1 (MQ=255)
                                                                                                         cGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCCCCGTCACCCGCCCGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcc                                     <  2:1879770/150‑1 (MQ=14)
                                                                                                         cGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATAGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcc                                     >  1:1038554/1‑150 (MQ=255)
                                                                                                                tgGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccagtcttccaaaagcccaaaagcccaaaagc                              <  1:2385530/150‑1 (MQ=17)
                                                                                                                tgGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaagc                              >  1:1379497/1‑150 (MQ=255)
                                                                                                                 gggAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcc                                                             <  1:3034756/118‑1 (MQ=255)
                                                                                                                 gggAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcc                                                             >  2:3034756/1‑118 (MQ=40)
                                                                                                                  ggAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaa                                                        <  2:2558079/122‑1 (MQ=255)
                                                                                                                  ggAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaa                                                        >  1:2558079/1‑122 (MQ=255)
                                                                                                                      ccGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaagcccaaaa                        <  1:2085934/150‑1 (MQ=255)
                                                                                                                      ccGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaagcccaaaa                        <  2:899819/150‑1 (MQ=255)
                                                                                                                        gACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGtgaat                                                                                            >  1:392021/1‑80 (MQ=255)
                                                                                                                        gACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGtgaat                                                                                            <  2:392021/80‑1 (MQ=255)
                                                                                                                        gACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaagcccaaaag                       >  2:2642472/1‑149 (MQ=255)
                                                                                                                        gACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaagcccaaaag                       <  1:2642472/149‑1 (MQ=255)
                                                                                                                        gACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACccaaaagcccaaaagcccaaaagcccaaaagcccaaaagc                      >  1:1048821/1‑150 (MQ=255)
                                                                                                                         aCGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCcaaaagcccaaaagcccaaaagcccaaaagcccaaaagcc                     >  1:282137/1‑150 (MQ=255)
                                                                                                                              tttGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaa                                        <  1:439184/126‑1 (MQ=255)
                                                                                                                              tttGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaa                                        >  2:439184/1‑126 (MQ=255)
                                                                                                                               ttGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaa                                        <  2:483941/125‑1 (MQ=255)
                                                                                                                               ttGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaa                                        >  1:483941/1‑125 (MQ=255)
                                                                                                                                 ggTTACTTTGGTCGTGGCCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCtt             >  2:1336382/1‑150 (MQ=18)
                                                                                                                                 ggTTACTTTGGTCGTGGCCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCtt             >  2:3053694/1‑150 (MQ=255)
                                                                                                                                   ttACTTTGGTCGTGGCCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaa                                        <  1:356048/121‑1 (MQ=14)
                                                                                                                                   ttACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTa           <  2:282137/150‑1 (MQ=255)
                                                                                                                                       tttGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGaaa                                                                       <  2:2612793/86‑1 (MQ=40)
                                                                                                                                            ccACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGtt  >  2:2645389/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                                     
GGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGTT  >  NC_002947/2087739‑2088027

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