Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,051,873 A→G 19.6% intergenic (+40/‑61) quiA → / → oprB‑III quinate dehydrogenase/carbohydrate‑selective porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,051,8730AG19.6% 204.1 / 44.2 103intergenic (+40/‑61)quiA/oprB‑IIIquinate dehydrogenase/carbohydrate‑selective porin
Reads supporting (aligned to +/- strand):  ref base A (37/42);  new base G (10/10);  total (47/56)
Fisher's exact test for biased strand distribution p-value = 8.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.24e-01

Only 100 of 107 total aligned reads displayed.

TGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCA  >  NC_002947/4051724‑4052021
                                                                                                                                                     |                                                                                                                                                    
tGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACaa                                                                                                                                                      <  2:2200923/150‑1 (MQ=255)
tGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACaa                                                                                                                                                      <  1:4184940/150‑1 (MQ=255)
 gcccTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAg                                                                                                                                                     <  2:4120651/148‑1 (MQ=255)
 gACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAg                                                                                                                                                     <  1:741013/150‑1 (MQ=255)
 gACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAg                                                                                                                                                     <  2:2055338/150‑1 (MQ=255)
 gACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAg                                                                                                                                                     <  1:3034014/150‑1 (MQ=255)
    cTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCCAGCCCGCTCCGACAAGGcc                                                                                                                                                  >  2:3884113/1‑150 (MQ=255)
    cTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAACCCCGCTCCGACAAGGcc                                                                                                                                                  >  2:637560/1‑150 (MQ=255)
     tACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCCAGCCCGCACCAACCAGGccc                                                                                                                                                 >  1:2809901/1‑150 (MQ=255)
       cGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCt                                                                                                                                               <  1:4085151/150‑1 (MQ=255)
         tctcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGt                                                                                                                                             <  1:3312785/150‑1 (MQ=255)
         tctcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGt                                                                                                                                             <  2:3318078/150‑1 (MQ=255)
          ctcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTa                                                                                                                                            <  1:4020228/150‑1 (MQ=255)
          ctcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTa                                                                                                                                            >  2:2435064/1‑150 (MQ=255)
          ctcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACAGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCGGTa                                                                                                                                            >  2:5006340/1‑150 (MQ=255)
           tcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCCACCCCGCTCCTACCAGGccccgggtg                                                                                                                                           >  1:144514/1‑145 (MQ=255)
           tcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAg                                                                                                                                           <  2:2212572/150‑1 (MQ=255)
           tcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAg                                                                                                                                           <  1:187297/150‑1 (MQ=255)
           tcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAg                                                                                                                                           <  1:187154/150‑1 (MQ=255)
             gCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCGTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGa                                                                                                                                         <  2:3892465/150‑1 (MQ=255)
               cGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAATCCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGc                                                                                                                                       <  1:4453883/150‑1 (MQ=255)
               cGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGc                                                                                                                                       <  1:3885887/150‑1 (MQ=255)
                gAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCg                                                                                                                                      <  2:752609/150‑1 (MQ=255)
                 aaGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGTCGCCCGCCAGTCGACCGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGTCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCgg                                                                                                                                     <  1:155341/150‑1 (MQ=255)
                 aaGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCgg                                                                                                                                     <  2:4854734/150‑1 (MQ=255)
                 aaGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCgg                                                                                                                                     <  1:1925186/150‑1 (MQ=255)
                  aGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCggg                                                                                                                                    <  2:1039627/150‑1 (MQ=255)
                   gACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAAGCGCGTCGCCTGCTTCGCGGGCTAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGc                                                                                                                                   <  2:107030/150‑1 (MQ=255)
                   gACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGc                                                                                                                                   <  2:1882409/150‑1 (MQ=255)
                       ggCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCACGCCCGCTCCTACAAGGCCCCGGTAGGAGCGGGCGTGc                                                                                                                               >  2:2996104/1‑150 (MQ=255)
                       ggCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGc                                                                                                                               <  1:84413/150‑1 (MQ=255)
                        gCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGcc                                                                                                                              <  1:3464904/150‑1 (MQ=255)
                             cAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCAAGCCcgctc                                                                                                                         <  2:388816/150‑3 (MQ=255)
                                 aCGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGc                                                                                                                     <  1:1614949/150‑1 (MQ=255)
                                   gtggtgATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCACGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGccc                                                                                                                   <  2:1960622/150‑1 (MQ=255)
                                    tggtgATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGAGATTTCTTACGCTCTGCCGAAAAACCCCGTCGCATGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGcccc                                                                                                                  <  2:3921796/150‑1 (MQ=255)
                                        gATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCaa                                                                                                              <  2:4833413/150‑1 (MQ=255)
                                         aTCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGTGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCaaa                                                                                                             <  1:1461020/150‑1 (MQ=255)
                                         aTCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAACCCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCaaa                                                                                                             <  1:2851229/150‑1 (MQ=255)
                                         aTCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCaaa                                                                                                             <  1:766892/150‑1 (MQ=255)
                                         aTCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCaaa                                                                                                             <  1:766885/150‑1 (MQ=255)
                                           cACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGa                                                                                                           <  2:2471741/150‑1 (MQ=255)
                                           cACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGa                                                                                                           <  2:2471800/150‑1 (MQ=255)
                                            aCTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCCAGCCCGCGCCTAAAAGGCCCCGGGAGGTGCGGGCGTGCCCGCGAAGGCCCCGCAAATAc                                                                                                          >  1:2689388/1‑150 (MQ=255)
                                             ctgctgGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCCCGCCGACCACGCCCCGGGTGGAGCGGGCTGGCCCGCGAATGCCCCGCAAAGACa                                                                                                         >  2:1834153/1‑150 (MQ=255)
                                              tgctgGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCAAGCCcgctcc                                                                                                                        >  1:4530133/1‑131 (MQ=255)
                                              tgctgGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCAAGCCcgctcc                                                                                                                        <  2:4530133/134‑4 (MQ=255)
                                              tgctgGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCATAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCCCGACAACAAGGCCCCAGTTTGAGCGGGCTTTCCCCCGAATGCCCCACAAAGaaaa                                                                                                        >  1:967669/1‑150 (MQ=255)
                                                ctgGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCGACCCCGCGCCAACAATGCCCCGGTAGGTGCGGGCGTGCCCGCGCAAGCCCCACAAAGACaaac                                                                                                      >  1:224286/1‑148 (MQ=255)
                                                   gcggcgCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCga                                                                                                   <  2:2502221/150‑1 (MQ=255)
                                                     ggcgCCCGCCTGTCTACTGACCGTGGCGACTTTGTGTTTTCTTACGCTCTGCCTTAAAACCGCGTCGCCTGCTTCGCGGGAACGCCCGCTCCTACACGGCCCCTGTTGGAGGGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCgaga                                                                                                 <  1:637560/150‑1 (MQ=255)
                                                            gccAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTa                                                                                          <  2:3324465/150‑1 (MQ=255)
                                                             ccAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTaa                                                                                         <  2:3408442/150‑1 (MQ=255)
                                                                   gACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCCCGCCAACAAGGCCCCCGTAGGAGCGGGCGTGCCCGCGAAAGCCCCCCAAAGACAACCCGAGATTATTTAAttgcat                                                                                   >  1:3578197/1‑145 (MQ=255)
                                                                            tGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTc                                                                          >  2:762453/1‑150 (MQ=255)
                                                                              gCGACTTTGTGATTTATTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACACGGCCCCTGTAGGTGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGc                                                                        <  1:1834153/150‑1 (MQ=255)
                                                                                   tATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAATCCCGCTCCTACAAGGCCCCGGTAGGAGCGGGCTTGCCCCCGAATTCCCCACAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCAc                                                                   >  2:838254/1‑150 (MQ=255)
                                                                                   tATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCAc                                                                   <  2:3121827/150‑1 (MQ=255)
                                                                                   tATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTAAAAGGCCCCGGGAGGAGCGGGCTTGCCCGCGAAAGCCCCTCAAAGACAACCCGAGAATCTGTAATGCCAACCGCAATTCGCATCAc                                                                   >  2:293114/1‑150 (MQ=39)
                                                                                   tATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCAACTAAAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCACAAAGACAACCCGAGATTCTGTAAAGCCAACCGCAAATAGCAacac                                                                   >  2:1314283/1‑150 (MQ=38)
                                                                                   tATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAACCCCGCTCCAACAACGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCACAAAGACAAACCGCGATTCAGTAAAGCCATCCGCAATTAGCAGCAc                                                                   >  2:195955/1‑150 (MQ=38)
                                                                                        gctttCTTACGCTCTGCCGTACAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTCCAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCa                                                              <  2:5076684/148‑1 (MQ=255)
                                                                                              ttACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTa                                                                                                                                            >  2:4616724/1‑66 (MQ=255)
                                                                                              ttACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTa                                                                                                                                            <  1:4616724/66‑1 (MQ=255)
                                                                                                  gCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCCCTCCTACAAGGCCCCCGTAGGCGCGGGCTTGCCCGCGAAAGCCCCGCAAAGACAACCCCAGATTCAGTAATGCCAACCGCAATTCGCATCACCCCCACCcagctgct                                                    >  2:4285745/1‑150 (MQ=255)
                                                                                                   ctctGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCctgctgctg                                                   <  2:5018842/150‑1 (MQ=255)
                                                                                                    tctGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTCCACGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGctgg                                                  <  1:3509872/150‑1 (MQ=255)
                                                                                                         cGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGATCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGctggccctg                                             <  1:600517/150‑1 (MQ=255)
                                                                                                          gTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGctggccctgg                                            <  2:1169048/150‑1 (MQ=255)
                                                                                                            aaaaCCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGctggccctggcc                                          <  2:3045916/150‑1 (MQ=255)
                                                                                                              aaCCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATTCCATCCGCAAATAGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAg                                        >  2:4484742/1‑150 (MQ=255)
                                                                                                                                   ggCAAGCCCGCTCCGACCAGGCCCCTGGAGGAGCGGGCGTGCCCGCGAATGCCCCTCACCGACAACCCGAGATTCTGTACTGCCATCCGCAATTCGCATCACCCCCACCCGGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcga                   <  2:4126160/150‑1 (MQ=255)
                                                                                                                                      aaGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcgacga                >  1:3139667/1‑150 (MQ=255)
                                                                                                                                        gCCCGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtc              >  1:4373609/1‑150 (MQ=255)
                                                                                                                                        gCCCGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCGTGCCAGCGACGAtc              >  1:352825/1‑150 (MQ=255)
                                                                                                                                         cccGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtct             >  1:1357695/1‑150 (MQ=255)
                                                                                                                                          ccGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCAt                                                                      <  2:4784233/92‑1 (MQ=255)
                                                                                                                                          ccGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtctc            >  1:4622054/1‑150 (MQ=255)
                                                                                                                                          ccGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtctc            >  1:3871377/1‑150 (MQ=255)
                                                                                                                                          ccGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtctc            >  2:2023192/1‑150 (MQ=255)
                                                                                                                                          ccGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtctc            >  1:3365521/1‑150 (MQ=255)
                                                                                                                                          ccGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtctc            >  1:3871383/1‑150 (MQ=255)
                                                                                                                                          ccGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtctc            >  1:3871552/1‑150 (MQ=255)
                                                                                                                                          ccGCTCCGACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCACAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtctc            <  2:4543861/150‑1 (MQ=255)
                                                                                                                                           cGCTCCTCCACGGCCCCTGTAGGAGGGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGATATTCTTTAATGCAATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCa           <  1:852472/150‑1 (MQ=255)
                                                                                                                                            gCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCAt          >  1:2186143/1‑150 (MQ=255)
                                                                                                                                            gCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCAt          >  2:1160075/1‑150 (MQ=255)
                                                                                                                                            gCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCAt          >  1:4425204/1‑150 (MQ=255)
                                                                                                                                            gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGTGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCAt          <  1:813187/150‑1 (MQ=255)
                                                                                                                                            gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCAt          >  1:1224874/1‑150 (MQ=255)
                                                                                                                                            gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCAt          >  1:4420865/1‑150 (MQ=255)
                                                                                                                                            gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCAt          >  1:3693612/1‑150 (MQ=255)
                                                                                                                                             cTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATg         >  2:1168494/1‑150 (MQ=255)
                                                                                                                                              tCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGa        >  1:1404612/1‑150 (MQ=255)
                                                                                                                                              tCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGa        >  2:2596850/1‑150 (MQ=255)
                                                                                                                                               ccTCCAAGCCCCCTGCAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAACACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAc       <  1:4458320/150‑1 (MQ=255)
                                                                                                                                                 tACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAccc     >  1:102149/1‑150 (MQ=255)
                                                                                                                                                 tACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAccc     >  2:873483/1‑150 (MQ=255)
                                                                                                                                                    aaGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCgca  >  1:3973076/1‑150 (MQ=255)
                                                                                                                                                    aaGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCgca  >  1:2229432/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCA  >  NC_002947/4051724‑4052021

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: