Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,116,703	T→G	K173N (AAA→AAC)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,116,703	0	T	G	100.0%	217.1 / NA	68	K173N (AAA→AAC)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (44/24);  total (44/24)

GAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGG  >  NC_000913/4116560‑4116847                                                                                                                                                 |                                                                                                                                                 gaaCGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCa                                                                                                                                             <  2:1303101/149‑1 (MQ=255) gaaCGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAAAGGCGTCGGTGTAAA‑CGGTCACATGGACATGGCCCTGTGTCATGTTCCa                                                                                                                                             <  2:128068/149‑1 (MQ=255)  aaCGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAg                                                                                                                                            <  1:1333906/149‑1 (MQ=255)      accgccTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACCACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACTCGGCCCTGAGTCATGTTCCAGAGaa                                                                                                                                        <  2:365866/144‑1 (MQ=255)          aCTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCa                                                                                                                                    >  1:1368531/1‑149 (MQ=255)                  cagcaTCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCaa                                                                                                                            >  2:1193859/1‑149 (MQ=255)                    gcaTCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAAcc                                                                                                                          <  1:692050/149‑1 (MQ=255)                       tctcGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTa                                                                                                                       <  2:982382/149‑1 (MQ=255)                       tctcGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGGATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTa                                                                                                                       >  1:885731/1‑149 (MQ=255)                       tctcGCGCGGAATATCCAGCACTTCCAACATTTTGTAGTCCCAGTCCAGGGTATGGATGTTGCACAACATGGTACGAGAGGCGTTGGTGTAAA‑CGGGCACATGGACACGGCCCTGAGTCATGTTCCAGATAAACCACGTCTCAcccctc                                                                                                                       >  2:16598/1‑146 (MQ=255)                             gcgGAATATCCAGCACTTCCAGCAATTTGTAGTCCCAGTCCAAGGTATAGATATTGAACAACAGCGTACACAACACGTTGCTGTACT‑CGGACACATGGAAACGACCCCGAGTGTTGTTCTAGATAATGCACCTGAAACCCGTACCAAAc                                                                                                                 >  2:596265/1‑149 (MQ=255)                              cgGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCGGATAAGCCACGTATCAACCGTACCAAACa                                                                                                                >  2:1222627/1‑149 (MQ=255)                              cgGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACa                                                                                                                >  1:359608/1‑149 (MQ=255)                               ggAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGCCACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCCACCGTACCACACAg                                                                                                               <  2:877047/149‑1 (MQ=255)                                gAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGc                                                                                                              >  1:844047/1‑149 (MQ=255)                                 aaTATCCCGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCa                                                                                                             <  1:787394/149‑1 (MQ=255)                                  atatCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCaa                                                                                                            <  2:136371/149‑1 (MQ=255)                                    atCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAAtt                                                                                                          <  1:169062/149‑1 (MQ=255)                                      ccAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGCAGGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTca                                                                                                        >  1:618820/1‑149 (MQ=255)                                          cACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCAGGTATTAACCGTACCAAAAAGCAATTcaccac                                                                                                    >  2:1236642/1‑149 (MQ=255)                                             ttCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGCACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCccgac                                                                                                 >  2:323449/1‑149 (MQ=255)                                             ttCCAGCATTTTGACGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCacgac                                                                                                 >  2:1205292/1‑149 (MQ=255)                                                cAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTg                                                                                              >  1:1143919/1‑149 (MQ=255)                                                cAGCATTTTGTCGTCCCAGTCCAGGGAATGCATTACAAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCCTGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTg                                                                                              >  1:1380651/1‑149 (MQ=255)                                                 aGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGc                                                                                             >  2:535820/1‑149 (MQ=255)                                                   cATTTTGTCGTCCCAGTACAGGGTATGGATGTTGAACAACATGGTACGAGGGGCGTTGGTGTAAT‑CGGTCACATGGAGACGGCCCTGAGTCATGTTCCAGATGAGCCGCGTATCAACCGTACCAACCAGCACTTCACACCGACGTGCAc                                                                                           >  1:713782/1‑149 (MQ=255)                                                            gtcCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCCCGACGTGCACGCTCCCgag                                                                                  >  1:88730/1‑149 (MQ=255)                                                            gtcCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCCCGCaag                                                                                  >  1:89300/1‑149 (MQ=255)                                                             tcCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCTCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACAAGGACGTGCACGCTCGCgaga                                                                                 >  1:231304/1‑149 (MQ=255)                                                             tcCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCacga                                                                                                  <  1:758252/132‑1 (MQ=255)                                                             tcCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCacga                                                                                                  >  2:758252/1‑132 (MQ=255)                                                             tcCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCgaga                                                                                 >  1:1248889/1‑149 (MQ=255)                                                             tcCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCgaga                                                                                 >  1:54658/1‑149 (MQ=255)                                                             tcCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCgaga                                                                                 >  1:1337697/1‑149 (MQ=255)                                                             tcCCAGCCCAGGTTATGGGTGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACGTTGAGAGGGCCCTGAGTCCTGTGCGAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCgaga                                                                                 <  1:1021315/149‑1 (MQ=255)                                                              ccgaggtcAGGGTATGGGTGTTGAGCAACATGGTACGAGAGGCGTGGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGGCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGagag                                                                                <  2:1301810/142‑1 (MQ=255)                                                                      aGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCCGAGTCACGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAACCCACCACGACGTGCACACCCCCCACAGCCCTCcac                                                                        >  2:353188/1‑149 (MQ=255)                                                                        ggTATGGATGTTGAACAACATGGTACGAGAGGTGTT‑GTGGACTGCGGTCATATGGGCACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCACTTCCCCAGGATGTGCACGATCGCGCGAGCCTTCCCCAt                                                                      <  2:1331314/149‑1 (MQ=255)                                                                           aTGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGt                                                                   <  1:1251854/149‑1 (MQ=255)                                                                             ggATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACAGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAAAGCCTTCCACATggacc                                                                 >  2:498915/1‑146 (MQ=255)                                                                              gATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGa                                                                >  1:715778/1‑149 (MQ=255)                                                                                     aacgacATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCAACGTATCAACCGTACCAAACCGCAATTCACCACGCCGTGCAGGCTCGCGAGAGCCTTCCACATAGTCCAGGATCca                                                         >  1:634498/1‑149 (MQ=255)                                                                                           aTGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCCACCGTACCAAACAGCAATTCACCACGACGTGCAGGCTAGCGAGAGCCTTCCACATGGTCGAGGATCcacttcaa                                                   >  1:351336/1‑148 (MQ=255)                                                                                           aTGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCcacttcac                                                   <  1:649816/149‑1 (MQ=255)                                                                                               tACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTg                                               >  1:884591/1‑149 (MQ=255)                                                                                                  gagaGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGtt                                            >  1:58576/1‑148 (MQ=255)                                                                                                  gagaGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTg                                            >  1:408528/1‑149 (MQ=255)                                                                                                     agGCGTTGGTGTAAT‑CGGTCACATGGGCACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCa                                         <  2:1339336/149‑1 (MQ=255)                                                                                                          ttGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGaaaa                                    <  1:1153359/149‑1 (MQ=255)                                                                                                          ttGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGaaaa                                    <  1:1193859/149‑1 (MQ=255)                                                                                                             gtgtAAT‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCCCTTCACTTTGGTGCCAGAAAAGTa                                 >  2:656074/1‑149 (MQ=255)                                                                                                                   t‑CGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCGTCCACATAGTCGAGGATCCCCCTCACTTCGGTGCCAGGAAACTACGTGTc                           >  1:599318/1‑149 (MQ=255)                                                                                                                       gTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAAt                        >  1:1219377/1‑149 (MQ=255)                                                                                                                         cacaTGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCa                      <  1:1205292/149‑1 (MQ=255)                                                                                                                            aTGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAGCAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGATGCCAGAAAAGTACGGGTCAATCACCa                   >  1:794505/1‑149 (MQ=255)                                                                                                                            aTGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAGCAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTGCCACATGGTCGGGGATCCACTTCACTTTGATGCCAGAAAAGTACGGGTCAATCGCCa                   >  1:794085/1‑149 (MQ=255)                                                                                                                            aTGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCa                   >  2:924944/1‑149 (MQ=255)                                                                                                                               gACACGGCCCGGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGAc                <  2:1175950/149‑1 (MQ=255)                                                                                                                                acacGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCCCCAGAcc               >  2:1018305/1‑149 (MQ=255)                                                                                                                                acacGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGAcc               >  2:66076/1‑149 (MQ=255)                                                                                                                                 cacGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCg              >  1:393957/1‑149 (MQ=255)                                                                                                                                  acGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCgg             >  2:112404/1‑149 (MQ=255)                                                                                                                                     gCCCTGAGGCAGGGTCCAGATAGGGCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTGAATCTCTAGACCGTTAt          <  2:167874/149‑1 (MQ=255)                                                                                                                                       ccTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATtg        >  2:1033453/1‑149 (MQ=255)                                                                                                                                          gAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGTGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAATAGTACGGGTCAATAACCAGACCGGTATtgctg     >  2:782557/1‑149 (MQ=255)                                                                                                                                           gctCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATtgctgc    <  1:315308/147‑1 (MQ=255)                                                                                                                                            gTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCg   <  1:353188/149‑1 (MQ=255)                                                                                                                                             tCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCgg  >  1:584568/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                 GAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAAT‑CGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGG  >  NC_000913/4116560‑4116847

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑