Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 404,117 G→T 6.3% G139C (GGC→TGC)  yaiC → diguanylate cyclase, cellulose regualtor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913404,1170GT 6.3% 296.1 / 7.1 95G139C (GGC→TGC) yaiCdiguanylate cyclase, cellulose regualtor
Reads supporting (aligned to +/- strand):  ref base G (43/45);  new base T (4/2);  total (48/47)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

GGCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTCTGCTG  >  NC_000913/403972‑404263
                                                                                                                                                 |                                                                                                                                                  
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                                                                                                                                     <  2:512848/150‑1 (MQ=255)
           cGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTg                                                                                                                                     >  1:1548639/1‑150 (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                                                                                  >  2:1425753/1‑150 (MQ=255)
                                                              tgtgGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGa                                                                                  <  1:701569/150‑1 (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                                           >  1:1567490/1‑150 (MQ=255)
                                                                                                     ccGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAAt                                           >  2:1689234/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                                      >  2:1163948/1‑150 (MQ=255)
                                                                                                          aTGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAAtggtgg                                      <  1:696271/150‑1 (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:1606960/150‑1 (MQ=255)
                                                                                                                                      tGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATc          >  1:201267/1‑150 (MQ=255)
                                                                                                                                      tGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATc          >  1:969868/1‑150 (MQ=255)
                                                                                                                                      cgggggCCGCCTGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATc          >  2:708850/2‑150 (MQ=255)
                                                                                                                                       gggggCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATcc         <  2:1457678/150‑1 (MQ=255)
                                                                                                                                       gggggCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATcc         <  1:533366/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)
                                                                                                                                                 |                                                                                                                                                  
GGCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTCTGCTG  >  NC_000913/403972‑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: