Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,876,476	T→C	E159E (GAA→GAG)	yidA ←	sugar phosphatase YidA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,876,476	0	T	C	100.0%	163.3 / NA	51	E159E (GAA→GAG)	yidA	sugar phosphatase YidA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (34/17);  total (34/17)

CCATCGCCACACCGACGCCTGCATATTCAATCATTGCGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCG  >  NC_000913/3876241‑3876722                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                       ccATCGCCACACCGACGCCTGCATATTCAATCATTGCGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAatcatcatc                                                                                                                                                                                                                                         >  1:103387/1‑251 (MQ=255)     agCCACACCGACGCCTGCATATTCAATCATTGCGATATCGTTTTCCTGATCGACAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTACCGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACtt                                                                                                                                                                                                                                     >  2:222347/2‑251 (MQ=255)       ccACACCGACGCCTGCATATTCAATCATTGCGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTc                                                                                                                                                                                                                                   <  2:368889/251‑1 (MQ=255)         acacCGACGCCTGCATATTCAATCATTGCGATATCGTTTTCCTGCTCGACAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTACCGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGCTAGCCTGGTCGCGGCTGGCGTGCTCCTCCATCATCATCACTTTc                                                                                                                                                                                                                                   >  2:234141/1‑249 (MQ=255)                     gCATATTCAATCATTGCGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATAACTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTAtt                                                                                                                                                                                                                     >  1:320320/1‑251 (MQ=255)                       atatTCAATCATTGCGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTgg                                                                                                                                                                                                                   >  2:365521/1‑251 (MQ=255)                              aTCATTGCGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTAGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTAGAGCAAGTACGGCGAACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCAt                                                                                                                                                                                                            >  1:192304/1‑251 (MQ=255)                                            gTTTTCCTGATCGCCAATCGCATTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGGATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCATTCATCATCATTTTCAGGCACTGGGTATTGGGGGCCATTTTCTCGGCTTCGc                                                                                                                                                                                              >  2:193920/1‑251 (MQ=255)                                            gTTTTCATGATCGCCAATCGCCATGACTGCTTCCGGTTTAATACCTTACACGGCGACCCGTGATTTCACCCCCGTACCTTTTTTAACGCGCTTATCGCGGGTTACGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGAGGTATACGCGCGATAGCCTGGTCGAGGATGGCCGGCTCCTCACTCCTCATCACTTTCAGGCACTGGGTATTGGGGTCCCTTTTTTCCggctcgc                                                                                                                                                                                              >  2:70100/1‑245 (MQ=255)                                                      tcgccattcgccatGATTTCCTCCGGATTAATCCCTATGACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTTACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCTCGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCATTCATCATCACTTTCAGGAACTGGGTATTGGGGCCCATTTTCTCCGCTTCGCAGAACAcc                                                                                                                                                                                      >  2:445968/1‑249 (MQ=255)                                                        gccaatcgccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCg                                                                                                                                                                                  >  2:396118/1‑251 (MQ=255)                                                            atcgccaTGATTTCGTCCGGTTTAATTCCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACGGGGTATTTGTGTCCATTTTTTCCGCTGCGCAGAACACCAGCGGAAt                                                                                                                                                                              <  1:410372/251‑1 (MQ=255)                                                                              ggTTTAAGACCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTc                                                                                                                                                            <  1:234141/251‑1 (MQ=255)                                                                                 ttAATACCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGCATTGGGGTCCATTTTCTCCGCTTCCCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATg                                                                                                                                                         >  1:178490/1‑251 (MQ=255)                                                                                   aaTACCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCa                                                                                                                                                       <  2:159659/251‑1 (MQ=255)                                                                                   aaTACCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGGCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGGGGCAACGAAGGATTCATGc                                                                                                                                                        >  2:335962/1‑250 (MQ=255)                                                                                   aaTACCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGGATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAGGGATTCATGcc                                                                                                                                                       >  2:351277/1‑250 (MQ=255)                                                                                            cacGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTgtagt                                                                                                                                              >  1:300264/1‑251 (MQ=255)                                                                                               gTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGACGGATTCATGCACCGTGTAGTAGc                                                                                                                                           >  2:185200/1‑251 (MQ=255)                                                                                                         gATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACg                                                                                                                                 >  1:49299/1‑251 (MQ=255)                                                                                                               aCCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGc                                                                                                                           >  1:51172/1‑251 (MQ=255)                                                                                                                 ccccGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAAtgg                                                                                                                                                                            <  2:92669/200‑1 (MQ=255)                                                                                                                 ccccGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAAtgg                                                                                                                                                                            >  1:92669/1‑200 (MQ=255)                                                                                                                                       ttatcaaggtcttCGAGGAAGTACGGAGAATTTTTCAGCACGGTCAATTTCTCTCTCACTTCATGCGGAATACGCGCGATCGCCTGGTCGAGGAAGGCGGGCTCATAAATCATCATCACTTTCAGGAACTGGGTATTGGGGCCCATTTTCTCCGATTTGCAGAAAACCAGCGGAATGGGTGCAACGAAGGATTCATGCACCGTGTCGTAGCTGCTATCCCGGTTGGCGGTGTACAGCGTGGTGAGGTcc                                                                                                     >  2:228652/12‑249 (MQ=255)                                                                                                                                         aTCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAggg                                                                                                 >  1:14240/1‑251 (MQ=255)                                                                                                                                            gAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGt                                                                                              >  1:348531/1‑251 (MQ=255)                                                                                                                                                     gAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATg                                                                                       >  2:112647/1‑249 (MQ=255)                                                                                                                                                          aGTACGGCGCACTTTTCAGCACGGTATATTGCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTTCAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCg                                                                                <  1:39833/251‑1 (MQ=255)                                                                                                                                                              cGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCAGCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTTATATCACGGTTGGCGGTGGACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACtt                                                                            <  1:436232/251‑1 (MQ=255)                                                                                                                                                                 cgcACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTcgc                                                                         <  1:177983/251‑1 (MQ=255)                                                                                                                                                                    aCTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGCGGATGGCGGGCTCATCAATCATCATCCCTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTcgcg                                                                        >  1:450498/1‑249 (MQ=255)                                                                                                                                                                     cTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTcgcg                                                                        >  1:147861/1‑248 (MQ=255)                                                                                                                                                                           aGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGttt                                                               >  1:62751/1‑251 (MQ=255)                                                                                                                                                                            gCACGGTATATTTCTCGGTCAGTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCGCTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATTGTGGCAACGAAGGATTCATGCGCCGTGTAGTAGCTGATAGCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGtttt                                                              <  1:421118/251‑1 (MQ=255)                                                                                                                                                                              aCGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCTGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTc                                                            <  1:193920/251‑1 (MQ=255)                                                                                                                                                                                ggTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTcc                                                           >  1:61711/1‑250 (MQ=255)                                                                                                                                                                                        tctcTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTcgatagcttctgaatcgcgaccgtt                                                   >  1:464723/1‑227 (MQ=255)                                                                                                                                                                                             ttCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCCATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGAAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAgtc                                             >  1:333446/1‑251 (MQ=255)                                                                                                                                                                                             ttCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATa                                                <  2:348531/248‑1 (MQ=255)                                                                                                                                                                                              tcccttgcggCGGAATACGCGCGAGAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGGATTGGGGTCCATTGTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGAGATCACGGTTGGCGGTGTACAGCGTGGGGCGGTCCAGGGCTTGGAAATGAGAGCCGACTTCGCGAGAGAGTTTTTCCAGGAAACGTTAgtc                                             <  1:74074/241‑1 (MQ=255)                                                                                                                                                                                               cACTTCCTGCGGAATACGAGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAAGTGGGGATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGCTTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGAGTGGTGCGGGCCAGGGCGTGGAAATGAGAACCGAATTCGCGAGAGAGTTTTTCCCGGAAACGATAggcgg                                           >  2:356938/1‑247 (MQ=255)                                                                                                                                                                                                ccTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAgtcgtc                                          <  2:14240/250‑1 (MQ=255)                                                                                                                                                                                                        cGGAATACGCGCGATCGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTg                                   >  1:361195/1‑250 (MQ=255)                                                                                                                                                                                                         ggAATACGCGCGATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTg                                   >  1:234379/1‑249 (MQ=255)                                                                                                                                                                                                                gcgcgATAGCCTGGTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGtt                          <  2:51172/251‑1 (MQ=255)                                                                                                                                                                                                                            ggTCGAGGATGGCGGGCTCATCAATCATCATCACTTTCACGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTg              >  1:470266/1‑251 (MQ=255)                                                                                                                                                                                                                                 gtgATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGGACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTGCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTAcc         <  1:286797/249‑1 (MQ=255)                                                                                                                                                                                                                                 aGGATGGCGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTAcc         >  1:125278/1‑251 (MQ=255)                                                                                                                                                                                                                                     tGGCGGGCTCATCACTCATCATCACTTTCAGGAACTGGGGATTGGGGTCCCTTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGCAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCCCGGTGCTACCATcc     >  2:335757/1‑250 (MQ=255)                                                                                                                                                                                                                                        cGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAGCACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCg  <  2:108992/251‑1 (MQ=255)                                                                                                                                                                                                                                        cGGGCTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCg  <  1:406958/251‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                       CCATCGCCACACCGACGCCTGCATATTCAATCATTGCGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCG  >  NC_000913/3876241‑3876722

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