Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,118,533 A→C 5.1% intergenic (‑129/+147) wcaK ← / ← wzxC colanic acid biosynthesis protein/putative colanic acid exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,118,5330AC 5.1% 288.4 / 4.4 97intergenic (‑129/+147)wcaK/wzxCcolanic acid biosynthesis protein/putative colanic acid exporter
Reads supporting (aligned to +/- strand):  ref base A (43/49);  new base C (2/3);  total (45/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

CAGAATAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGTTGTCGGAACCGAAAACAGCAAC  >  NC_000913/2118387‑2118679
                                                                                                                                                  |                                                                                                                                                  
cAGAATAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGc                                                                                                                                                 >  1:998436/1‑150 (MQ=255)
cAGAATAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGc                                                                                                                                                 >  2:1712876/1‑150 (MQ=255)
cAGAATAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGc                                                                                                                                                 >  2:1001002/1‑150 (MQ=255)
 aGAATAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCg                                                                                                                                                <  1:1428784/150‑1 (MQ=255)
        gTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCAcgc                                                                                                                                         <  2:208121/150‑1 (MQ=255)
        gTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCAcgc                                                                                                                                         <  2:82482/150‑1 (MQ=255)
        gTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCAcgc                                                                                                                                         <  2:1225781/150‑1 (MQ=255)
         tAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCAcgcc                                                                                                                                        >  2:567936/1‑150 (MQ=255)
         tAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCAcgcc                                                                                                                                        <  1:1457166/150‑1 (MQ=255)
           aTTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATGCGGCGTGAAcgcctt                                                                                                                                      >  1:296119/1‑148 (MQ=255)
               cATTTTTTCCTCCTAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATcc                                                                                                                                  <  2:350017/150‑1 (MQ=255)
                 ttttttCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                >  2:629691/1‑150 (MQ=255)
                  tttttCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATGCGGCGTGAACGCCTTATCCgtc                                                                                                                               <  2:296119/150‑3 (MQ=25)
                  tttttCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                               <  2:431605/150‑1 (MQ=255)
                  tttttCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                               <  2:224383/150‑1 (MQ=255)
                    tttCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATGCGGCGTGAACGCCTTATCCgtcct                                                                                                                             <  1:997188/150‑5 (MQ=25)
                     ttCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAtt                                                                                                                            <  2:1276854/150‑1 (MQ=255)
                       ccTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATGCGGCGtgaacg                                                                                                                                          >  2:258157/1‑129 (MQ=255)
                       ccTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATGCGGCGtgaacg                                                                                                                                          <  1:258157/134‑6 (MQ=255)
                        cTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAg                                                                                                                         >  1:1152749/1‑150 (MQ=255)
                        cTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAg                                                                                                                         >  1:527482/1‑150 (MQ=255)
                        cTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAg                                                                                                                         <  1:865095/150‑1 (MQ=255)
                          cATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTg                                                                                                                       >  1:829057/1‑150 (MQ=255)
                               aTTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGa                                                                                                                  >  1:1242498/1‑150 (MQ=255)
                                  tGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGc                                                                                                               <  2:499690/150‑1 (MQ=255)
                                  tGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGc                                                                                                               <  2:1249166/150‑1 (MQ=255)
                                     tGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGAc                                                                                                            >  1:371551/1‑150 (MQ=255)
                                      gCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACg                                                                                                           >  2:212769/1‑150 (MQ=255)
                                            tGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGAcg                                                                                                     <  1:1576053/150‑1 (MQ=255)
                                               ggCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGt                                                                                                  >  2:560587/1‑150 (MQ=255)
                                               ggCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCCGACGCGTCTTATCAGGCCTACAAGACCCGTGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGt                                                                                                  >  2:1476291/1‑150 (MQ=255)
                                                gCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACTGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTc                                                                                                 <  2:1172628/150‑1 (MQ=255)
                                                gCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTc                                                                                                 <  1:1104036/150‑1 (MQ=255)
                                                 cTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCt                                                                                                >  2:568060/1‑150 (MQ=255)
                                                  tGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCtt                                                                                               <  1:702408/150‑1 (MQ=255)
                                                         gttttgttCGTTGCCGCTTTTTGCCTGCTGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATAAGACGCGTCAGCGTCGCATCAGGc                                                                                        <  1:1001002/150‑1 (MQ=21)
                                                          ttttgttCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGcc                                                                                       <  1:629691/150‑1 (MQ=255)
                                                          ttttgttCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGcc                                                                                       <  2:1389952/150‑1 (MQ=255)
                                                          ttttgttCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGcc                                                                                       <  1:731464/150‑1 (MQ=255)
                                                          ttttgttCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGcc                                                                                       <  2:1593092/150‑1 (MQ=255)
                                                          tttggttCGTTGCCGCGTTTTGCCTGATGCGACGCTGACGCGTCTTAGCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGcc                                                                                       <  1:271179/150‑1 (MQ=255)
                                                            ttgttCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTa                                                                                     >  1:1026439/1‑150 (MQ=255)
                                                               ttCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAg                                                                                  <  1:633721/150‑1 (MQ=255)
                                                               ttCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAg                                                                                  <  2:829057/150‑1 (MQ=255)
                                                                tCGTTGCCGCTATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAgg                                                                                 <  1:47927/150‑1 (MQ=255)
                                                                   ttGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTcc                                                                              <  1:1665461/150‑1 (MQ=255)
                                                                    tGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTccc                                                                             >  2:593926/1‑150 (MQ=255)
                                                                      ccGCTTTTTGCCTGATGCGACGCCGACGCGTCTTATCAGGCCTACAAGACCCGTGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGa                                                                           <  1:1476291/150‑1 (MQ=255)
                                                                       cGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAg                                                                          >  1:358205/1‑150 (MQ=255)
                                                                       cGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACACGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATGCGGCATTCAGTGCCTGATGCGACGCTGACGCGGGTTATCAGGCCGACAGGTCCagag                                                                          >  1:262038/1‑150 (MQ=25)
                                                                         cTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGca                                                                        <  1:1166734/150‑1 (MQ=255)
                                                                          tttttGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGcac                                                                       <  1:212769/150‑1 (MQ=255)
                                                                           ttttGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGcaca                                                                      >  2:1561462/1‑150 (MQ=255)
                                                                               gCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAAc                                                                  >  2:1639409/1‑150 (MQ=255)
                                                                                  tGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGGATTCAGTGCCTGATGCGACGCTGACGCGTCGTATCAGGCCTACAGGTCCAGAGCACAGAACCGt                                                               >  2:1663319/1‑150 (MQ=255)
                                                                                  tGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGc                                                                                                               >  1:1346444/1‑102 (MQ=255)
                                                                                  tGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGc                                                                                                               <  2:1346444/102‑1 (MQ=255)
                                                                                  tGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGt                                                               >  1:783115/1‑150 (MQ=255)
                                                                                        gACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACACGTCCCGAGCACAGAACCGTAGGACg                                                         >  1:1659371/1‑150 (MQ=255)
                                                                                         aCGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACgg                                                        <  2:527482/150‑1 (MQ=255)
                                                                                          cGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGa                                                       >  2:1717045/1‑150 (MQ=255)
                                                                                           gCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTGCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGAt                                                      <  2:70670/150‑1 (MQ=255)
                                                                                            cTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATa                                                     <  1:568060/150‑1 (MQ=255)
                                                                                             tGACTCGTCTTATCAGGCCTACAAGACCCGAGCACAGACCCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATaa                                                    <  1:346276/150‑1 (MQ=255)
                                                                                              gACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAg                                                   <  1:1604458/150‑1 (MQ=255)
                                                                                              gACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAg                                                   >  1:64187/1‑150 (MQ=255)
                                                                                               aCGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAgg                                                  <  2:1026439/150‑1 (MQ=255)
                                                                                                      ttATGAGGTCTTCAAGACCCTAGCACAGAAACGTTGGAGGGATAAGGCGTTCACGCCGCATCGGGCATTGAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTa                                           <  1:1529787/150‑1 (MQ=21)
                                                                                                      ttATCAGGCCTACAAGACCCGAGCACAGACCCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGttttt                                           <  2:1256165/150‑2 (MQ=255)
                                                                                                      ttATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTAACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCACGAGCACAGAAACGTAGGACGGATACGGCGTTTTa                                           >  1:669526/1‑150 (MQ=255)
                                                                                                          aaGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTAcgcc                                       <  2:719484/149‑1 (MQ=255)
                                                                                                           aGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTAcgccg                                      <  2:562408/150‑1 (MQ=255)
                                                                                                            ggCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTAcgccgc                                     >  2:1245808/1‑150 (MQ=255)
                                                                                                             gCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCa                                    >  2:1309796/1‑150 (MQ=255)
                                                                                                              ccTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCaa                                   >  1:1002172/1‑149 (MQ=255)
                                                                                                              ccTACAACACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACCCCGCATCCGGCATTCAGTGCCTCAAGCGACGCTTACGCGTCTTATCAGGCCTACAGCTCACGAGCACACAACCGTAGCACCGATAACGCGTTTTACGCCGCaa                                   >  2:430976/1‑149 (MQ=25)
                                                                                                                  cAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCgg                               <  1:956717/150‑1 (MQ=255)
                                                                                                                    aGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCa                             <  2:218848/150‑1 (MQ=255)
                                                                                                                         caaGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGGAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCgg                        >  1:845111/3‑149 (MQ=255)
                                                                                                                         cGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGt                        <  1:584710/150‑1 (MQ=255)
                                                                                                                              acaGAACCGTAGGACGGATAAGGCGTTCATGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGTTGTCg                   <  2:477370/150‑1 (MQ=255)
                                                                                                                               caGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGCCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGTTGTCgg                  >  2:799102/1‑150 (MQ=255)
                                                                                                                               caGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGACTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGTTGTCgg                  >  2:799100/1‑150 (MQ=255)
                                                                                                                                 gAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGTTGTCGGaa                >  1:1071348/1‑150 (MQ=255)
                                                                                                                                  aaCCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCACTGCCTGATCCGACGCTGACGCGTCTTATCAGGCCTACAGGACCCGAGCACAGAACCGCAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGGTGTCGGAAc               <  1:1058399/150‑1 (MQ=25)
                                                                                                                                    ccGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGTTGTCGGAACCg             >  1:23581/1‑150 (MQ=255)
                                                                                                                                      gTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGTTGTCGGAACCGaa           >  2:1099821/1‑150 (MQ=255)
                                                                                                                                         ggACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGTTGTCGGAACCGAAAAc        >  1:351289/1‑150 (MQ=255)
                                                                                                                                         ggACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGTTGTCGGAACCGAAAAc        >  1:1092433/1‑150 (MQ=255)
                                                                                                                                           aCGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGTTGTCGGAACCGAAAACAg      <  1:1309796/150‑1 (MQ=255)
                                                                                                                                            cGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGTTGTCGGAACCGAAAACAGc     >  1:221054/1‑150 (MQ=255)
                                                                                                                                            cGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGTTGTCGGAACCGAAAACAGc     <  1:1401041/150‑1 (MQ=255)
                                                                                                                                            cGGATAAGGCGTTCACGCCGCAGCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGTTGTCGGAACCGAAAACAGc     >  1:850209/1‑150 (MQ=255)
                                                                                                                                              gATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGTTGTCGGAACCGAAAACAGCaa   >  1:1368254/1‑150 (MQ=255)
                                                                                                                                              gATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGTTGTCGGAACCGAAAACAGCaa   <  2:468740/150‑1 (MQ=255)
                                                                                                                                               aTAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGTTGTCGGAACCGAAAACAGCAAc  >  2:1304163/1‑150 (MQ=255)
                                                                                                                                               aTAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGTTGTCGGAACCGAAAACAGCAAc  <  2:1659371/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CAGAATAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCAACCGTTGTCGGAACCGAAAACAGCAAC  >  NC_000913/2118387‑2118679

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: