Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,115,849 C→G 10.4% S16S (TCG→TCC) ‡ wcaM ← putative colanic acid biosynthesis protein WcaM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,115,8490CG10.4% 292.4 / 21.7 96S16S (TCG→TCC) ‡wcaMputative colanic acid biosynthesis protein WcaM
Reads supporting (aligned to +/- strand):  ref base C (38/47);  new base G (6/4);  total (45/51)
Fisher's exact test for biased strand distribution p-value = 5.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

CGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCG  >  NC_000913/2115707‑2115998
                                                                                                                                              |                                                                                                                                                     
cGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCTCGGCGCTTg                                                                                                                                                >  1:494868/1‑150 (MQ=255)
 ggcaggaaCGAAGACTGATTGCCATTCGCTGAAGGCTTGTTTAAATGAGGCGATACAGTCGTGAAGGTTGTACTCGTTAATGTTAAAGATTAGTAGGGCGGGAAGCGCGAGGGCGAAAGGGGTATGGAGGAAGGAAAGCCCACAGCttac                                                                                                                                               >  2:806744/1‑148 (MQ=255)
   caggcaCGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCTCGGCGCTTGcct                                                                                                                                             <  1:1170487/150‑2 (MQ=255)
   caggcaCGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCg                                                                                                                                             <  2:415361/150‑1 (MQ=255)
    aggcaCGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCTCGGCGCTTGcctt                                                                                                                                            <  1:313672/150‑3 (MQ=255)
    aggcaCGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCTCGGCGCTTGcctt                                                                                                                                            <  2:494868/150‑3 (MQ=255)
     ggcaCGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTc                                                                                                                                           >  2:4908/1‑150 (MQ=255)
       caCGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAg                                                                                                                                         <  2:203528/150‑1 (MQ=255)
         cgacgaCTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCTCGGCGCTTGCCTTcctcc                                                                                                                                       <  2:622475/150‑5 (MQ=255)
           acgacTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAAc                                                                                                                                     >  2:1314912/1‑150 (MQ=255)
             gacTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGt                                                                                                                                   >  1:639129/1‑150 (MQ=255)
              acTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTg                                                                                                                                  <  1:741351/150‑1 (MQ=255)
                tGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCg                                                                                                                                <  1:403703/150‑1 (MQ=255)
                 gTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGt                                                                                                                               >  2:325976/1‑150 (MQ=255)
                 gTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGt                                                                                                                               >  1:479659/1‑150 (MQ=255)
                      cccTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGga                                                                                                                          >  2:943347/1‑150 (MQ=255)
                       ccTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGgag                                                                                                                         >  1:1305715/1‑150 (MQ=255)
                          tCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAGGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGt                                                                                                                      >  2:853382/1‑150 (MQ=255)
                            gCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGttt                                                                                                                    >  2:96265/1‑150 (MQ=255)
                             cTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACCCTTTGTCGGGAGGGAAGCGCCCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGtttt                                                                                                                   >  2:307907/1‑150 (MQ=255)
                                  ggCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTaaa                                                                                                              >  2:1250754/1‑150 (MQ=255)
                                   gCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAAt                                                                                                             >  1:410879/1‑150 (MQ=255)
                                   gCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAAt                                                                                                             >  1:416955/1‑150 (MQ=255)
                                   gCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAAt                                                                                                             <  1:461682/150‑1 (MQ=255)
                                    gttGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATg                                                                                                            <  1:635131/149‑1 (MQ=255)
                                    cTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATg                                                                                                            <  2:254060/150‑1 (MQ=255)
                                    cTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATg                                                                                                            <  2:619238/150‑1 (MQ=255)
                                    cTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATg                                                                                                            >  2:489979/1‑150 (MQ=255)
                                    cTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATg                                                                                                            <  1:664080/150‑1 (MQ=255)
                                     ttGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATgg                                                                                                           <  1:71471/150‑1 (MQ=255)
                                          tAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGc                                                                                                      >  1:1532574/1‑150 (MQ=255)
                                            aaTGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGt                                                                                                    >  2:280508/1‑150 (MQ=255)
                                             aTGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCTAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGtt                                                                                                   >  1:1587596/1‑150 (MQ=255)
                                              tGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTtc                                                                                                  >  2:951431/1‑150 (MQ=255)
                                                 ggCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTctcc                                                                                               <  1:1314912/150‑1 (MQ=255)
                                                 ggCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTctcc                                                                                               <  2:491372/150‑1 (MQ=255)
                                                    gATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCctct                                                                                            <  2:491626/150‑1 (MQ=255)
                                                     aTCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCta                                                                                           <  2:1233159/150‑1 (MQ=255)
                                                       ccAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGTGGTATGGATGAAGGCAAGCGCCGAGCTTGCCGTCAGTAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCtata                                                                                         <  2:1524151/150‑1 (MQ=255)
                                                          gTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAg                                                                                      <  2:416955/150‑1 (MQ=255)
                                                          gTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAg                                                                                      <  2:410879/150‑1 (MQ=255)
                                                            cGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGcc                                                                                    <  1:33275/150‑1 (MQ=255)
                                                                tgGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTgca                                                                                <  1:318115/150‑1 (MQ=255)
                                                                 gggTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTgcag                                                                               <  2:639129/150‑1 (MQ=255)
                                                                     tGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAg                                                                           >  2:357260/1‑150 (MQ=255)
                                                                       tAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCt                                                                         >  2:413159/1‑150 (MQ=255)
                                                                         gTCGTTAATGTTCACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTgg                                                                       <  1:1502413/150‑1 (MQ=255)
                                                                         gTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTgg                                                                       <  1:1508109/150‑1 (MQ=255)
                                                                           cGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCg                                                                     <  2:1479184/150‑1 (MQ=255)
                                                                           cGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGAGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCg                                                                     >  2:1287157/1‑150 (MQ=255)
                                                                             ttAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAg                                                                   <  1:1357319/150‑1 (MQ=255)
                                                                                aTGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTtc                                                                <  1:1250754/150‑1 (MQ=255)
                                                                                aTGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTtc                                                                <  2:1606706/150‑1 (MQ=255)
                                                                                 tGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTtct                                                               >  2:1573498/1‑150 (MQ=255)
                                                                                 tGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTtct                                                               >  1:1015838/1‑150 (MQ=255)
                                                                                  gTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTtctc                                                              <  2:301965/150‑1 (MQ=255)
                                                                                  gTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGGTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTtctc                                                              <  2:1119891/150‑1 (MQ=255)
                                                                                  gTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGGCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTtctc                                                              <  1:1570490/150‑1 (MQ=255)
                                                                                   ttAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCg                                                             >  2:784754/1‑150 (MQ=255)
                                                                                      aCGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGAtt                                                          <  1:495568/150‑1 (MQ=255)
                                                                                       cGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGAtt                                                          >  2:412917/1‑149 (MQ=255)
                                                                                       cGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGAtt                                                          <  1:412917/149‑1 (MQ=255)
                                                                                         cTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGAt                                                       <  1:280508/150‑1 (MQ=255)
                                                                                          tttGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATc                                                      <  1:916504/150‑1 (MQ=255)
                                                                                           ttGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCa                                                     <  1:131637/150‑1 (MQ=255)
                                                                                               cgggcgggAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctg                                                 >  1:893359/1‑150 (MQ=255)
                                                                                                 ggcgggAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGGGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctgct                                               <  2:505248/150‑1 (MQ=255)
                                                                                                    gggAAGCGCGCGGGCGAAGGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGt                                            <  1:853382/150‑1 (MQ=255)
                                                                                                     ggAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGtt                                           <  1:1511871/150‑1 (MQ=255)
                                                                                                     ggAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGtt                                           <  2:718925/150‑1 (MQ=255)
                                                                                                        aGCGGGCTGGCGAAAGGGGTATGGAGGATGGCAAGCGCCTAGCTTGCCGTCAGGAACGTGCGTCGGGATAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTaaa                                        <  2:1627615/150‑1 (MQ=255)
                                                                                                          cgcgcgGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAAt                                      >  1:870628/1‑150 (MQ=255)
                                                                                                           gcgcgGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATc                                     >  1:1276472/1‑150 (MQ=255)
                                                                                                              cgGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCgtg                                  <  1:307907/150‑1 (MQ=255)
                                                                                                                ggCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGtt                                >  2:364061/1‑150 (MQ=255)
                                                                                                                ggCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGtt                                <  2:1532574/150‑1 (MQ=255)
                                                                                                                 gCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTc                               <  1:951431/150‑1 (MQ=255)
                                                                                                                  cGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCa                              <  1:96265/150‑1 (MQ=255)
                                                                                                                        gggTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACtttt                        >  1:608428/1‑150 (MQ=255)
                                                                                                                         ggTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACttttt                       <  2:679222/150‑1 (MQ=255)
                                                                                                                         ggTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACttttt                       <  2:1587596/150‑1 (MQ=255)
                                                                                                                         ggTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTACAACttttt                       >  2:1349760/1‑150 (MQ=255)
                                                                                                                              ggaggaAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTcgcgc                  >  2:1237893/1‑150 (MQ=255)
                                                                                                                              ggaggaAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTcgcgc                  >  2:1624844/1‑150 (MQ=255)
                                                                                                                              ggaggaAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCAAGCGTTCTCCTATATAAAGCCTGCAGCAATCTGGCGAGTTCTAGATTTATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTcgcgc                  >  2:435385/1‑150 (MQ=255)
                                                                                                                              ggaggaAGACAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTcgcgc                  >  1:1215379/1‑150 (MQ=255)
                                                                                                                                ctgacggCAAGCTCGGCGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTcgcgcgc                >  1:1116015/6‑150 (MQ=255)
                                                                                                                                ctgacggCAAGCTCGGCGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGCGTTCAACTTTTTcgcgcgc                >  2:441730/6‑150 (MQ=255)
                                                                                                                                 tgacggCAAGCTCGGCGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTcgcgcgcg               >  1:619477/5‑150 (MQ=255)
                                                                                                                                 tgacggCAAGCTCGGCGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTcgcgcgcg               >  2:142974/5‑150 (MQ=255)
                                                                                                                                   aaGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTcgcgcgcgcg             <  2:479659/150‑1 (MQ=255)
                                                                                                                                    cggCAAGCTCGGCGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGt            >  2:701585/2‑150 (MQ=255)
                                                                                                                                    aGGCAAGCTCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGTCATGCGTTCTCCTCTATACAGCCTGCAGCAAGCTGGCGAGTTCTCGATTTATAACATGCTGGTTAAAATCGTCTTAAACTTTTTAGCGCGCGCGt            >  2:1148530/1‑150 (MQ=255)
                                                                                                                                     ggCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGtt           >  1:1516766/1‑150 (MQ=255)
                                                                                                                                      gCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCAAGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGttt          >  2:700791/1‑150 (MQ=255)
                                                                                                                                          gcgcCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTgacg      <  2:1305715/150‑1 (MQ=255)
                                                                                                                                              cGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCg  <  1:489979/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
CGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCG  >  NC_000913/2115707‑2115998

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