Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,089,299 C→G 7.5% intergenic (+147/+771) 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,2990CG 7.5% 507.3 / 26.9 159intergenic (+147/+771)PP_1865/aroHtransposase/3‑deoxy‑7‑phosphoheptulonate synthase
Reads supporting (aligned to +/- strand):  ref base C (65/82);  new base G (10/2);  total (75/84)
Fisher's exact test for biased strand distribution p-value = 1.36e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-01

Only 100 of 172 total aligned reads displayed.

GAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGC  >  NC_002947/2089151‑2089448
                                                                                                                                                    |                                                                                                                                                       
 aaaGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCGCGCGCCGGCCGTGCGCCCCGGCGCACGCGCGGCAGGCCccctccca                                                                                                                                                       >  1:3838110/1‑148 (MQ=17)
  aaGGGATGGCCCTTACGGCGGGTCCCTATTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCCGCCCCTACGCTGCGCTCCGGGGTCCCCACGCTCCCCCCGAGCTACCGGGCGCCCCCCCACCAGGCCccagcccag                                                                                                                                                      >  1:2343358/1‑147 (MQ=14)
   aGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATgg                                                                                                                                                    <  2:475175/150‑1 (MQ=255)
   aGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATgg                                                                                                                                                    <  2:2276202/150‑1 (MQ=255)
    gggATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATggg                                                                                                                                                   <  1:3882502/150‑1 (MQ=255)
    gggATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATggg                                                                                                                                                   <  2:1892025/150‑1 (MQ=255)
    gggATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGAGGCCCCATGGATgggg                                                                                                                                                  <  1:2353952/150‑1 (MQ=255)
          gCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCAGCCCTGGGGCCTCa                                                                                                                                            >  2:1445355/1‑150 (MQ=255)
          gCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCGTTGTCCCGGGTGCACCCGCTGGAGGGCCCATCCAATgggggtga                                                                                                                                              >  1:846414/1‑146 (MQ=14)
            ccTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGCGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAgc                                                                                                                                          >  2:820763/1‑150 (MQ=255)
            ccTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCCGGAGCACGCGCTGCAGGCCCCAGCCAAGGGGCCGCAgc                                                                                                                                          >  1:1473566/1‑150 (MQ=255)
             cTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAgcg                                                                                                                                         <  2:4225046/150‑1 (MQ=255)
              ttACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAgcgc                                                                                                                                        <  2:2327709/150‑1 (MQ=255)
              ttACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCGCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCCAGGGGCCTCAgcgc                                                                                                                                        >  1:2629311/1‑150 (MQ=255)
               tACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCt                                                                                                                                       <  1:2268400/150‑1 (MQ=255)
                 cggcggGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATGGATGGGGCCTGCAgcgcgt                                                                                                                                      <  2:3326480/150‑3 (MQ=255)
                   gcggGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGAc                                                                                                                                   <  1:2054993/150‑1 (MQ=255)
                    cggGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACg                                                                                                                                  <  1:3465746/150‑1 (MQ=255)
                       gttccTTTTTGCCACGACAAAAAGGGACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGc                                                                                                                               >  2:3626748/4‑150 (MQ=17)
                        ttccTTTTTGCCACGACAAAAAGGGACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCa                                                                                                                              >  2:2524627/3‑150 (MQ=17)
                           cTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCACCCAACCGCCATGCACCCCGGAGCAACCGCAGCAGGCCCCAACCAAGGGGCCGCAGCCCGTGAAGGGCAGcc                                                                                                                           >  1:2906997/1‑150 (MQ=17)
                            tttttGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCa                                                                                                                          >  2:1764188/1‑150 (MQ=255)
                             ttttGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCAt                                                                                                                         <  2:657411/150‑1 (MQ=255)
                              tttGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCACGCTCCGGCCGTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCCTCCATg                                                                                                                        >  2:2552602/1‑150 (MQ=255)
                               ttGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGc                                                                                                                       <  1:1057033/150‑1 (MQ=255)
                                 gTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATGGATGccc                                                                                                                     <  1:3328455/150‑1 (MQ=255)
                                    gTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTc                                                                                                                  <  2:1693403/150‑1 (MQ=255)
                                     tGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCGCCCGCCGTGCTCCCGGGGGCACGCCCGGCAGGCCCCCAACAAGGGGCCGCAGCCCGTGTGGGGCATCCCAGccccgcc                                                                                                                 >  1:3223542/1‑146 (MQ=14)
                                              ggAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATctccctcc                                                                                                        >  2:2591651/1‑150 (MQ=255)
                                                  ccAAAAACCGTCGGCTCCCATCATCCAGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCaa                                                                                                    <  2:2214327/150‑1 (MQ=255)
                                                   cAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCGCCGGCCTTGCTACCGGGAGCACGCGCTGCAGGCCCCATCCAAGGGGCCGCAGCGCGTGAAGGGCATCCATGCCCGTCATCTCCcgccgccag                                                                                                   >  2:1537634/1‑147 (MQ=17)
                                                       aaCCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCt                                                                                               <  1:458429/150‑1 (MQ=255)
                                                        aCCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTc                                                                                              <  2:4238598/150‑1 (MQ=255)
                                                        aCCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTc                                                                                              <  2:1970099/150‑1 (MQ=255)
                                                          cGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCa                                                                                            <  2:4472605/150‑1 (MQ=255)
                                                            tCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACt                                                                                          <  2:2685684/150‑1 (MQ=255)
                                                                 tCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCAGCTCCCTCCGCCATGCCTCCACTCGGcc                                                                                     >  2:3756963/1‑150 (MQ=255)
                                                                   ccATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCACCCGGGGGCACGCGCGGCAGGCCCCAACCATGGGGCCGCAGCGCATGACCGGCCTCCAAGCCCGTCATATCCCTCCGCAAAGCCACCACACGGccac                                                                                   >  1:2928193/1‑148 (MQ=18)
                                                                    catcatCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGCGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGGCATCTCCCTCCGCACTGCCTCCACTCGGcctcc                                                                                  >  2:4025459/1‑150 (MQ=255)
                                                                    catcatCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCACCGGCCGTGCTCCCGGGCGCACGCCCGGCAGGCCCCAGCCAAGGGGCCGCACCGCATGGCGGGCCACCCAGCCCCTCAACACCCCCCCCAAAGCCACCCCACGGccacc                                                                                  >  1:4486697/1‑150 (MQ=12)
                                                                     atcatcCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGcctcct                                                                                 <  1:1419543/150‑1 (MQ=255)
                                                                        atcCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGaa                                                                              <  2:3868948/150‑1 (MQ=255)
                                                                        atcCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGaa                                                                              <  2:307706/150‑1 (MQ=255)
                                                                           cGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTc                                                                           <  2:3993902/150‑1 (MQ=255)
                                                                             gCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGc                                                                         <  2:3094764/150‑1 (MQ=255)
                                                                                      cTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCTTCTCCCTGCGCAATGCCTCCAGTCGGCCTCCTGAAGTCGCAATCTGCGTc                                                               <  1:689307/150‑1 (MQ=255)
                                                                                            tCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTg                                                          <  2:1436401/150‑1 (MQ=255)
                                                                                                     cccTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACcgc                                                 >  1:1854465/1‑150 (MQ=255)
                                                                                                     cccTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACcgc                                                 >  1:1854551/1‑150 (MQ=255)
                                                                                                     cccTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACcgc                                                 >  2:3009721/1‑150 (MQ=255)
                                                                                                        tCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACcgcgcg                                              >  2:123679/1‑150 (MQ=255)
                                                                                                         cGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCTTCCATGGGGCCTCAGCGCTTGACGGGCATCAATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCgcgcgc                                             <  1:4037956/150‑1 (MQ=255)
                                                                                                            tCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTa                                          >  2:1219361/1‑150 (MQ=255)
                                                                                                            tCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTa                                          <  1:408798/150‑1 (MQ=255)
                                                                                                             ccGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAg                                         <  2:167040/150‑1 (MQ=255)
                                                                                                               ggCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGaa                                       >  1:3885390/1‑150 (MQ=255)
                                                                                                               ggCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCATCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGaa                                       <  2:2844921/150‑1 (MQ=255)
                                                                                                                gCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAg                                      >  1:223828/1‑150 (MQ=255)
                                                                                                                   ttGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCaa                                   <  2:4479482/150‑1 (MQ=255)
                                                                                                                      cTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGAt                                >  1:1955218/1‑150 (MQ=255)
                                                                                                                      cTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGAt                                <  2:2906997/150‑1 (MQ=255)
                                                                                                                      cTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGAt                                <  2:2398571/150‑1 (MQ=255)
                                                                                                                         ccGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCaa                             >  1:2561640/1‑150 (MQ=255)
                                                                                                                         ccGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCaa                             >  1:2823600/1‑150 (MQ=255)
                                                                                                                         ccGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCaa                             >  1:4113647/1‑150 (MQ=255)
                                                                                                                         ccGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCaaaagcca                             >  1:1920603/1‑144 (MQ=255)
                                                                                                                           gggAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCaaaa                           <  1:2524627/150‑1 (MQ=255)
                                                                                                                           gggAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCaaaa                           >  2:2096552/1‑150 (MQ=255)
                                                                                                                             gAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGc                         <  2:812946/150‑1 (MQ=255)
                                                                                                                             gAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGc                         <  2:1875424/150‑1 (MQ=255)
                                                                                                                             gAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGc                         <  2:843009/150‑1 (MQ=255)
                                                                                                                               tcaagcgcTGAGGCCCCATGGATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCgag                      >  2:881204/5‑150 (MQ=17)
                                                                                                                               gCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCga                       <  2:3451044/150‑1 (MQ=255)
                                                                                                                                caagcgcTGAGGCCCCATGGATGGGGCCTGCAGCGCGTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCgag                      >  2:1215293/4‑150 (MQ=17)
                                                                                                                                  cgcgCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGaa                    <  2:3513584/150‑1 (MQ=255)
                                                                                                                                   gcgcTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAg                   >  2:3969432/1‑150 (MQ=255)
                                                                                                                                   gcgcTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAg                   >  2:189505/1‑150 (MQ=255)
                                                                                                                                    cgcTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGc                  <  2:2933205/150‑1 (MQ=255)
                                                                                                                                     gcTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCg                 >  2:4254638/1‑150 (MQ=255)
                                                                                                                                      gctgaGGCCCCATGGATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATAAAAAGCGAGAAGCGa                >  1:4238047/5‑150 (MQ=255)
                                                                                                                                       tGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGaa               >  2:3251950/1‑150 (MQ=255)
                                                                                                                                       tGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGaa               <  2:3895061/150‑1 (MQ=255)
                                                                                                                                        tgaGGCCCCATGGATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAAt              >  1:726544/3‑150 (MQ=255)
                                                                                                                                        tgaGGCCCCATGGATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAAt              >  1:3726216/3‑150 (MQ=255)
                                                                                                                                        gCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAAAGCCAAAGCaagagcaagag              >  1:3656601/1‑143 (MQ=11)
                                                                                                                                         gaGGCCCCATGGATGGGGCCTGCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAAt              >  1:3262239/2‑150 (MQ=255)
                                                                                                                                          aGGCCCCATGGATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCGGGACTACCGCGGGCGGGGACGCAAGAACAAAAGCGAGAAGCGAAtac            >  1:3105022/1‑148 (MQ=18)
                                                                                                                                          aGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTc            >  2:4112034/1‑150 (MQ=255)
                                                                                                                                           ggCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCg           >  1:2951610/1‑150 (MQ=255)
                                                                                                                                            gCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCgg          >  2:1284554/1‑150 (MQ=255)
                                                                                                                                            gCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCgg          <  1:900708/150‑1 (MQ=255)
                                                                                                                                             ccccATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGt         <  2:4172117/150‑1 (MQ=255)
                                                                                                                                             ccccATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGt         <  2:2858261/150‑1 (MQ=255)
                                                                                                                                               ccatccatGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCg           <  1:1025855/146‑1 (MQ=255)
                                                                                                                                               ccatccatGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCg           >  2:1025855/1‑146 (MQ=255)
                                                                                                                                               ccatccatGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGa       >  1:1278277/1‑150 (MQ=255)
                                                                                                                                               ccatccatGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGa       >  1:1207919/1‑150 (MQ=255)
                                                                                                                                                 atccatGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATc     <  2:1067242/150‑1 (MQ=255)
                                                                                                                                                  tccatGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATcg    >  1:231949/1‑150 (MQ=255)
                                                                                                                                                    catGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATcggc  <  1:1288137/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                       
GAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGC  >  NC_002947/2089151‑2089448

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: