Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 404,115 G→C 6.3% 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 6.3% 313.2 / 9.9 95G138A (GGC→GCC) yaiCdiguanylate cyclase, cellulose regualtor
Reads supporting (aligned to +/- strand):  ref base G (44/44);  new base C (4/2);  total (48/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 = 4.64e-01

AGCAGGGCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTCTGCTG  >  NC_000913/403967‑404263
                                                                                                                                                    |                                                                                                                                                    
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                          >  2:876280/1‑150 (MQ=255)
                                                                                                                           aTGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTccc                          <  1:1352447/150‑1 (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:1606960/150‑1 (MQ=255)
                                                                                                                                           tGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATc          >  1:201267/1‑150 (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)
                                                                                                                                                   ggcggcCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTctgctg  >  2:1444913/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGCAGGGCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTCTGCTG  >  NC_000913/403967‑404263

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: