Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,976,027:1 +G 6.3% coding (916/1104 nt) pheA → bifunctional chorismate mutase/prephenate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,976,0271.G 6.3% 447.9 / 22.7 111coding (916/1104 nt)pheAbifunctional chorismate mutase/prephenate dehydratase
Reads supporting (aligned to +/- strand):  ref base . (48/56);  new base G (4/3);  total (52/59)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.41e-01

Only 100 of 118 total aligned reads displayed.

CCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTG  >  NC_002947/1975878‑1976173
                                                                                                                                                      |                                                                                                                                                  
 cGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTa                                                                                                                                                   >  2:923253/1‑150 (MQ=255)
  gagaAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCt                                                                                                                                                     >  1:1989562/1‑148 (MQ=255)
  gagaAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCt                                                                                                                                                     <  2:1989562/148‑1 (MQ=255)
  gagaAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTAt                                                                                                                                                  <  2:787931/150‑1 (MQ=255)
  gagaAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTAt                                                                                                                                                  <  2:2210440/150‑1 (MQ=255)
    gaagaTCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCa                                                                                                                                                <  1:798942/150‑1 (MQ=255)
          cGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACgg                                                                                                                                          >  1:2279919/1‑150 (MQ=255)
           gAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGc                                                                                                                                         >  2:1112141/1‑150 (MQ=255)
             aGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTGATAGAACGGCac                                                                                                                                       <  1:2146634/150‑2 (MQ=255)
              gACCTCCCGGACAACGCCACGAGCTTCATGCTGATCGGTAACCAGTAATTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATc                                                                                                                                      <  2:1592341/150‑1 (MQ=255)
              gACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATc                                                                                                                                      <  2:1993549/150‑1 (MQ=255)
                ccGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGa                                                                                                                                    <  1:656479/150‑1 (MQ=255)
                   cccGGACAGCTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCt                                                                                                                                 <  2:1308726/150‑1 (MQ=255)
                   cccGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCt                                                                                                                                 <  2:569891/150‑1 (MQ=255)
                     cGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGa                                                                                                                               <  1:29048/150‑1 (MQ=255)
                        aCAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACGGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCt                                                                                                                                 <  1:2289509/145‑1 (MQ=255)
                        aCAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATc                                                                                                                                      >  1:1660836/1‑140 (MQ=255)
                            cTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCAt                                                                                                                        >  1:2430437/1‑150 (MQ=255)
                             tCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGCCCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATc                                                                                                                       >  1:1035673/1‑150 (MQ=255)
                              ccACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCg                                                                                                                      >  2:1575700/1‑150 (MQ=255)
                               cACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCga                                                                                                                     >  1:1060267/1‑150 (MQ=255)
                               cACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGAACTGACCCGCATCga                                                                                                                     <  2:2279919/150‑1 (MQ=255)
                                 cgcgCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCgaga                                                                                                                   >  1:1792575/1‑150 (MQ=255)
                                  gcgcTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGAc                                                                                                                  <  2:1060267/150‑1 (MQ=255)
                                      ttCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGAcccgc                                                                                                              <  1:1406418/150‑1 (MQ=255)
                                        ccTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGAcccgccc                                                                                                            <  2:1114033/150‑1 (MQ=255)
                                          tgatgaTCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGt                                                                                                          >  1:2218692/1‑150 (MQ=255)
                                          tgatgaTCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGt                                                                                                          >  2:931287/1‑150 (MQ=255)
                                           gatgatCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTc                                                                                                         >  1:199048/1‑150 (MQ=255)
                                             tgatCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGATGCCGTTCTGATAGAACGGCacca                                                                                                                                     <  2:1001097/120‑4 (MQ=255)
                                             tgatCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTGGCATGAGTTGCTGATGCCGTTCTGATAGAACGGCacca                                                                                                                                     >  1:1001097/1‑117 (MQ=255)
                                             tgatCGGTAACAAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTcgc                                                                                                       <  1:937198/150‑1 (MQ=255)
                                               atCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCgcgc                                                                                                     >  1:2142113/1‑150 (MQ=255)
                                                 cGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAg                                                                                                   <  1:923253/150‑1 (MQ=255)
                                                  ggTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGc                                                                                                  <  1:2359363/150‑1 (MQ=255)
                                                  ggTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGc                                                                                                  <  1:2080148/150‑1 (MQ=255)
                                                   gTAATCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTCCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCg                                                                                                 <  2:1192205/150‑1 (MQ=255)
                                                       ccAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCaa                                                                                             <  1:392372/150‑1 (MQ=255)
                                                        cAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTGATAGAACGGCa                                                                                                                                        <  2:757874/106‑1 (MQ=255)
                                                        cAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTGATAGAACGGCa                                                                                                                                        >  1:757874/1‑106 (MQ=255)
                                                        cAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAg                                                                                            >  2:1291572/1‑150 (MQ=255)
                                                        cAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAg                                                                                            >  2:2280227/1‑150 (MQ=255)
                                                          ggAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTg                                                                                          <  2:1091658/150‑1 (MQ=255)
                                                              gTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGAcc                                                                                      <  1:1112141/150‑1 (MQ=255)
                                                               tACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCAt                                                                                                                        <  2:490975/115‑1 (MQ=255)
                                                               tACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCAt                                                                                                                        >  1:490975/1‑115 (MQ=255)
                                                               tACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCt                                                                                     >  2:2027822/1‑150 (MQ=255)
                                                               tACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCt                                                                                     >  2:980075/1‑150 (MQ=255)
                                                                 ccgccgACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTAc                                                                                   <  1:1193360/150‑1 (MQ=255)
                                                                   gccgACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTcgc                                                                                                       >  2:1811940/1‑128 (MQ=255)
                                                                   gccgACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTcgc                                                                                                       <  1:1811940/128‑1 (MQ=255)
                                                                   gccgACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACgt                                                                                 <  1:1291572/150‑1 (MQ=255)
                                                                   gccgACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGATCCTAGCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACgt                                                                                 <  2:1796522/150‑1 (MQ=255)
                                                                    ccgACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACgtg                                                                                <  1:1110150/150‑1 (MQ=255)
                                                                    ccgACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACgtg                                                                                >  2:767396/1‑150 (MQ=255)
                                                                         cGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGAccc                                                                                                                >  1:609084/1‑113 (MQ=255)
                                                                         cGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGAccc                                                                                                                <  2:609084/113‑1 (MQ=255)
                                                                           gCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCa                                                                         <  2:199048/150‑1 (MQ=255)
                                                                            cgacgaCAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCAt                                                                        <  1:95761/150‑1 (MQ=255)
                                                                             gacgacAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATc                                                                       <  1:980075/150‑1 (MQ=255)
                                                                              acgacAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCg                                                                      >  1:1070387/1‑150 (MQ=255)
                                                                                 acAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGa                                                                                        >  2:1498469/1‑129 (MQ=255)
                                                                                 acAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGa                                                                                        <  1:1498469/129‑1 (MQ=255)
                                                                                 acAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACt                                                                   >  1:1951722/1‑150 (MQ=255)
                                                                                        cTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCgg                                                            >  1:1118204/1‑150 (MQ=255)
                                                                                               atcGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCcaccac                                                     >  2:2330666/1‑150 (MQ=255)
                                                                                                tcGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGccacca                                                      <  2:2153222/148‑1 (MQ=255)
                                                                                                tcGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGccacca                                                      >  1:2153222/1‑148 (MQ=255)
                                                                                                  gTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACcgc                                                  >  1:1429517/1‑150 (MQ=255)
                                                                                                   tttCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTAc                                                                                   <  2:911257/116‑1 (MQ=255)
                                                                                                   tttCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTAc                                                                                   >  1:911257/1‑116 (MQ=255)
                                                                                                         tGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGc                                           <  2:1777114/150‑1 (MQ=255)
                                                                                                          gAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCGACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCt                                          <  2:1035673/150‑1 (MQ=255)
                                                                                                            gCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGa                                        >  1:1389068/1‑150 (MQ=255)
                                                                                                             caacaaGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTTAt                                       >  2:307631/1‑150 (MQ=255)
                                                                                                             caacaaGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGAt                                       <  2:2430437/150‑1 (MQ=255)
                                                                                                              aacaaGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATc                                      >  1:590001/1‑150 (MQ=255)
                                                                                                              aacaaGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATc                                      <  2:1747053/150‑1 (MQ=255)
                                                                                                                  aGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAgg                                  <  2:2280822/150‑1 (MQ=255)
                                                                                                                  aGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAgg                                  <  2:308270/150‑1 (MQ=255)
                                                                                                                   gCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGc                                 <  2:2447492/150‑1 (MQ=255)
                                                                                                                    ccAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCg                                <  1:2330666/150‑1 (MQ=255)
                                                                                                                     cAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCgg                               <  1:1575700/150‑1 (MQ=255)
                                                                                                                      aGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGt                              >  1:1595628/1‑150 (MQ=255)
                                                                                                                       ggTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGt                              <  1:906765/149‑1 (MQ=255)
                                                                                                                       ggTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGt                              >  2:906765/1‑149 (MQ=255)
                                                                                                                            gTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGa                        <  2:590001/150‑1 (MQ=255)
                                                                                                                             ttGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAg                       <  2:1792575/150‑1 (MQ=255)
                                                                                                                               gCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCg                                                                      >  1:94039/1‑101 (MQ=255)
                                                                                                                               gCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCg                                                                      <  2:94039/101‑1 (MQ=255)
                                                                                                                                   gAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGt                                                                                           >  1:2437913/1‑76 (MQ=255)
                                                                                                                                   gAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGt                                                                                           <  2:2437913/76‑1 (MQ=255)
                                                                                                                                        gCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCa            <  2:1951722/150‑1 (MQ=255)
                                                                                                                                            gTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCaggag        <  1:753125/150‑1 (MQ=255)
                                                                                                                                            atGCCGTTCTGATAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTAGATTTTCTTCATCGCCTTCGTCGGCCACAAGACTAGCCCGGTGGTCAAGGCGGTGCTCGAGCAGATCAGCCaggag        >  1:1116436/2‑150 (MQ=255)
                                                                                                                                              gCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGccac                                                        >  1:941248/1‑100 (MQ=255)
                                                                                                                                              gCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGccac                                                        <  2:941248/100‑1 (MQ=255)
                                                                                                                                                 gTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCACGGCGGGGCTCGAGCAGCTCCGCCACGGGGCggt   >  2:822439/1‑150 (MQ=255)
                                                                                                                                                  ttCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCtgcg  >  2:2057064/1‑148 (MQ=255)
                                                                                                                                                     tATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTtc                                                                 <  2:698720/84‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
CCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTG  >  NC_002947/1975878‑1976173

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