Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 404,113 A→C 5.8% 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 5.8% 275.8 / 6.6 97A137A (GCA→GCCyaiCdiguanylate cyclase, cellulose regualtor
Reads supporting (aligned to +/- strand):  ref base A (46/45);  new base C (4/2);  total (50/47)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01

CGAGCAGGGCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTCTGC  >  NC_000913/403965‑404261
                                                                                                                                                    |                                                                                                                                                    
cGAGCAGGGCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAg                                                                                                                                                     >  2:1634455/1‑150 (MQ=255)
  aGCAGGGCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAggc                                                                                                                                                   >  2:435746/1‑150 (MQ=255)
       ggCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGccc                                                                                                                                              <  2:967788/150‑1 (MQ=255)
         ccGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGccccc                                                                                                                                            >  2:383971/1‑150 (MQ=255)
            tCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCCGCCTGCCCCCATc                                                                                                                                         <  1:1461372/150‑1 (MQ=255)
            gcGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCCGCCTGCCCCCATc                                                                                                                                         <  2:750930/149‑1 (MQ=255)
               aTCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGt                                                                                                                                      <  1:959251/150‑1 (MQ=255)
                  cGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTg                                                                                                                                   >  2:998690/1‑150 (MQ=255)
                  cGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTg                                                                                                                                   >  1:1548639/1‑150 (MQ=255)
                  cGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTg                                                                                                                                   <  2:512848/150‑1 (MQ=255)
                     ttAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTcg                                                                                                                                <  2:789439/150‑1 (MQ=255)
                        gCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCgg                                                                                                                             <  1:998690/150‑1 (MQ=255)
                        gCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCgg                                                                                                                             >  2:159384/1‑150 (MQ=255)
                         ccGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCggg                                                                                                                            <  2:439223/150‑1 (MQ=255)
                         ccGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCggg                                                                                                                            <  1:1577081/150‑1 (MQ=255)
                          cGGTAAATTTACAACGTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCCCCGCGATGTTGAGGATTATGTTTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCgggg                                                                                                                           >  2:1085262/1‑149 (MQ=255)
                            ggAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCt                                                                                                                         >  2:1352447/1‑150 (MQ=255)
                             gAAATTTACAACTTAAAAACCGATGCAGTATTCGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTg                                                                                                                        <  1:635152/150‑1 (MQ=255)
                             gAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTg                                                                                                                        <  2:1493153/150‑1 (MQ=255)
                              aaaTTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCATCTGTTTGTCGCGGGTCTgg                                                                                                                       <  2:507320/150‑1 (MQ=255)
                              aaaTTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGGAATGGGAGTAAACGTGCTGCCTTCCACCGCGAAGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGACTGTTTGTCGCGGGTCTgg                                                                                                                       >  2:758715/1‑150 (MQ=255)
                               aaTTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGgt                                                                                                                      <  1:383971/150‑1 (MQ=255)
                                  ttACAACTTAATACCCGATGCAGATTTAGGGGGACTGTGGGTAGGCGTAATGGGCGTTACCGTGCTGCCTTGCACCGCGAAGTTGATGATTATGTGTCTGAATTTGCTGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGtt                                                                                                                   <  2:1078108/150‑1 (MQ=255)
                                   tACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTg                                                                                                                  >  2:416320/1‑150 (MQ=255)
                                    aCAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGa                                                                                                                 >  2:1320963/1‑150 (MQ=255)
                                     cAACTTAATAACCGATGGAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCGGTTTGTCGCGGGTCGGGTGTTGAt                                                                                                                >  2:264834/1‑150 (MQ=255)
                                     cAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGAt                                                                                                                >  2:332302/1‑150 (MQ=255)
                                      aaCTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGAtg                                                                                                               >  2:1606960/1‑150 (MQ=255)
                                      aaCTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGAtg                                                                                                               <  1:992412/150‑1 (MQ=255)
                                       aCTTAAAAACCGATGCAGTATTAGCGGAAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGAtgg                                                                                                              >  2:833294/1‑150 (MQ=255)
                                          tAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGAtggtgg                                                                                                           >  2:533366/1‑150 (MQ=255)
                                            aaaaCCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGtt                                                                                                         >  2:1179117/1‑150 (MQ=255)
                                              aaCCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTGGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTc                                                                                                       <  1:103372/150‑1 (MQ=255)
                                              aaCCGATGCAGGATTAGCGGGAATGTGGGTCGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTCATGCGGGCAGGCGGCCCCCGGCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTc                                                                                                       <  2:1159804/150‑1 (MQ=255)
                                               aCCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTcc                                                                                                      <  1:990772/150‑1 (MQ=255)
                                                 cGATGCCGTATTAGCGGGACTGTGGGGAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTg                                                                                                    <  2:1292124/150‑1 (MQ=255)
                                                  gATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGGCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGc                                                                                                   <  1:1691080/150‑1 (MQ=255)
                                                       aGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGt                                                                                              >  1:1457678/1‑150 (MQ=255)
                                                       aGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGGCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGt                                                                                              >  2:175941/1‑150 (MQ=255)
                                                              gCGGGTATGTGGGGAGGCGTACTGGGCGTAAACCTGCGGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGTTGGGGGCAGTCGGCCCGCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTGTCCTGCCTTGTCACCCTc                                                                                       <  1:331503/150‑1 (MQ=255)
                                                               cGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCg                                                                                      <  1:1333574/150‑1 (MQ=255)
                                                                gggAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGa                                                                                     >  2:696271/1‑150 (MQ=255)
                                                                 ggAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAg                                                                                    >  1:750759/1‑150 (MQ=255)
                                                                   aaTGTGGGTAGGCGAAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGACTTTGCTGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCt                                                                                  <  2:669430/150‑1 (MQ=255)
                                                                     tgtgGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGa                                                                                <  1:701569/150‑1 (MQ=255)
                                                                     tgtgGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGa                                                                                >  2:1425753/1‑150 (MQ=255)
                                                                      gtgGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGAc                                                                               <  2:342067/150‑1 (MQ=255)
                                                                        gggTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACgg                                                                             <  1:1425753/150‑1 (MQ=255)
                                                                         ggTAGGCGTAATTGGCGTAAACGCGCTGCCTTCCCCCGCGATGTTGATGATTATGTGGCTGACTTTGCTGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTTGTTTCCTGCCTTGTCACCCTCGAGCTGACggg                                                                            <  1:1085262/150‑1 (MQ=255)
                                                                             ggCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCAtt                                                                        >  1:1455743/1‑150 (MQ=255)
                                                                              gCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTAGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTa                                                                       <  1:159384/150‑1 (MQ=255)
                                                                               cGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCGGTGTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCTTCGAGCTGACGGGCATTAc                                                                      >  2:843790/1‑150 (MQ=255)
                                                                                 tactcggCGTAAACGCGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATTGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCg                                                                    <  2:873187/145‑1 (MQ=255)
                                                                                   aTGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGCCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCgtg                                                                  >  1:1309781/1‑150 (MQ=255)
                                                                                     gggCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGCCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTc                                                                <  2:1309781/150‑1 (MQ=255)
                                                                                       gCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGt                                                              >  2:771871/1‑150 (MQ=255)
                                                                                         gTACACGTGCGGCCTTCCACCGCGATGTTGATGATTCTGTGTCTGAATTTGATGGTGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGGTGATGGTGGTTTCCTGCCTTGTCACCCTCGCGCTGACGGGCATTACCGTGTCGTTc                                                            <  2:1058306/150‑1 (MQ=255)
                                                                                         gTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGTATTTGATGGGGGCAGGCGGCCCCCGTCTGTTGGTCGCGGGTCTTGTGGTGATGGGGGTGTCCTGGCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTc                                                            <  1:684405/150‑1 (MQ=255)
                                                                                          tAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCGGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGCCGGGCATTACCGTGTCGTTCa                                                           >  2:913853/1‑150 (MQ=255)
                                                                                              cGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAg                                                       >  1:711362/1‑150 (MQ=255)
                                                                                                tgctgcCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTg                                                     >  2:550230/1‑150 (MQ=255)
                                                                                                  ctgcCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTgcg                                                   >  1:701068/1‑150 (MQ=255)
                                                                                                     ccTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGcgccg                                                >  1:1718951/1‑150 (MQ=255)
                                                                                                        tCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTg                                             >  2:781173/1‑150 (MQ=255)
                                                                                                            ccGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAAt                                         >  2:1689234/1‑150 (MQ=255)
                                                                                                            ccGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAAt                                         >  1:1567490/1‑150 (MQ=255)
                                                                                                              gcgATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAAtgg                                       >  1:1095553/1‑150 (MQ=255)
                                                                                                                gATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACGCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAAtggtg                                     <  1:1440363/150‑1 (MQ=255)
                                                                                                                gATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAAtggtg                                     <  1:416320/150‑1 (MQ=255)
                                                                                                                 aTGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAAtggtgg                                    <  1:696271/150‑1 (MQ=255)
                                                                                                                 aTGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAAtggtgg                                    >  2:1163948/1‑150 (MQ=255)
                                                                                                                     tgatgaTTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGctct                                >  1:981979/1‑150 (MQ=255)
                                                                                                                      gatgatTATGTGACTGAATTTGCTGGGGGCATGCGGCCCCCGTCTGTGTGTCGCGGGTCTGGTGTTGATGGTGGTTGCCTGCCTTGTCACCCTCGAGATGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCtctc                               <  1:1036777/150‑1 (MQ=255)
                                                                                                                       atgatTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTcc                              <  1:781173/150‑1 (MQ=255)
                                                                                                                          atTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCtt                           <  1:1163948/150‑1 (MQ=255)
                                                                                                                           ttctgtgTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGAGGGTGGTTTCCTGCCTTGGCACCCTCGCGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTc                          <  1:913853/147‑1 (MQ=255)
                                                                                                                             aTGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTccc                        <  1:1352447/150‑1 (MQ=255)
                                                                                                                             aTGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTccc                        >  2:876280/1‑150 (MQ=255)
                                                                                                                                 gtCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATTGTGGTTTCCTGCCTTGTCACCCTCTAGCTGACGGGCATTACCGTGTCGATCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCaata                    >  2:644006/1‑150 (MQ=255)
                                                                                                                                   cTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCattatt                  >  2:979547/1‑150 (MQ=255)
                                                                                                                                     gAATTTGATGGGGGCAGGCGGCCCGCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCAGGCTTGTCACCGTCCAGCAGTCGGGGATTTCGGTGTCGTTGAATAGTGGGCCGCTGGTATGGTGGGTGTCGCGTCCGATTTTTGt                >  1:1251399/1‑150 (MQ=255)
                                                                                                                                       aTTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCTATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCa              <  1:332302/150‑1 (MQ=255)
                                                                                                                                       aTTTGATGGGGGCAGGCCGCCCCCGACTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCAACCGCGAGCTGACGGGCATTACCGTGTCGATCAATAGTGCGCCGCTTGCATGCTGGCTCTCCCTTCCCATAATTGTCa              >  1:259277/1‑150 (MQ=255)
                                                                                                                                        tttGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCAt             >  1:548812/1‑150 (MQ=255)
                                                                                                                                        acagacgggggCCGCCTGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTGCCATTATTGTCAt             >  1:1015347/7‑150 (MQ=255)
                                                                                                                                         cagacgggggCCGCCTGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCAtt            >  1:1663377/6‑150 (MQ=255)
                                                                                                                                           gatcgcggCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGCTGTTGATGGTGGATTCCAGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTa          <  2:729135/144‑1 (MQ=255)
                                                                                                                                           gacgggggCCGCCTGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGAGGGTGGGTTCCTGCCTTGTCCCCCTCGAGCTGACGGGCATTACCGTGGCGTTCACTCGTGCGCCGCTGGAATGGTGGGTCTCCCTTCGCAGTATTGTCAgttt          >  1:1545022/4‑149 (MQ=255)
                                                                                                                                            aTGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCCCCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCACTGGATTGGTGGCTATCCCTTCCCATTATTGTATTTTaa         >  2:320992/1‑149 (MQ=255)
                                                                                                                                             tGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATc        >  1:969868/1‑150 (MQ=255)
                                                                                                                                             tGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATc        >  1:201267/1‑150 (MQ=255)
                                                                                                                                             tGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATc        <  1:1606960/150‑1 (MQ=255)
                                                                                                                                             cgggggCCGCCTGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATc        >  2:708850/2‑150 (MQ=255)
                                                                                                                                              gggggCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATcc       <  1:533366/150‑1 (MQ=255)
                                                                                                                                              gggggCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATcc       <  2:1457678/150‑1 (MQ=255)
                                                                                                                                                gggCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGGTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCctc     <  2:1246679/150‑1 (MQ=255)
                                                                                                                                                   cAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTctgc  <  1:771871/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: