Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,311,337	C→T	A236T (GCA→ACA)	alsB ←	D‑allose ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,311,337	0	C	T	100.0%	134.8 / NA	38	A236T (GCA→ACA)	alsB	D‑allose ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/14);  total (24/14)

CATTCAGTTTTATTGAGTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGCGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAGGCGTTACCGGCTTTACCCTCAATGATTGCGACTTCAC  >  NC_000913/4311098‑4311573                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             catTCAGTTTTATTGAGTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGAGTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGa                                                                                                                                                                                                                                    <  1:195243/250‑1 (MQ=255)  atTCAGTTTTATTGAGTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTGTTTCCGGCGTTTGTGACTGCCTGAGc                                                                                                                                                                                                                                  <  1:54279/251‑1 (MQ=255)           tATTGAGTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCAt                                                                                                                                                                                                                         >  2:277756/1‑251 (MQ=255)             ttGAGTGACCAGGCGTGGATCGACCAGTTTAAATTCCGGCGCTTTATCCAGCGGGATAACCTTGCCGGATTTATCAACGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCAATTTTCCCGTTTTTCCGACGTTTGTGACTGCCTGAGCAACACCCATTg                                                                                                                                                                                                                       >  2:33667/1‑251 (MQ=255)              tGAGTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTg                                                                                                                                                                                                                       <  1:70933/250‑1 (MQ=255)              tGAGTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACAGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGc                                                                                                                                                                                                                      >  1:130573/1‑251 (MQ=255)                             gAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTCCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTcgc                                                                                                                                                                                                       >  1:27748/1‑251 (MQ=255)                                  gACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAAt                                                                                                                                                                                                  <  1:336166/251‑1 (MQ=255)                                               tCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTaa                                                                                                                                                                                      <  1:327841/250‑1 (MQ=255)                                               tCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAat                                                                                                                                                                                     >  1:60970/1‑251 (MQ=255)                                                          aTCCAGCGGGATAACCTTTCCGGATTTCACAGCGTCAACCATCAGCTTCAGACCCGTTGCGGTGATGTCCTCCTGGTTCTGGGCAACCGTCGCGGTCATTTTTCCGGCTTCCACCATTTTGCGGTCTTCCGGAATGCCATCTGTCCCGACGACTAGCACTTTTCCCGTGTTTCCGGCGTTGGTGACTGCCTGTGCAACACCCATTTCCATGGCGTCATTCGCGCAATAGATCCCTTTATTTTTCTGATTAc                                                                                                                                                                          <  2:374033/251‑1 (MQ=255)                                                                 gggATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCCGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCa                                                                                                                                                                    >  2:164755/1‑250 (MQ=255)                                                                        ccTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTa                                                                                                                                                            <  1:164755/251‑1 (MQ=255)                                                                                gATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGAc                                                                                                                                                    >  1:279552/1‑251 (MQ=255)                                                                                         gCGTCAACAATCAGCTTCAGACCAGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTACCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCCTCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGc                                                                                                                                           >  2:45689/1‑251 (MQ=255)                                                                                             cAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTa                                                                                                                                       <  2:225009/251‑1 (MQ=255)                                                                                                  aTCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCg                                                                                                                                  >  1:278627/1‑251 (MQ=255)                                                                                                      gCTTCAGACCCGTTGCGCCGATATCCGCCGGTTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTGCCGACGACCAGCACTTTTCCCGTTTTTCCGTCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAata                                                                                                                                                                                    <  1:180648/197‑1 (MQ=255)                                                                                                      gCTTCACACACGTTGCGCCGATATCCCCCGGGCTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGCCTGCCTGAGCAACACCCATTGCCCTCGTGTCATTCGCGCAATAGATCGCTTTAata                                                                                                                                                                                    >  2:180648/1‑197 (MQ=255)                                                                                                        ttCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCCTTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATTCGGTccc                                                                                                                             >  1:237779/1‑250 (MQ=255)                                                                                                                     cgcCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCCCTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAggctgg                                                                                                               >  1:261781/1‑251 (MQ=255)                                                                                                                            atCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCCTTCTGCGGGCTTCCGGAATGCAATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCg                                                                                                         >  1:202981/1‑250 (MQ=255)                                                                                                                            atCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCg                                                                                                         >  2:351519/1‑250 (MQ=255)                                                                                                                                    ggTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTg                                                                                                >  2:178148/1‑251 (MQ=255)                                                                                                                                                    cgcgGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCttttt                                                                                <  1:178148/251‑1 (MQ=255)                                                                                                                                                    cgcgGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCttttt                                                                                >  1:150870/1‑251 (MQ=255)                                                                                                                                                        gacatatGTCCGGATTCCACCATTTTGCGGGCTTCCGGCATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGCTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTCATGCGGTCCCAGTCGGCAGTCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGa                                                                             >  2:66214/7‑250 (MQ=255)                                                                                                                                                                  cGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGCCAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTg                                                                  >  1:322968/1‑251 (MQ=255)                                                                                                                                                                   ggCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCTACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTgg                                                                 <  1:241072/251‑1 (MQ=255)                                                                                                                                                                       tCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTGTTTGAAGGCTTCGGTGGCAc                                                              <  2:251165/250‑1 (MQ=255)                                                                                                                                                                         caccaTTTTGCGGGCTTCCGGAATGCCATCTGTACCGTCGACCAGCTCTTTGCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCTGTCCCAGTCGGCAGGTTTGCTGGCGTCAAGCTTGATCTGGCTTGCTGTTTTGAAGGCTTCGGTGGGACCAt                                                           <  2:198971/251‑1 (MQ=255)                                                                                                                                                                                           cGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCgg                                         >  1:56760/1‑251 (MQ=255)                                                                                                                                                                                            ggAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCgg                                         >  2:123662/1‑250 (MQ=255)                                                                                                                                                                                             gAATGCCATCTGTACCGACGGCCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAg                                       <  2:278627/251‑1 (MQ=255)                                                                                                                                                                                                               cgacCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGCTTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAGGCGTTACCGGCTTTAccc                     >  2:187380/1‑251 (MQ=255)                                                                                                                                                                                                                              ccGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGCGGCGTTACCGGCTTTACCCTCAATGATTGCGACt      >  1:206678/1‑251 (MQ=255)                                                                                                                                                                                                                                gTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAGGCGTTACCGGCTTTACCCTCAATGATTGCGACTTc    >  1:224675/1‑251 (MQ=255)                                                                                                                                                                                                                                  ttttCCGGCGTTTGTGACTGCCTGAGCAAAACACATTGCCATCGTTTCATGCGCGCAATAGCTCGCTCTCATATTCGGCTTACGTTGCAACACGTTAGTGGCGACATCTAGTGCTTTAATGCGGTCCCAGTAGGCAGGCTGGCAGGCGACAAGCTTGATCAGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGCGGCGTTACCGGTTTTACCCTCCATGCTTGAGACTTcac  >  2:143145/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             CATTCAGTTTTATTGAGTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGCGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAGGCGTTACCGGCTTTACCCTCAATGATTGCGACTTCAC  >  NC_000913/4311098‑4311573

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