Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 404,113 A→C 6.4% A137A (GCA→GCC yaiC → diguanylate cyclase, cellulose regualtor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913404,1130AC 6.4% 300.0 / 6.6 108A137A (GCA→GCCyaiCdiguanylate cyclase, cellulose regualtor
Reads supporting (aligned to +/- strand):  ref base A (56/42);  new base C (2/6);  total (58/50)
Fisher's exact test for biased strand distribution p-value = 1.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

CGAGCAGGGCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTCTGC  >  NC_000913/403965‑404261
                                                                                                                                                    |                                                                                                                                                    
cGAGCAGGGCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGGCGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAg                                                                                                                                                     >  1:1406149/1‑150 (MQ=255)
cGAGCAGGGCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGAAAACGTGATGCCTTCCACCGCTATGTTGATGATTATGTGTCTGAATTTGATCGGGGCAg                                                                                                                                                     >  1:75040/1‑150 (MQ=255)
      ggggCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGcc                                                                                                                                               >  1:1824514/1‑150 (MQ=255)
       ggCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGccc                                                                                                                                              <  1:422118/150‑1 (MQ=255)
       ggCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGccc                                                                                                                                              >  1:30389/1‑150 (MQ=255)
       ggCCGTCGATCCGCTTAGCAGGTAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCTCGATGTTGTTGATTATGTGTCTGAATTTGATGGGGGCAGGCGGccc                                                                                                                                              <  1:188370/150‑1 (MQ=255)
        gCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGTGGCAGGCGGcccc                                                                                                                                             >  1:86988/1‑150 (MQ=255)
        gCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGcccc                                                                                                                                             >  1:1889486/1‑150 (MQ=255)
          cGTCGATCCGCTTATCCTGGAAATTTACAACATAAAAACCGATGCAGTATTAGCGGGTATGTTGGTAGTCGTAATGTGCGTACACGTGCTGCCTTCCACCGGGATGTTGATGATTATGTGTCTGAATTTGATGGGGTCAGGCGGCCCCCg                                                                                                                                           <  2:1018937/150‑1 (MQ=255)
          cGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCg                                                                                                                                           >  1:1849768/1‑150 (MQ=255)
            tCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCCGCCTGCCCCCATc                                                                                                                                         <  1:694636/150‑1 (MQ=255)
             cGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCt                                                                                                                                        >  1:7549/1‑150 (MQ=255)
             cGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCt                                                                                                                                        >  2:1290155/1‑150 (MQ=255)
              gATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTg                                                                                                                                       <  2:983500/150‑1 (MQ=255)
               aTCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGt                                                                                                                                      >  2:1466327/1‑150 (MQ=255)
                tCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGCTGCAGTATTAGCGGGAATGTGGGCAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGtt                                                                                                                                     <  1:1010203/150‑1 (MQ=255)
                tCCGATTAGCCGGGATATTTACAACTTGAAAACCGATGCAGTATTAGCGGGAAGGAGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTACACCGCGAGGTTGATGATTATGTGTCTGAATTTGATGGGGGCTGGCGGCCCCCGTCTGtt                                                                                                                                     <  2:731160/150‑1 (MQ=255)
                      tAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTcgg                                                                                                                               >  1:437804/1‑149 (MQ=255)
                        gCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCgg                                                                                                                             >  2:854099/1‑150 (MQ=255)
                        gCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCgg                                                                                                                             >  2:605553/1‑150 (MQ=255)
                        gCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCgg                                                                                                                             >  2:814463/1‑150 (MQ=255)
                        gCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCACGTCTGTTTGTCGCgg                                                                                                                             >  1:1703663/1‑150 (MQ=255)
                         ccGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCggg                                                                                                                            >  2:560848/1‑150 (MQ=255)
                         ccGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCggg                                                                                                                            >  1:887138/1‑150 (MQ=255)
                          cGGGAAATTTACAACTTAAAACCCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGt                                                                                                                           <  1:1312198/150‑1 (MQ=255)
                            ggAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCt                                                                                                                         <  2:1424677/150‑1 (MQ=255)
                              aaaTTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGGGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTgg                                                                                                                       <  2:1499462/150‑1 (MQ=255)
                               aaTTTACACCTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGgt                                                                                                                      >  1:1616551/1‑150 (MQ=255)
                                   tacttcTTATAAAGCGATGCAGTATTAGCGGGATGGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTTCACCGCGATGTTGATGATTATGTGTCTGTATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTg                                                                                                                  <  1:42780/145‑1 (MQ=255)
                                   tacctcTTAAAAACCGATGCACTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTg                                                                                                                  <  1:1789354/145‑1 (MQ=255)
                                   tACATCTTAAAAACCGATGCAGTATTAGCGGTAATGTGGGTAGGCGTAATGTGCGTAATCGTGCTGCCTTCCACGGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGGCGCTGGTCTGGTGTTg                                                                                                                  <  2:1330141/150‑1 (MQ=255)
                                     cAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGAt                                                                                                                >  1:961108/1‑150 (MQ=255)
                                     cAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGAt                                                                                                                <  2:1493714/150‑1 (MQ=255)
                                     cAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGAt                                                                                                                >  1:1013036/1‑150 (MQ=255)
                                     cAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGAt                                                                                                                >  1:1306656/1‑150 (MQ=255)
                                       aCTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGAtgg                                                                                                              >  1:1442429/1‑150 (MQ=255)
                                       aCTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGAtgg                                                                                                              <  2:1446641/150‑1 (MQ=255)
                                         tttaaacCCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCATCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGAtggtg                                                                                                            >  1:220446/4‑150 (MQ=255)
                                          tAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGAtggtgg                                                                                                           <  1:878023/150‑1 (MQ=255)
                                              aaCCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTc                                                                                                       >  1:654463/1‑150 (MQ=255)
                                               aCCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTcc                                                                                                      <  2:994860/150‑1 (MQ=255)
                                                 cGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTCGTTTCCTg                                                                                                    <  2:1778399/150‑1 (MQ=255)
                                                    tGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGGGGTTTCCTGCCt                                                                                                 <  2:975452/150‑1 (MQ=255)
                                                     gCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCtt                                                                                                <  2:1889486/150‑1 (MQ=255)
                                                      cAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTg                                                                                               >  2:1593179/1‑150 (MQ=255)
                                                      cAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTg                                                                                               <  1:927533/150‑1 (MQ=255)
                                                          aTTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCAc                                                                                           >  2:8776/1‑150 (MQ=255)
                                                          aTTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCAc                                                                                           >  2:439843/1‑150 (MQ=255)
                                                            tAGCGGGAATGGGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCAccc                                                                                         <  1:324592/150‑1 (MQ=255)
                                                               cGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCg                                                                                      <  2:866027/150‑1 (MQ=255)
                                                                 ggAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAg                                                                                    <  1:1884228/150‑1 (MQ=255)
                                                                       tgGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACg                                                                              >  2:14887/1‑150 (MQ=255)
                                                                        gggTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTAGGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCCCCCTCGAGCTGACgg                                                                             >  2:132386/1‑150 (MQ=255)
                                                                         ggTAGGCGTCATGGGCGTAACCGTGCTGCCTTCCACCGCGATGTTGATGATTAGGTGGCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACggg                                                                            <  1:47137/150‑1 (MQ=255)
                                                                          gttggCGTAATGGGCGTATACGTGCTGCCTTCCTCCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTATGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGc                                                                           <  2:774632/147‑1 (MQ=255)
                                                                          gTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGc                                                                           >  2:1308056/1‑150 (MQ=255)
                                                                           tAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCa                                                                          <  2:415965/150‑1 (MQ=255)
                                                                            aGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCAt                                                                         >  2:1379413/1‑150 (MQ=255)
                                                                                     gggCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTAGGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTc                                                                <  1:1310351/150‑1 (MQ=255)
                                                                                         gTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTc                                                            <  1:598472/150‑1 (MQ=255)
                                                                                           aaaCGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCaa                                                          >  2:1739392/1‑150 (MQ=255)
                                                                                            accgggctgcCTTCCACCGGGTTGTTGATGATTATGTGTCTGACTTTGATGGGGGCCGGCGGCCCCCGTCTGTTTGGCGCGGGTCTGGTGTTGATGGTTGTTTCCTGCCTTGTCTCCATCGAGCTGACGGGCATTACCGTGTCGTTCAAt                                                         <  2:1183127/145‑1 (MQ=255)
                                                                                            accGTGCTGCCGTCCACCGCGAGTGTGATGATTCTGTGTCTGAATTTGATGGGGGCTGGCGGCCCCCGTCTGTTTGTCGCGGGGCTGGTGTTGAAGGTGGTTTCATGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAAt                                                         <  1:115640/148‑1 (MQ=255)
                                                                                            aaCGGGCTGCCTTCCCCCGCGCTGGGGCGGACTCTCTGGCCGCCTATGATGGGCGCCGCCGGCCCCCGTCTGTTTGCCGCGGGTCGGGGGTTGATGGTGGGTTCCTGCCTTGTCACCCTCGAGCGGACGGGCATTACCGTGTCGTTCAAt                                                         <  2:1241530/150‑1 (MQ=255)
                                                                                                tgctgcCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTg                                                     >  2:22653/1‑150 (MQ=255)
                                                                                                ggctgcGTTGCACCGCGATGTGGATTATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGGGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTg                                                     <  1:377272/149‑1 (MQ=255)
                                                                                                ggctgcCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTg                                                     <  1:1393729/149‑1 (MQ=255)
                                                                                                  ctgcCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTgcg                                                   >  2:982718/1‑150 (MQ=255)
                                                                                                    gcCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGcgcc                                                 >  1:1118270/1‑150 (MQ=255)
                                                                                                    gcCTTCCACCGATATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGATGGTCTGGTGTTTATTTTGTTTTCCTTCCTTGTCACCCTAGAGCTGACGGGCATTACAGTTTCTTTCAATAGTGctcc                                                 >  1:235070/1‑150 (MQ=255)
                                                                                                      cTTCCACGGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTCCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGcgccgc                                               >  1:40061/1‑150 (MQ=255)
                                                                                                            ccGCGATGTTGATGATTATGTGTCTGAATTTGTTGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAAt                                         >  2:353672/1‑150 (MQ=255)
                                                                                                                 aTGTTGCTGATTATGTGTCTGAATTTGATGGGGGCAGTCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGGTCAATAGTGCGCCGCTGGAAtggtgg                                    <  2:1019945/150‑1 (MQ=255)
                                                                                                                  tGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGc                                   <  1:8776/150‑1 (MQ=255)
                                                                                                                  tGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGc                                   <  2:30389/150‑1 (MQ=255)
                                                                                                                  tGTTGATGATTATGTGTCGGAATTTGATGGGGGCAGGCGGCCCCCGTCGGTTTGTCGCGGGTCTGGTGTTGATGGTTGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGc                                   <  2:1406812/150‑1 (MQ=255)
                                                                                                                     tgatgaTTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTAACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGctct                                >  1:531136/1‑150 (MQ=255)
                                                                                                                      gatgatTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCtctc                               >  2:559462/1‑150 (MQ=255)
                                                                                                                      gatgatTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCtctc                               >  1:939093/1‑150 (MQ=255)
                                                                                                                       atgatTATGTGTCTGAAGGTGATGGGGGCAGGCGGCCCCCGGCTGTTTGTCGCGGTTCTGGTGTTGATGGTGGTTTCCTGCCTTGGCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTcc                              <  2:303693/150‑1 (MQ=255)
                                                                                                                            tATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTcc                         <  2:990717/150‑1 (MQ=255)
                                                                                                                            aaTGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTcc                         <  2:454759/149‑1 (MQ=255)
                                                                                                                                tgtCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTTCCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCatt                     >  1:739402/1‑150 (MQ=255)
                                                                                                                                 gtCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCattt                    >  2:439129/1‑149 (MQ=255)
                                                                                                                                 gtCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCattt                    >  1:1487395/1‑149 (MQ=255)
                                                                                                                                 gtCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCatta                    >  2:1812387/1‑150 (MQ=255)
                                                                                                                                 gtCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCatta                    >  1:898520/1‑150 (MQ=255)
                                                                                                                                  tCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCattat                   >  2:1544432/1‑150 (MQ=255)
                                                                                                                                  tCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGAGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGCCGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCagtat                   >  1:1445761/1‑150 (MQ=255)
                                                                                                                                    tGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCGATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTg                 >  2:656010/1‑150 (MQ=255)
                                                                                                                                     gAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCTCCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGt                >  1:1897070/1‑150 (MQ=255)
                                                                                                                                       aTTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGTAGGTTTCCTCCCTTTTCACCCTCTAGCTAACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCa              >  2:2784/1‑150 (MQ=255)
                                                                                                                                       aTTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCa              >  1:1314290/1‑150 (MQ=255)
                                                                                                                                           gacgggggCCGCCTGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCAtttt          >  2:1575257/4‑149 (MQ=255)
                                                                                                                                           gATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTa          <  1:494369/150‑1 (MQ=255)
                                                                                                                                              gggggCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATcc       <  1:1673584/150‑1 (MQ=255)
                                                                                                                                                gggCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGCGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCctc     <  2:53289/150‑1 (MQ=255)
                                                                                                                                                 ggCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGGGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCctct    >  2:23885/1‑150 (MQ=255)
                                                                                                                                                  gCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGGGTGGATGGTGGGTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTctg   <  1:1699114/150‑1 (MQ=255)
                                                                                                                                                   cAGGCGGCCCCCGTCTGTATGTCGCGGGTCTGGTGTTGATGGGGGTGTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGGGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCGCTCCCTTCCCATTATTGGCATTTATCCTctgc  <  2:1255904/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGAGCAGGGCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTCTGC  >  NC_000913/403965‑404261

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