Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 404,115 G→C 5.4% G138A (GGC→GCC)  yaiC → diguanylate cyclase, cellulose regualtor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913404,1150GC 5.4% 332.7 / 6.2 106G138A (GGC→GCC) yaiCdiguanylate cyclase, cellulose regualtor
Reads supporting (aligned to +/- strand):  ref base G (54/45);  new base C (2/4);  total (56/50)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

GGGCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTCTGC  >  NC_000913/403971‑404261
                                                                                                                                                |                                                                                                                                                  
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)
   ccGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGCGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGccccc                                                                                                                                            >  2:627855/1‑150 (MQ=255)
    cGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCg                                                                                                                                           >  1:1849768/1‑150 (MQ=255)
     gTCGATCCGCTTAGCCTGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATTTGTCTGAATTTGATGGGGGCCGCCTGcccccat                                                                                                                                          <  1:232603/150‑3 (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)
               ttAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTcg                                                                                                                                <  1:676502/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)
                             tACATCTTAAAAACCGATGCAGTATTAGCGGTAATGTGGGTAGGCGTAATGTGCGTAATCGTGCTGCCTTCCACGGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGGCGCTGGTCTGGTGTTg                                                                                                                  <  2:1330141/150‑1 (MQ=255)
                               cAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGAt                                                                                                                >  1:1306656/1‑150 (MQ=255)
                               cAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGAt                                                                                                                >  1:1013036/1‑150 (MQ=255)
                               cAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGAt                                                                                                                >  1:961108/1‑150 (MQ=255)
                               cAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGAt                                                                                                                <  2:1493714/150‑1 (MQ=255)
                                 aCTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGAtgg                                                                                                              <  2:1446641/150‑1 (MQ=255)
                                 aCTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGAtgg                                                                                                              >  1:1442429/1‑150 (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                                                                                               <  1:927533/150‑1 (MQ=255)
                                                cAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTg                                                                                               >  2:1593179/1‑150 (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)
                                                                                           gctgcCTTCCACCGCGATTGTGATGATTATGTGTCTGAATTTGATGGGGGCCGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGGTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTgc                                                    <  2:146149/150‑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)
                                                                                                cTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGcgccgc                                               >  1:1562996/1‑150 (MQ=255)
                                                                                                cTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGcgccgc                                               >  2:1014657/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)
                                                                                                                            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)
                                                                                                                                     gacgggggCCGCCTGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTa          >  2:235627/4‑150 (MQ=255)
                                                                                                                                        gggggCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATcc       <  1:1673584/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)
                                                                                                                                                |                                                                                                                                                  
GGGCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTCTGC  >  NC_000913/403971‑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: