Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 354,608 G→A 8.5% intergenic (+16/‑314) prpE → / → codB propionate‑‑CoA ligase/cytosine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913354,6080GA 8.5% 533.6 / 27.6 164intergenic (+16/‑314)prpE/codBpropionate‑‑CoA ligase/cytosine transporter
Reads supporting (aligned to +/- strand):  ref base G (74/76);  new base A (3/11);  total (77/87)
Fisher's exact test for biased strand distribution p-value = 5.32e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

Only 100 of 164 total aligned reads displayed.

GCCAAAAACGCGATCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGT  >  NC_000913/354460‑354755
                                                                                                                                                    |                                                                                                                                                   
gCCAAAAACGCGATCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGt                                                                                                                                                    <  2:1751067/150‑1 (MQ=255)
gCCAAAAACGCGATCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGt                                                                                                                                                    <  2:2644798/150‑1 (MQ=255)
   aaaaaCGCGATCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTcgc                                                                                                                                                 >  1:606751/1‑150 (MQ=255)
   aaaaaCGCGATCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTcgc                                                                                                                                                 >  1:606759/1‑150 (MQ=255)
   aaaaaCGCGATCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTcgc                                                                                                                                                 >  2:382464/1‑150 (MQ=255)
    aaaaCGCGATCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTcgcg                                                                                                                                                >  2:1756232/1‑150 (MQ=255)
    aaaaCGCGATCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTcgcg                                                                                                                                                >  2:1355503/1‑150 (MQ=255)
      aaCGCGATCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGAGGGAAGAGTAGGTTATTGTCGGATGCGTcgcgcg                                                                                                                                              >  2:2227743/1‑150 (MQ=255)
        cgcgATCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGt                                                                                                                                            <  2:410/150‑1 (MQ=255)
        cgcgATCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGt                                                                                                                                            <  1:1439310/150‑1 (MQ=255)
             tCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATc                                                                                                                                       >  1:2171355/1‑150 (MQ=255)
             tCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATc                                                                                                                                       <  1:1296026/150‑1 (MQ=255)
             tCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCACCGCGCGACGCATc                                                                                                                                       <  2:294105/150‑1 (MQ=255)
             tCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCACCGCGCGACGCATc                                                                                                                                       <  2:1728486/150‑1 (MQ=255)
             tCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCACCGCGCGACGCATc                                                                                                                                       <  2:2460997/150‑1 (MQ=255)
             tCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCACCGCGCGACGCATc                                                                                                                                       <  2:436560/150‑1 (MQ=255)
             tCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCACCGCGCGACGCATc                                                                                                                                       <  2:2373495/150‑1 (MQ=255)
             tCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGACGCGTCGCGCGGTGCATc                                                                                                                                       <  1:38706/150‑1 (MQ=255)
             tCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGCCGGATGCACCGCGCGACGCATc                                                                                                                                       <  2:37845/150‑1 (MQ=255)
             tCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGCCGGATGCACCGCGCGACGCATc                                                                                                                                       <  2:915784/150‑1 (MQ=255)
               cGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCg                                                                                                                                     >  1:1173391/1‑150 (MQ=255)
                ggAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCgg                                                                                                                                    <  1:2273490/150‑1 (MQ=255)
                ggAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCgg                                                                                                                                    <  1:1642078/150‑1 (MQ=255)
                ggAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCgg                                                                                                                                    <  1:1642166/150‑1 (MQ=255)
                ggAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCgg                                                                                                                                    <  1:379889/150‑1 (MQ=255)
                ggAAAAATGCTGCGCCGCACAATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCgg                                                                                                                                    <  1:1464402/150‑1 (MQ=255)
                      aTGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACtgt                                                                                                                              <  2:2160766/150‑1 (MQ=255)
                       tgctgcGCCGCACGATCCTGGCGATTGGCGAAGGACGCTATCCGGGAGATCTGACGACCATTGTTGATCCTGCGTGTTTGGATCAGATCCGCCTGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACtgtg                                                                                                                             <  1:1570736/150‑1 (MQ=255)
                       tgctgcGCCGCACGATCCAGGCGTTTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGTTGGAAGAGTAGGTTATTGTCGGTTGCGTCGCGCGGTGCATCCGGCACtgtg                                                                                                                             <  1:1570721/150‑1 (MQ=255)
                       tgctgcGCCGCACGATCCAGGCGATTTGCGAAGGCCGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGTTGGAAGAGTAGGTTATTGTCGGTTGCGTCGCGCGGTGCAACCGGCACtgtg                                                                                                                             <  1:1589951/150‑1 (MQ=255)
                         ctgcGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACtgtgtg                                                                                                                           <  2:639121/150‑1 (MQ=255)
                          tgcGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGc                                                                                                                          <  2:856986/150‑1 (MQ=255)
                                  aCGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCt                                                                                                                  <  2:1177337/150‑1 (MQ=255)
                                           gCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACg                                                                                                         >  2:1475289/1‑150 (MQ=255)
                                           gCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACg                                                                                                         >  2:1285384/1‑150 (MQ=255)
                                                ttGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGAc                                                                                                    <  2:504318/150‑1 (MQ=255)
                                                ttGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGAc                                                                                                    <  2:2113776/150‑1 (MQ=255)
                                                ttGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGAc                                                                                                    <  2:148430/150‑1 (MQ=255)
                                                ttGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGAc                                                                                                    <  2:1321899/150‑1 (MQ=255)
                                                ttGCGAAGGACGCGATCCGGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGACCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGTATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGAc                                                                                                    <  1:844521/150‑1 (MQ=255)
                                                      aGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGCCGCGttt                                                                                              >  2:28008/1‑150 (MQ=255)
                                                      aGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGttt                                                                                              >  2:1654520/1‑150 (MQ=255)
                                                      aGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGttt                                                                                              >  2:2338971/1‑150 (MQ=255)
                                                      aGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGttt                                                                                              >  2:1008943/1‑150 (MQ=255)
                                                      aGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGttt                                                                                              >  2:1902415/1‑150 (MQ=255)
                                                           gcgATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCa                                                                                         >  1:1826185/1‑150 (MQ=255)
                                                           gcgATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCa                                                                                         >  1:1826221/1‑150 (MQ=255)
                                                             gATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATg                                                                                       <  2:1895786/150‑1 (MQ=255)
                                                                ccTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCt                                                                                    <  1:1331639/150‑1 (MQ=255)
                                                                        tCTGACGACCATTTATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGTTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCt                                                                            <  1:1520036/150‑1 (MQ=255)
                                                                        tCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCt                                                                            <  2:58745/150‑1 (MQ=255)
                                                                             cgacCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAAcc                                                                       <  2:1890047/150‑1 (MQ=255)
                                                                             cgacCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAAcc                                                                       <  2:735894/150‑1 (MQ=255)
                                                                             cgacCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAAcc                                                                       <  2:621580/150‑1 (MQ=255)
                                                                              gacCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCg                                                                      >  1:2582789/1‑150 (MQ=255)
                                                                              gacCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCg                                                                      >  1:1276049/1‑150 (MQ=255)
                                                                               acCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTTTTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGt                                                                     <  1:2211582/150‑1 (MQ=255)
                                                                                  aTTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAgg                                                                  <  2:78158/150‑1 (MQ=255)
                                                                                  aTTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGAGGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAgg                                                                  <  2:479203/150‑1 (MQ=255)
                                                                                   ttGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGt                                                                 >  2:1635752/1‑150 (MQ=255)
                                                                                     gatgatCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTTCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCg                                                               <  1:964471/150‑1 (MQ=255)
                                                                                     gatgatCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCg                                                               <  1:2503753/150‑1 (MQ=255)
                                                                                           ccTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAg                                                         <  2:1138447/150‑1 (MQ=255)
                                                                                           ccTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAg                                                         <  2:1766613/150‑1 (MQ=255)
                                                                                                gtcgtTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCt                                                    >  2:1089213/1‑150 (MQ=255)
                                                                                                 tcgtTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTc                                                   >  2:1013111/1‑150 (MQ=255)
                                                                                                 tcgtTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTc                                                   >  2:852570/1‑150 (MQ=255)
                                                                                                 tcgtTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGAGAAGGCGCTc                                                   >  2:2695716/1‑150 (MQ=255)
                                                                                                     tGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGt                                               <  2:1903623/150‑1 (MQ=255)
                                                                                                       gATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCg                                             >  2:1571366/1‑150 (MQ=255)
                                                                                                       gATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCg                                             >  2:2179780/1‑150 (MQ=255)
                                                                                                       gATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCg                                             >  1:1826394/1‑150 (MQ=255)
                                                                                                       gATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCg                                             >  1:1847635/1‑150 (MQ=255)
                                                                                                       gATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCg                                             >  1:769159/1‑150 (MQ=255)
                                                                                                       gATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCg                                             >  2:149455/1‑150 (MQ=255)
                                                                                                       gATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCg                                             >  2:2180062/1‑150 (MQ=255)
                                                                                                         tCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATTCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAATGCGCTCGCGTCGCa                                           >  1:668853/1‑150 (MQ=255)
                                                                                                         tCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCa                                           >  2:2075611/1‑150 (MQ=255)
                                                                                                         tCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCa                                           >  1:2029181/1‑150 (MQ=255)
                                                                                                         tCAGATCCGCCAGGCGATGGAAGAGGAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGAAGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCa                                           >  2:1574992/1‑150 (MQ=255)
                                                                                                          cAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCAt                                          >  2:1898991/1‑150 (MQ=255)
                                                                                                          cAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCAt                                          >  2:31227/1‑150 (MQ=255)
                                                                                                           aGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATc                                         <  2:358726/150‑1 (MQ=255)
                                                                                                           aGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATc                                         <  1:2661621/150‑1 (MQ=255)
                                                                                                           aGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATc                                         >  2:1694156/1‑150 (MQ=255)
                                                                                                           aGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGGCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATc                                         >  2:2504654/1‑150 (MQ=255)
                                                                                                            gATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATcc                                        >  2:1887565/1‑150 (MQ=255)
                                                                                                                 gCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGacac                                   <  2:2582789/150‑1 (MQ=255)
                                                                                                                  ccAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACAcc                                  >  2:260988/1‑150 (MQ=255)
                                                                                                                     ggCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATg                               <  2:1627864/150‑1 (MQ=255)
                                                                                                                     ggCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATg                               <  2:1173391/150‑1 (MQ=255)
                                                                                                                         aTGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCa                           >  1:870973/1‑150 (MQ=255)
                                                                                                                           ggAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGa                         <  1:1974437/150‑1 (MQ=255)
                                                                                                                             aaGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATg                       <  2:2651578/150‑1 (MQ=255)
                                                                                                                                       tATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGa             <  1:79650/150‑1 (MQ=255)
                                                                                                                                           gTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCt         >  2:1224453/1‑150 (MQ=255)
                                                                                                                                            tCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTg        >  2:2547915/1‑150 (MQ=255)
                                                                                                                                            ccGGATGCACCGCGCGACGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTg        >  2:1569089/2‑150 (MQ=255)
                                                                                                                                            ccGGATGCACCGCGCGACGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTg        >  1:1204843/2‑150 (MQ=255)
                                                                                                                                                  gCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGt  >  1:860007/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCCAAAAACGCGATCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGT  >  NC_000913/354460‑354755

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