Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,089,805 A→G 5.0% intergenic (+653/+265) 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,089,8050AG 5.0% 631.9 / 15.3 222intergenic (+653/+265)PP_1865/aroHtransposase/3‑deoxy‑7‑phosphoheptulonate synthase
Reads supporting (aligned to +/- strand):  ref base A (143/62);  new base G (6/5);  total (151/71)
Fisher's exact test for biased strand distribution p-value = 3.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

Only 100 of 225 total aligned reads displayed.

TGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGTTCTGAT  >  NC_002947/2089656‑2089954
                                                                                                                                                     |                                                                                                                                                     
   cGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCCGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGGGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGAGGATGGGAGCCGACGGTTGTGGTTACTTTGGCCACGACCa                                                                                                                                                    >  2:5782162/1‑150 (MQ=255)
       cATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGTCGTGGCCAAAGt                                                                                                                                                <  2:1598531/150‑1 (MQ=255)
        aTGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTa                                                                                                                                               >  2:2376601/1‑150 (MQ=34)
         tGGAGGCCCGTCAAGCGCTGGGCCCCCGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTaa                                                                                                                                              <  2:591138/150‑1 (MQ=21)
           gATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCaaa                                                                                                                                                  >  1:7309270/1‑144 (MQ=34)
           gATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCaaa                                                                                                                                                  <  2:7309270/144‑1 (MQ=34)
            aTGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAAccc                                                                                                                                           >  1:226546/1‑150 (MQ=34)
             tGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGTCGTGGCCAAAGTAAccaa                                                                                                                                          <  2:3860458/150‑3 (MQ=255)
             tGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACccg                                                                                                                                          >  2:4000843/1‑150 (MQ=34)
              gCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACccgc                                                                                                                                         >  2:8040493/1‑150 (MQ=34)
                 cGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCgg                                                                                                                                      >  1:5434306/1‑150 (MQ=34)
                  gTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGa                                                                                                                                     >  1:329704/1‑150 (MQ=34)
                   tCACGCGCTGGGCCCCTGGTTGGGGCCTGCTGCGCGGTCCTCCCGGGAGCACGCCGGGTGTGAGGGGACCCCGGAGCGTCGGGTAGGGGCCGGATGTTGGGTGCCGACGGTTTTGGTTACTTTGGCCACGACCAATGTAACCCGCCGGaa                                                                                                                                    <  2:1419418/150‑1 (MQ=11)
                     aaGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAgg                                                                                                                                  >  1:1390243/1‑150 (MQ=34)
                     aaGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAgg                                                                                                                                  >  1:1027707/1‑150 (MQ=34)
                       gcgcTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGc                                                                                                                                >  1:3564066/1‑150 (MQ=34)
                       gcgcTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGc                                                                                                                                >  1:6671464/1‑150 (MQ=34)
                         tgcaggCCCCATCCTGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCgg                                                                                                                              >  1:5686891/5‑150 (MQ=255)
                               ccccTGGTTGGGGCCTGCCGGGGGGTCCCCCGGGTAGCACGCCGGGCGTGCGGGGACCCGGCCGCGTCGGGTCGGGGCGGGTTGTTGGGCGCCGAGGTTTTTGTTTAATTGGGCCACGACAAAAGTAACCCGCCGGAAGGGCGGAAAGGt                                                                                                                        <  2:4845752/150‑1 (MQ=2)
                                cccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTg                                                                                                                       >  2:7187996/1‑150 (MQ=34)
                                cccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTg                                                                                                                       >  1:1244167/1‑150 (MQ=34)
                                cccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTg                                                                                                                       >  2:1031910/1‑150 (MQ=34)
                                 cccggcgggggCGTGCTGGGCGGTCCTCCGGGTAGCACGCCGGGTTTGTGGGGACCCGGGTGGGTTGGGTTGGGGCCGGATGTTGGGAGCCGAGGTTTTTGGTTACTTGGGCCACGACAAAATTAACCCGCCGGAAGGGCGGAAAGGTGa                                                                                                                      <  2:8491320/143‑1 (MQ=2)
                                 ccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGa                                                                                                                      >  1:5475828/1‑150 (MQ=34)
                                  cAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGaa                                                                                                                     >  2:6833257/1‑150 (MQ=34)
                                  cAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGaa                                                                                                                     >  1:7192261/1‑150 (MQ=34)
                                       tGGGGCCTGCAGCGCGGTCCTCCCGGGTGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCGGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGgc                                                                                                                <  2:8320871/150‑1 (MQ=21)
                                       tGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTaa                                                                                                                                              >  1:7416930/1‑120 (MQ=35)
                                       tGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTaa                                                                                                                                              <  2:7416930/120‑1 (MQ=35)
                                            ccTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTa                                                                                                           >  2:4172613/1‑150 (MQ=35)
                                             cTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTAc                                                                                                          >  2:2512195/1‑150 (MQ=35)
                                              tGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACa                                                                                                         >  1:6083081/1‑150 (MQ=37)
                                               gCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACAt                                                                                                        >  1:8907088/1‑150 (MQ=37)
                                                cAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATc                                                                                                       >  1:1796955/1‑150 (MQ=37)
                                                  gcgcGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGAGGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTCCATCGt                                                                                                     >  1:7700504/1‑150 (MQ=21)
                                                   cgcgGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATCGGCGCCTCCATCGtg                                                                                                    >  1:8069326/1‑150 (MQ=21)
                                                     cgGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAggg                                                                                                                                 <  2:9120915/119‑1 (MQ=35)
                                                          ctccGGGGAGCACGCCGGGTGTGTGGGGACCCGGGAGCGTAGCGTTGGGGCGGTTTTTTGGGAGCCGAGGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTgaatgaat                                                                                             <  2:4335203/150‑1 (MQ=11)
                                                          ctccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTgaatgaat                                                                                             >  1:5110161/1‑150 (MQ=37)
                                                           tccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGaatgaatg                                                                                            >  1:3820051/1‑150 (MQ=37)
                                                            ccgggtaGCACGCCGGGTGTGTGGGGACCCCGGCGGGTCGGGTTGGGGCGGGATTATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACAAAATTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGt                                                                                           <  2:4719715/144‑1 (MQ=11)
                                                             ccGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTc                                                                                          >  1:9033152/1‑150 (MQ=37)
                                                              cGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACGTCGTGAATGACTGTCt                                                                                         >  2:8865444/1‑150 (MQ=255)
                                                               gggAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTCGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTa                                                                                        <  2:1474224/150‑1 (MQ=37)
                                                                ggAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTAc                                                                                       >  1:786266/1‑150 (MQ=37)
                                                                    caagcccGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGa                                                                                   >  1:2979879‑M2/8‑150 (MQ=255)
                                                                    cAAGCCCGTTGTGTGGGCACCCCGGACCTTCGGTTGGGGGCGGTATTTTGGGAGCCGCGGTTTTTGGTTAATTTCCCCACCAAAAAACAAACCCGCCGGAAGGGCGTAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGa                                                                                   <  1:4342823/150‑1 (MQ=2)
                                                                     aagcccGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGAt                                                                                  >  1:7827943‑M2/7‑150 (MQ=255)
                                                                     aagcccGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGAt                                                                                  >  1:1628096‑M2/7‑150 (MQ=255)
                                                                     aagcccGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGAt                                                                                  >  2:2180847‑M2/7‑150 (MQ=255)
                                                                       gccgGGGGTGGGGGCACCCGGCTGCGTAGGGTTGGGGCCGTTTGATGGGAGCCGCGGGTTTTGGTTACTTTGGCCACGACCAAAGAAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGc                                                                                <  2:1236654‑M2/146‑1 (MQ=255)
                                                                          cGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGaa                                                                             >  1:2881905‑M2/2‑150 (MQ=255)
                                                                           ggAGTGAGGGGACCCCGGAGGGTAGCGTTGGGGCCGGATGAGGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGaaa                                                                            <  2:8191272/150‑1 (MQ=255)
                                                                           ggAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGaaa                                                                            >  2:6384702/1‑150 (MQ=39)
                                                                             tgtgtGGGGACCCGGTAGCTTTGCTTTGGGGCCGGATTTTGGGAGCCGCGGTTTTTGTTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAgc                                                                          <  2:523282/149‑1 (MQ=14)
                                                                                  ggggACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGcaca                                                                     >  2:1344584/1‑150 (MQ=39)
                                                                                     gACCCCGGAGCGTAGCGTAGGGGCCGGATGTTGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcc                                                                  <  2:2210996/150‑1 (MQ=255)
                                                                                        cccGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCTGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaa                                                               >  2:3471454/1‑150 (MQ=255)
                                                                                        cccGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaa                                                               >  1:8331328/1‑150 (MQ=39)
                                                                                         ccGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaa                                                              >  1:7038187/1‑150 (MQ=39)
                                                                                          cGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaa                                                             >  1:91431/1‑150 (MQ=40)
                                                                                            gcgcgtagcgttGGGGCCGGTTGATGGGAGCCGACGTTTTTGGTTACTTTGGCCACGACCAATGAAACCCGCCGGAAGGGCGGAAAGGTGAGTAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagc                                                           <  1:8150002/148‑1 (MQ=14)
                                                                                              gcgtagcgtcGGGGCCGGTTGTTGGGAGCCGCGGTTTTTGGTTACTTTGGCCACGACCAAATTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagccc                                                         <  2:2692316/150‑1 (MQ=14)
                                                                                                  agcgtagGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGACTAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaa                                                     >  2:7897964/1‑150 (MQ=255)
                                                                                                  agcgtagGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaa                                                     <  1:6833257/150‑1 (MQ=255)
                                                                                                    cgtagGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagc                                                   >  1:6350857/1‑150 (MQ=255)
                                                                                                     gtagGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcc                                                  >  1:7001001/1‑150 (MQ=255)
                                                                                                          ggCCGGTTGTTGGGAGCCGTCGGTTGTGGTTAGTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaacc                                             <  2:3711041/150‑3 (MQ=14)
                                                                                                               gatgatGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagccca                                        >  2:4040340/1‑150 (MQ=255)
                                                                                                                atgatgGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaa                                       >  2:2286684/1‑150 (MQ=255)
                                                                                                                  gatgGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaa                                     >  2:7512664/1‑150 (MQ=255)
                                                                                                                     gggAGCCGACGGTTTTGGTTACTTTGGCCACGAGCAAAGTAACCCGCCGGAAGGGCGGAATGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaatgcccaaaagaccaatagcccataagcc                                  >  2:7874095/1‑150 (MQ=14)
                                                                                                                        tgCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaagcccaa                               <  2:2125268/149‑1 (MQ=255)
                                                                                                                          ccGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGt                                                                                                     <  2:6579787/78‑1 (MQ=39)
                                                                                                                          ccGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGt                                                                                                     >  1:6579787/1‑78 (MQ=39)
                                                                                                                           cGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaagcccaaaag                            >  1:7627713/1‑150 (MQ=255)
                                                                                                                           cGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaagcccaaaag                            >  2:7250069/1‑150 (MQ=255)
                                                                                                                           cGACGGTTTAGGATACTTAGGCCACGACAAAATAAAACCGCCGGAAGGTCTGAAAGGTGAATATGCGCATACATCGTGAATGAATGTCTACATGATGCGAAAGAGCACAAccaaaaatcaaaaaatacaaaaagccaaaaatcccaaaat                            >  2:5358114/1‑149 (MQ=9)
                                                                                                                            gACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACccaaaagcccaaaagcccaaaagcccaaaagcccaaaagc                           >  1:884951/1‑150 (MQ=255)
                                                                                                                            gACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACccaaaagcccaaaagcccaaaagcccaaaagcccaaaagc                           >  1:6613390/1‑150 (MQ=255)
                                                                                                                                 ttttGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAagcccaaaagcccaaaagcccaaaagcccaaaagcccaaa                      >  1:5374935/1‑150 (MQ=255)
                                                                                                                                  tttGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAg                     >  2:1319493/1‑150 (MQ=255)
                                                                                                                                  tttGGTTACTTTGGCCACGACAAAATTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAACGCCCAAAg                     <  1:4440661/150‑1 (MQ=255)
                                                                                                                                   ttGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGACGGGCGGTACGGTGACTAGGCGCCTACATCGTGAATGAAGGTCTACTTGATGCGCAAGCGCACAACCCACAAGCCCAAAAGCCCAAAAGCCCACAAGCCCAAAAGCCCACAGc                    >  2:5990947/1‑150 (MQ=17)
                                                                                                                                      gTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCttt                 <  2:2968202/150‑1 (MQ=255)
                                                                                                                                       ttACTTTGGTCGTGGCCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaa                                     >  2:3889203/1‑129 (MQ=255)
                                                                                                                                       ttACTTTGGTCGTGGCCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTa                >  1:6086455/1‑150 (MQ=255)
                                                                                                                                       ttACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaa                                                             <  1:4619664/105‑1 (MQ=38)
                                                                                                                                       ttACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaa                                                             >  2:4619664/1‑105 (MQ=38)
                                                                                                                                       ttACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTa                >  2:6093635/1‑150 (MQ=255)
                                                                                                                                       ttACTTTGCCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTa                <  1:5197593/150‑1 (MQ=255)
                                                                                                                                           tttGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAgcgc            >  1:4504931/1‑150 (MQ=255)
                                                                                                                                           tttGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAgcgc            >  1:2759997/1‑150 (MQ=255)
                                                                                                                                           tttGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAgcgc            >  2:8536045/1‑150 (MQ=255)
                                                                                                                                               gCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGt        >  2:8084047/1‑150 (MQ=255)
                                                                                                                                               gCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGt        >  2:9024737/1‑150 (MQ=255)
                                                                                                                                                ccACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAATAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGtt       >  2:2944230/1‑150 (MQ=255)
                                                                                                                                                 cACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGTTc      >  2:5056162/1‑150 (MQ=255)
                                                                                                                                                    gACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGTTCTGa   >  2:19126/1‑150 (MQ=255)
                                                                                                                                                    gACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGTTCTGa   >  2:8933947/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGTTCTGAT  >  NC_002947/2089656‑2089954

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: