Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,087,377 C→G 13.6% 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,087,3770CG13.6% 299.2 / 33.2 103intergenic (+147/+771)PP_1865/aroHtransposase/3‑deoxy‑7‑phosphoheptulonate synthase
Reads supporting (aligned to +/- strand):  ref base C (40/49);  new base G (6/8);  total (46/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 108 total aligned reads displayed.

GAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGC  >  NC_002947/2087229‑2087526
                                                                                                                                                    |                                                                                                                                                       
  aaGGGATGGCCCTTACGGCGGGACCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATg                                                                                                                                                      >  1:1829435/1‑150 (MQ=255)
       aTGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATGGATGGGGcc                                                                                                                                                 <  2:501390/150‑1 (MQ=255)
           cccTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAg                                                                                                                                            >  1:1405718/1‑150 (MQ=255)
              ttACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGAGGCCCCATGGATGGGGCCTGCAgcgc                                                                                                                                         <  2:320921/150‑1 (MQ=17)
                aCGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCtt                                                                                                                                       <  2:471344/150‑1 (MQ=255)
                  ggcggGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCt                                                                                                                                                >  2:1069968/1‑140 (MQ=255)
                  ggcggGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCt                                                                                                                                                <  1:1069968/140‑1 (MQ=255)
                  ggcggGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGa                                                                                                                                     <  1:929151/150‑1 (MQ=255)
                  ggcggGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGa                                                                                                                                     <  1:15272/150‑1 (MQ=255)
                     gggTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCACCTGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACgg                                                                                                                                  <  1:23584/150‑1 (MQ=255)
                          ccTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATgg                                                                                                                                                     >  2:250100/1‑127 (MQ=255)
                           cTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCCTGGGGCCTCAGCGCTTGACGGGCAtgt                                                                                                                            <  2:1842356/150‑3 (MQ=255)
                           cTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCAtgg                                                                                                                            >  2:1709854/1‑148 (MQ=255)
                             ttttGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCAt                                                                                                                          <  2:1405718/150‑1 (MQ=255)
                               ttGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGc                                                                                                                        <  2:1746046/150‑1 (MQ=255)
                                  tCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGAGGCCCCATGGATGGGGCCtgcag                                                                                                                                            <  2:1633720/127‑5 (MQ=17)
                                  tCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGAGGCCCCATGGATGGGGCCtgcag                                                                                                                                            >  1:1633720/1‑123 (MQ=17)
                                    gTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGc                                                                                                                        <  1:756065/145‑1 (MQ=255)
                                    gTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGc                                                                                                                        >  2:756065/1‑145 (MQ=255)
                                    gTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTc                                                                                                                   <  1:2231728/150‑1 (MQ=255)
                                     tGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATGGATGCCCGTCa                                                                                                                  <  1:530324/150‑1 (MQ=255)
                                      ggCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATGGATGCCCGTCaa                                                                                                                 <  1:1709854/150‑2 (MQ=255)
                                      ggCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCAt                                                                                                                 <  1:940498/150‑1 (MQ=255)
                                        cAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCAtct                                                                                                               >  1:1213889/1‑150 (MQ=255)
                                        cAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCAtct                                                                                                               <  1:1897327/150‑1 (MQ=255)
                                        cAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCAtct                                                                                                               >  2:2092965/1‑150 (MQ=255)
                                         aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATgg                                                                                                                                                     <  2:1562278/112‑1 (MQ=255)
                                          aaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATctcc                                                                                                             >  1:131959/1‑150 (MQ=255)
                                              ggAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATctccctcc                                                                                                         <  1:2092965/150‑1 (MQ=255)
                                                aaCCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGc                                                                                                       >  2:1255357/1‑150 (MQ=255)
                                                  ccAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCa                                                                                                                  >  1:777176/1‑137 (MQ=255)
                                                   cAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCACCCATGGGTCCTCCGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAt                                                                                                    <  2:2393389/150‑1 (MQ=255)
                                                    aaaaaCCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATg                                                                                                   >  2:595881/1‑150 (MQ=255)
                                                        aCCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTc                                                                                               <  1:96099/150‑1 (MQ=255)
                                                          cGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATgg                                                                                                                                                     >  1:1528140/1‑95 (MQ=255)
                                                          cGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATgg                                                                                                                                                     <  2:1528140/95‑1 (MQ=255)
                                                           gTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCAc                                                                                            <  1:397582/150‑1 (MQ=255)
                                                             cGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTc                                                                                          <  1:1841152/150‑1 (MQ=255)
                                                              ggCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCg                                                                                         <  2:295358/150‑1 (MQ=255)
                                                                ctgccATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTGGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCCTGGGGCCTCAGCGCTTGACGGGCAGCCATGCCCGTCATCTCCGTCCGCAATGCGTCCACTCGGc                                                                                       >  1:61204/4‑150 (MQ=17)
                                                                  cccATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGcct                                                                                     <  2:20829/150‑1 (MQ=255)
                                                                   ccATCATTCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTGCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCAATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGcctc                                                                                    <  1:2301361/150‑1 (MQ=255)
                                                                   ccATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGcctc                                                                                    >  1:1262258/1‑150 (MQ=255)
                                                                    catcatCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGcctcc                                                                                   >  1:1607098/1‑150 (MQ=255)
                                                                    catcatCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGcctcc                                                                                   >  2:2262881/1‑150 (MQ=255)
                                                                    catcatCCGGCCCCTACGCTACGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGcctcc                                                                                   >  1:2047720/1‑150 (MQ=255)
                                                                     atcatcCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGa                                                                                                                                     <  2:1157847/99‑1 (MQ=255)
                                                                     atcatcCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGa                                                                                                                                     >  1:1157847/1‑99 (MQ=255)
                                                                     atcatcCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGcctcct                                                                                  >  2:2341772/1‑150 (MQ=255)
                                                                        atcCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGaa                                                                               <  1:2025463/150‑1 (MQ=255)
                                                                        atcCGGCCCCTACGCTGCGCTCCGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGAGGCCCCATGGATGGGGCCTGCAgcgc                                                                                                                                         >  1:244185/1‑91 (MQ=25)
                                                                        atcCGGCCCCTACGCTGCGCTCCGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGAGGCCCCATGGATGGGGCCTGCAgcgc                                                                                                                                         <  2:244185/91‑1 (MQ=21)
                                                                           cGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTc                                                                            >  1:1430947/1‑150 (MQ=255)
                                                                                   aCGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTg                                                                   <  2:1213889/150‑1 (MQ=255)
                                                                                    cGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGc                                                                  <  1:595881/150‑1 (MQ=255)
                                                                                    cGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGc                                                                  <  1:573826/150‑1 (MQ=255)
                                                                                    cGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGc                                                                  <  2:1262258/150‑1 (MQ=255)
                                                                                       tGCGCTCCGGGGTCCCCTCGCTCCGGCCTTTCTCCCGGGAGCACGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTgc                                                                                                                                     <  2:1791613/81‑2 (MQ=18)
                                                                                       tGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGAGGCCCCATGGATGGGGCCTGCAgcgcgtgc                                                                                                                                     >  1:1791613/1‑77 (MQ=18)
                                                                                          gcTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGcc                                                            <  2:2047717/150‑1 (MQ=255)
                                                                                          gcTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGcc                                                            >  1:2169917/1‑150 (MQ=255)
                                                                                              cGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACgg                                                                                                                                  >  1:2227855/1‑77 (MQ=255)
                                                                                              cGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACgg                                                                                                                                  <  2:2227855/77‑1 (MQ=255)
                                                                                              cGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCt                                                                    >  2:575274/1‑138 (MQ=255)
                                                                                              cGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCt                                                                    <  1:575274/138‑1 (MQ=255)
                                                                                                     cccTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGACTAATGAAGTCGCAATCTTCCTCACCTGTACTACcgc                                                 >  2:1344011/1‑150 (MQ=17)
                                                                                                       cTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTCAcgcgt                                               >  2:2177394/1‑149 (MQ=255)
                                                                                                       cTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACcgcgc                                               >  1:1517530/1‑150 (MQ=255)
                                                                                                          gCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCt                                            >  1:895716/1‑150 (MQ=255)
                                                                                                            tCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATGTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTa                                          >  1:762006/1‑150 (MQ=255)
                                                                                                            tCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTa                                          >  2:438724/1‑150 (MQ=255)
                                                                                                              cGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGa                                        >  2:1436668/1‑150 (MQ=255)
                                                                                                                  cTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGAc                                                       >  1:2216325/1‑131 (MQ=255)
                                                                                                                  cTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGAc                                                       <  2:2216325/131‑1 (MQ=255)
                                                                                                                   ttGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCaa                                   <  2:131959/150‑1 (MQ=255)
                                                                                                                     gCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGa                                 <  1:275681/150‑1 (MQ=255)
                                                                                                                         ccGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGCATACCGAGCGCt                                            <  2:1844259/135‑1 (MQ=255)
                                                                                                                         ccGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCt                                            >  1:1844259/1‑135 (MQ=255)
                                                                                                                             gAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGc                         >  2:1962274/1‑150 (MQ=255)
                                                                                                                                       tGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGaa               <  2:1607098/150‑1 (MQ=255)
                                                                                                                                        tgaGGCCCCATGGATGGGGCCTGCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGaa               >  1:322514/3‑150 (MQ=18)
                                                                                                                                        tgaGGCCCCATGGATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAg                                  <  2:2348254/128‑1 (MQ=255)
                                                                                                                                        tgaGGCCCCATGGATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAg                                  >  1:2348254/3‑130 (MQ=255)
                                                                                                                                         cAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGc                  >  1:2168907/1‑145 (MQ=255)
                                                                                                                                         cAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGc                  <  2:2168907/145‑1 (MQ=255)
                                                                                                                                         cAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAAtt             <  2:1829435/150‑1 (MQ=255)
                                                                                                                                          aGGCCCCATGGATGGGGACTGCAGCGCGTGACAAGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAAtt             <  2:1831879/150‑1 (MQ=21)
                                                                                                                                           ggCCCCATGGATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACcgcgcg                                              <  2:1703631/115‑1 (MQ=255)
                                                                                                                                           ggCCCCATGGATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACcgcgcg                                              >  1:1703631/1‑115 (MQ=255)
                                                                                                                                           ggCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCg           <  1:363556/150‑1 (MQ=255)
                                                                                                                                           ggCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCg           >  2:363556/1‑150 (MQ=255)
                                                                                                                                           ggCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCg           <  1:1436668/150‑1 (MQ=255)
                                                                                                                                            gCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGaa               >  1:2127615/1‑145 (MQ=255)
                                                                                                                                            gCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGaa               <  2:2127615/145‑1 (MQ=255)
                                                                                                                                            gCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCACGAGATCCCCCGCCGCAATGCCTCCACTCGGCCGCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCgg          <  1:2341772/150‑1 (MQ=21)
                                                                                                                                              cccATCCATGGGGCGTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGAAATGCCTCCCCTCGGCCTCCTGAAGTCGCACTCTGCGTCGCCTGGACTACCCCGCGCCTAGAAGCACGAGCAAAAGCGAGAAGCGAATTCGGTg        >  2:607789/1‑150 (MQ=18)
                                                                                                                                               ccatccatGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGa       >  1:1615034/1‑150 (MQ=255)
                                                                                                                                                 atccatGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCAACTCGGCCGCCTAAAGAAGCTAATGAGATCGGCAGAACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATGCGGTGATc     >  1:1376324/1‑150 (MQ=18)
                                                                                                                                                    catGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATcggc  <  1:2397585/150‑1 (MQ=255)
                                                                                                                                                    catGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATcggc  <  1:1255357/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                       
GAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGC  >  NC_002947/2087229‑2087526

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