Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,992,970	T→G	R606R (CGT→CGG)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,992,970	0	T	G	98.4%	182.1 / ‑3.7	63	R606R (CGT→CGG)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (33/29);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.71e-01

TCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTG  >  NC_000913/3992723‑3993209                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                gcGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACg                                                                                                                                                                                                                                              <  1:475365/250‑1 (MQ=255)    ccGAAGTGTCGCACCATTTCCCGCTGCGCTTACCTGCGCCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCAGCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAAGTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTg                                                                                                                                                                                                                                           <  2:212446/251‑1 (MQ=255)    acgccGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCTTCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGAATTCCGTAAGATGGATGTCTTCAGCTTTGGCGAGAGTCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGAGGAACGACGTGCGGACGCTg                                                                                                                                                                                                                                           >  2:342336/5‑251 (MQ=255)        cGTGTCGAACCATTTCCCGCTGCGCTTACCTGCACCGACCCCGAAGGCGCTCTACAGCCCGAGCGAGATCCGCCATCTGGCGCTTATCGTTCACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGGGAGAATCAGAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACt                                                                                                                                                                                                                                       >  1:311054/1‑251 (MQ=255)         gtgtCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGACTATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACtt                                                                                                                                                                                                                                      >  2:517662/1‑251 (MQ=255)           gtCGCACCATTTCCCGCGGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTCACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAATATTGCCTAGTAGGTAGCGTTGACCTGCTGGACCGCACCTCGTGGCACGAGGAGCGGACGCTTCACTTc                                                                                                                                                                                                                                     >  2:252224/1‑250 (MQ=255)             cgcccctttgccctctccgcTTACCGGCCCCGACACCGAAGTCGTTCTACAGCCCGTGTGAGATCCGCCAGCGGGCGATGATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCTGAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCCCTTCAAc                                                                                                                                                                                                                                  <  2:517320/234‑1 (MQ=255)             ccaaCCATTCCCCGCTTCGATGACCTGCACCGACACTGAAGGCGCTCTACAGCCCGTGTGAGACCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGAGTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCTGTTCAAc                                                                                                                                                                                                                                  <  2:203059/248‑1 (MQ=255)                   aTTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAg                                                                                                                                                                                                                            <  1:506892/251‑1 (MQ=255)                       cacgcgccgcTTACCTGATCCGACACCGAAGGCGATCTACAGACCGTGTGAGATTCGCCATTTGGCGATTATCGTTAACCTGGAATATGCCCCGACAGCGGCGCTCCGCAATCAGGTGGTGCATTTCTATTTCTGGAAGCTGGAGGTCTTCAGCTTTGGCGAGAATCAAAAGTGCCTGTTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGGGCGGACGCTGCACGGCATCGGCGAGCAAt                                                                                                                                                                                                                        <  2:377344/244‑1 (MQ=255)                        ccGCTGCGCTTACATGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCCTCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATAAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGAGGAACGAAGTGCGGAAGCTGCACTTAAACGGCGAGCactc                                                                                                                                                                                                                       >  2:448245/1‑248 (MQ=255)                          gCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGAGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGGAGCGTTGACCTGCTGTACCGCAACTCGTGGAAAGAAGTGCGGACGCAGCCCTTCAACGACGAGCAATCg                                                                                                                                                                                                                      >  2:209454/1‑250 (MQ=255)                             gcgcTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCgat                                                                                                                                                                                                                    <  1:275285/249‑1 (MQ=255)                              cgcgtaccgtccccgacgcccaAGGCGCCCTACAGCCCTTGTGTGATCCGCCTACTGGGGAATACCGGTAACCTGTAATATGACCCGACAGCGGAGTTCAGCAATCAGGTGGCGCTTTTCGATTTCCGTAAGATTGATGGATTCAGCTTTGTCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTCCACTTCAACGGCGAGCAAGCgatgat                                                                                                                                                                                                                 <  2:305897/239‑1 (MQ=255)                                  tACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGa                                                                                                                                                                                                              <  1:389300/250‑1 (MQ=255)                                        ccccgacaccgaAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTAGCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTg                                                                                                                                                                                                       <  2:50691/249‑1 (MQ=255)                                              caccgaAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAAc                                                                                                                                                                                                  >  1:477072/1‑250 (MQ=255)                                                ccgaAGGCGCTCTACAGCCCGTGTGAGATCCGCCAGCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAAc                                                                                                                                                                                                  <  2:477072/248‑1 (MQ=255)                                                    aggctccctccAGCCCGTGTTAGATCCGCCACTTGGCGATTATCGTTAACCTGGAATATGCCCCGACAGCGGCGGTCCGCAAGCAGGTGGTGCATGTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATtctc                                                                                                                                                                                           <  2:249165/241‑1 (MQ=255)                                                              cAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGc                                                                                                                                                                                 >  1:456539/1‑251 (MQ=255)                                                               aGCCAGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGAACCGACAGCGGCGTTCCGCAATCAGGTAGTGCATTTCGATTTCCGGAAGCTGGATGTCTTCAGCTTTGGAGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGTGGAAGCTGCACTTCCACGGCGAGCAATCGGTGATCGAAGACCGGAACACTATTATCGGCAAAAtac                                                                                                                                                                                 >  2:160447/1‑248 (MQ=255)                                                                      gtgAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGAc                                                                                                                                                                         >  1:444549/1‑251 (MQ=255)                                                                       tgAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGAcg                                                                                                                                                                        >  1:242710/1‑251 (MQ=255)                                                                          gaTCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGAcgccg                                                                                                                                                                     >  1:387411/1‑251 (MQ=255)                                                                            cccGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCa                                                                                                                                                                   <  1:517662/250‑1 (MQ=255)                                                                              cGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCAcc                                                                                                                                                                 <  1:263319/251‑1 (MQ=255)                                                                                    cTGGCGATTATCGTTAACCTGGAATATGACCCGACATCGTCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTACGCTGGATGTCTTCAGATTTGGCGAGACTCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTTGAACGAAGTGCGGACGCTGCACTTCAACGGCGATCAATCGATGATCGCGACCCTGAACACTATTCCCGGCACAATGCATCAGCACCCCGCACCGcctg                                                                                                                                                            >  2:464444/1‑248 (MQ=255)                                                                                         gATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCg                                                                                                                                                      >  1:139332/1‑251 (MQ=255)                                                                                                     ccTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGtcttct                                                                                                                                          >  1:276333/1‑251 (MQ=255)                                                                                                        ggAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCACTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGACCCGCCAGATAGCGTGGAAGTCTTATGtt                                                                                                                                       >  2:428872/1‑251 (MQ=255)                                                                                                           atatGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTtata                                                                                                                                    <  1:26133/251‑1 (MQ=255)                                                                                                            tatGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTgg                                                                                                                                                   >  1:262521/1‑235 (MQ=255)                                                                                                            tatGACCCGACAGCGGCGTTCCGCAAGCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTgg                                                                                                                                                   <  2:262521/235‑1 (MQ=255)                                                                                                                 cccGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGc                                                                                                                              >  1:118980/1‑251 (MQ=255)                                                                                                                  ccGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGc                                                                                                                              >  1:567661/1‑250 (MQ=255)                                                                                                                    gACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGACGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGTTCGCGTGGACGTCTTCTGTTATAGCCAGCATc                                                                                                                           >  2:152626/1‑251 (MQ=255)                                                                                                                          ggctgtCCGCAATCAGGTGGTGTAGTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGcgcg                                                                                                                     <  2:415277/246‑1 (MQ=255)                                                                                                                            cgTTCCGCAATCAGGTGGTGCATTTCGATTTCCGAAAGCTGGATGTCTTCAGCTTTAGCGAGAGTCATAATTGCCTGGTAGGTAGCGTTGACATGCTGTACCGCAACTCGTGGAACGAAGTGCGGCCGCTGCACTTCAACGGCGGGCAATCGATGATCGAAGCCCTGACAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATGGCCAGCATCTGCGCgg                                                                                                                    >  2:532059/1‑250 (MQ=255)                                                                                                                                ccGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTaa                                                                                                               <  1:340825/251‑1 (MQ=255)                                                                                                                                ccGCAAACAGGGGGTGCATTTCGTTTTCCGGAATCTGGATGTCTGCAGACGTGGCGAGAATCAAAATTGCCTGGTAGGTAGTGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGGGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAACGCATCAGGACGCCGCACCGCCAGTAAGCGTGGAAGTCGGCTGTTATAGTCAGCATCTGCGCGGCTTaa                                                                                                               <  2:31855/251‑1 (MQ=255)                                                                                                                                 cGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAAt                                                                                                              <  1:107760/251‑1 (MQ=255)                                                                                                                                  gCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCCGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAACCTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGACGTCTTCTGTTATCGCCAGCATCTGCGCGGCTTAAtt                                                                                                             >  1:208550/1‑251 (MQ=255)                                                                                                                                  gCAATCAGATGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACATGCTGTACCGCAACTCGTGGAACGATGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGACAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTATGTTATAGCAGGCATATGCGCGGCTTAAtt                                                                                                             >  2:344119/1‑251 (MQ=255)                                                                                                                                         ggtggtGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTACTTCGTACTc                                                                                                      >  2:153278/1‑251 (MQ=255)                                                                                                                                          gtggtgCATTTCGATTTCCGTAAGCTGGATGTCTGCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTcg                                                                                                     <  2:36090/251‑1 (MQ=255)                                                                                                                                              tgCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCAAGCATCTGCGCGGCTTAATTCGTACTCGCGTg                                                                                                 >  1:532805/1‑251 (MQ=255)                                                                                                                                               gCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTg                                                                                                 <  2:81875/250‑1 (MQ=255)                                                                                                                                                 aTTTCGATTTCCGTAAGCTGTATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTTCCGCAACTCGTTGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTgcag                                                                                              <  2:444956/251‑1 (MQ=255)                                                                                                                                                 aTTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTgcag                                                                                              <  2:526432/251‑1 (MQ=255)                                                                                                                                                       aTTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCTGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTg                                                                                        <  1:212615/251‑1 (MQ=255)                                                                                                                                                            cGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAAAGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTc                                                                                   >  2:506151/1‑251 (MQ=255)                                                                                                                                                             gTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAATATTGCCTGGTAGGTAGCGTTGGCCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCt                                                                                  <  1:326068/251‑1 (MQ=255)                                                                                                                                                               aaGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGGACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTg                                                                                 >  2:45670/1‑250 (MQ=255)                                                                                                                                                                    ggATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAAtt                                                                                                             <  1:506057/217‑1 (MQ=255)                                                                                                                                                                    ggATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAAtt                                                                                                             >  2:506057/1‑217 (MQ=255)                                                                                                                                                                      aTGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTAtt                                                                         >  1:101117/1‑251 (MQ=255)                                                                                                                                                                                  ttGGCGAGAATCAAAATTTCTTGGGAGGGAGCGTGGACCTGCTGGACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCTCCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACGCGCGTGCAGCAACTGGTTTGTGAGTGTGTTGAATTGCGTc                                                               <  2:54898/249‑1 (MQ=255)                                                                                                                                                                                           aTCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCTTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCAcc                                                    <  1:263883/251‑1 (MQ=255)                                                                                                                                                                                                          taggtagCGTTGACCTGCTGTACCGCAACTCGAGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGCGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCggg                                     >  2:322863/1‑251 (MQ=255)                                                                                                                                                                                                                      tccTGCTGTACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCACCGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGCCGCCGCACCGCCAGATAGCGTGGAAGGCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCATCACCCGCCAGGAAACCGGGCGTTTTAAGgcg                         >  2:24869/2‑251 (MQ=255)                                                                                                                                                                                                                              tACCGCAACTCGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGt                 >  2:404981/1‑251 (MQ=255)                                                                                                                                                                                                                                        cGTGGAACGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGCGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCa        >  1:122926/1‑250 (MQ=255)                                                                                                                                                                                                                                             aaCGAAGTGCGGACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGCCGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGACTTTCCAGCACCCGCCAGGAAACCGGGCGTTTAAAGGCGCTGCGAGTTTCTGGTCAAACCTg  >  2:451903/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                TCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTG  >  NC_000913/3992723‑3993209

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