Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,776,715	A→T	S377C (AGT→TGT)	ydiF →	putative acetate‑CoA transferase YdiF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,776,715	0	A	T	98.7%	223.7 / NA	75	S377C (AGT→TGT)	ydiF	putative acetate‑CoA transferase YdiF
Reads supporting (aligned to +/- strand):  ref base A (0/1);  major base T (33/40);  minor base G (1/0);  total (34/41)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GCATTATTCGAAATGCGTAAAGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCAGGAAGGACGGGTGAAGAAATTTATTCGGGAACTACCGGAAATTACTTTCAGCGGAAAA  >  NC_000913/1776469‑1776951                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             gcATTATTCGAAATGCGTAAAGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTGCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGttt                                                                                                                                                                                                                                          <  1:100919/251‑1 (MQ=255) gcATTATTCGAAATGCGTAAAGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGttt                                                                                                                                                                                                                                          >  1:63250/1‑251 (MQ=255) gcATTATTCGAAATGCGTAAAGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGttt                                                                                                                                                                                                                                          <  1:737203/251‑1 (MQ=255) gcATTATTCGAAATGCGTAAAGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGttt                                                                                                                                                                                                                                          >  1:692749/1‑251 (MQ=255)     tattCGAAATGCGTAAAGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGCACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCt                                                                                                                                                                                                                                      >  1:948594/1‑251 (MQ=255)            aaTGCGTAAAGTCGCGGGGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGTGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCg                                                                                                                                                                                                                               <  2:257494/251‑1 (MQ=255)              tGCGTAAAGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGAc                                                                                                                                                                                                                             >  1:897494/1‑251 (MQ=255)               gCGTAAAGGCGCGGTGGGGAATGTCGGCGGCGGTATTTCTGACGGCAATGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTAGTCTGACGGTAGAAACATGTCCGATTGTCGGAAGTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCTCATTTTGATTTTTATCACGGTGGCG‑GTCTGGATGGTTGGTATGTGTGTTTTGCTGAAGTCGAcc                                                                                                                                                                                                                            <  2:159180/251‑1 (MQ=255)               gCGTAAAGGCGCAGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GCCTGGATGTTTGTTCTTTGTGTTTTGCTGAAGTCGAcc                                                                                                                                                                                                                            >  1:223775/1‑251 (MQ=255)                  tAAAGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGc                                                                                                                                                                                                                         >  1:171057/1‑251 (MQ=255)                     aGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGTTCGCCTTTGGAGCGAACGTGAATACCCGTGCCATTTTGGGTAATACGTACCAGTTTGATTTTTATCACGGGGGCG‑GTCTGGATGATAGGTATTTGGGTTTTGCTGAAGttagccagcacg                                                                                                                                                                                                                      >  2:402577/1‑240 (MQ=255)                      ggCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGAGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GACTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACgg                                                                                                                                                                                                                     >  1:466050/1‑251 (MQ=255)                        cgcgGGGGGGAATGTCGGCGTCGGTATTGCTGACGGCCTTGGCCTGGTCGCCCGAGAAGCAGGTTGTGATGATGACTTTATTCTGCAGGTAGAACCAGGTCCGATTGGCGGAATTACTTCACAGGGGCTCCCCTTTGGCGCGAACGTGAATAACCGTGCCCTTCTGGATATGACGTCCCCGTTTGATTTTTATCACGG‑GGCGCGGCTGGATGTTTGTTTGTTGTGTTTTTCTGGAGACGACCCGCACgact                                                                                                                                                                                                                   >  2:748132/1‑248 (MQ=255)                                          cGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAAt                                                                                                                                                                                                 <  1:801036/251‑1 (MQ=255)                                              tgTATTGCTGACGGCATTGGCCTGTTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGGTTGGCGGAATTACTTCACAGGGGGTCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATAGGACGTACCAGATTGATTTTTATCAAGGAGGCG‑GTCTGGTTGTTTGTTCTTTGTTTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCCTAAATTCa                                                                                                                                                                                              >  2:550937/2‑250 (MQ=255)                                                 aTTGCTGACGGCATTGGCCTGGTCGCTCGAGAAGAAGGTTGGGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGTAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATgg                                                                                                                                                                                          <  1:346935/251‑1 (MQ=255)                                                               tGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGc                                                                                                                                                                             <  1:16512/250‑1 (MQ=255)                                                               gggCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATCCTGACGGTAGAAACAGGTCCGGTCGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGGGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCa                                                                                                                                                                            <  2:103192/250‑1 (MQ=255)                                                                 gCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCAcc                                                                                                                                                                          <  1:260282/251‑1 (MQ=255)                                                                      gTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCCTAAATTCAATGGTAAAATCATGGGCACCGGTgg                                                                                                                                                                     >  2:501417/1‑251 (MQ=255)                                                                        cGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATGCCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTGTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGa                                                                                                                                                                    <  2:119772/250‑1 (MQ=255)                                                                             gagaAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTAtt                                                                                                                                                              <  1:292863/251‑1 (MQ=255)                                                                              agaagaAGGTTGTGCGGATGGCTTTAGTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTg                                                                                                                                                             <  2:499515/251‑1 (MQ=255)                                                                               gaagaaGGTTGTGCGGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGa                                                                                                                                                            <  2:736946/251‑1 (MQ=255)                                                                                    agggtgtgTTGATGGCTTTATGCTGACGGTAGAAACAGGTTCGATTGGCGGAATTACTTCACAGGGGATAGGCTTTGGCGCGAACGTGAATACCCGTGCCATTCTTGATATGACGTTCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGataaca                                                                                                                                                       <  2:688668/247‑4 (MQ=255)                                                                                    agGTTGTGTTGATGACTTTATTCTGACGGTCGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATAGCCTTTGGCTCGAACGTGAATACCCGTGCCATTCTGGCTATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGGTGTTTGGTATTTGTGTTTTGCTTAATACGACCAGCACGGTAACATGGGCGAGCATAAATTTCATGGTAAAATCATGGGCAACGGGGGAGTTATTGatta                                                                                                                                                         >  2:794138/1‑247 (MQ=255)                                                                                    agGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATTCCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCa                                                                                                                                                       <  1:287007/251‑1 (MQ=255)                                                                                      gtgggtcTGATGACTTGATTCTGACGGGAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCTCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATAGCAGt                                                                                                                                                     <  2:69293/245‑1 (MQ=255)                                                                                      gTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGt                                                                                                                                                     <  1:501417/251‑1 (MQ=255)                                                                                      gTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGt                                                                                                                                                     >  1:144487/1‑251 (MQ=255)                                                                                        tgtgCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGc                                                                                                                                                   >  1:605093/1‑251 (MQ=255)                                                                                        tgtgATGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGGATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTTTTCTTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTGTTTATATCAGTGc                                                                                                                                                   >  2:20482/1‑251 (MQ=255)                                                                                          tgCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGcc                                                                                                                                                  >  1:340817/1‑250 (MQ=255)                                                                                           gCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCAc                                                                                                                                                <  1:2595/251‑1 (MQ=255)                                                                                               atgaATTTATTCTGACGGTAGAAACAGGTTCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTc                                                                                                                                             <  1:577585/250‑1 (MQ=255)                                                                                                  aCTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGCAAGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGGTTGATTTTTATCCCGGGGGCG‑GGCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGAATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTTATATCAGTGCAACttata                                                                                                                                           >  2:346444/1‑246 (MQ=255)                                                                                                     ttATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATGTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCgaagaa                                                                                                                                       <  2:948594/250‑1 (MQ=255)                                                                                                                gTAGAAACAGGTTCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCTGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTCTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTg                                                                                                                           <  2:948445/251‑1 (MQ=255)                                                                                                                    aaaCAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGc                                                                                                                       <  1:51017/251‑1 (MQ=255)                                                                                                                        aGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACAt                                                                                                                   <  1:275188/251‑1 (MQ=255)                                                                                                                         ggTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACAtt                                                                                                                  <  1:244710/251‑1 (MQ=255)                                                                                                                          gTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCCCCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACAtt                                                                                                                  >  2:813124/1‑250 (MQ=255)                                                                                                                               aTTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATTCCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGc                                                                                                            <  2:229263/251‑1 (MQ=255)                                                                                                                                     ggAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGc                                                                                                        >  1:847414/1‑249 (MQ=255)                                                                                                                                        aTTACTTCACAGGGGATCGCCCTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCTCTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGttt                                                                                                   <  2:249968/251‑1 (MQ=255)                                                                                                                                           aCTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTaaa                                                                                                <  1:813124/251‑1 (MQ=255)                                                                                                                                             ttCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAACATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTaaaa                                                                                               >  2:792662/1‑250 (MQ=255)                                                                                                                                                acaGGGGATCGCCTTTGGCGCGAACGTGCATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAAc                                                                                             >  2:464110/1‑249 (MQ=255)                                                                                                                                                 caGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGa                                                                                          <  1:502584/251‑1 (MQ=255)                                                                                                                                                         ggCCTTTGTCGCGAACGTGAATGCCCGTGCCATTCTGGAGATGACGTCCCAGGTTGATTTTTATCACGGGGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGGCGGCGGGCATAAAGTCAATGGTAAAATCATGGGCAGCGGTTGATTTATTGATATCAGTGCCACTTCGAAGAAAGTCATTGTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCg                                                                                  <  1:398869/250‑1 (MQ=255)                                                                                                                                                         cGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAATAACAGAAATTACCg                                                                                  >  2:207805/1‑251 (MQ=255)                                                                                                                                                           ccTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGTTGTTTGTTCTTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGCTTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTTCGGGCAGTTTAAACACAGAAATTACCGAc                                                                                >  2:139421/1‑251 (MQ=255)                                                                                                                                                           ccTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGAc                                                                                >  2:841551/1‑251 (MQ=255)                                                                                                                                                           ccTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGAc                                                                                >  2:79982/1‑251 (MQ=255)                                                                                                                                                            cTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACg                                                                               >  1:57949/1‑251 (MQ=255)                                                                                                                                                              ttGGCGCGAACGTGAATGCCCGGGCCATTTTGGATAGTACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGGCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGTCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACg                                                                               <  2:650957/249‑1 (MQ=255)                                                                                                                                                                ggCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGTAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGCAAAATAATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCAAAGAAAATCCTTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCa                                                                            >  2:718658/1‑250 (MQ=255)                                                                                                                                                                  cgcgAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAAt                                                                         <  1:691408/251‑1 (MQ=255)                                                                                                                                                                    cgAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATAATTTTCTGCGGCACATTAACTGCGGGCAGTTGAAAAACAAAAATTACCGAAGGCaacat                                                                        >  2:279198/1‑247 (MQ=255)                                                                                                                                                                     gAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTa                                                                       >  1:361880/1‑250 (MQ=255)                                                                                                                                                                            aaTACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAAt                                                                                                                                                                                            >  1:163802/1‑126 (MQ=255)                                                                                                                                                                                     ggcatgcTGGAGATGACGTCCCAGTTTGATTTTTAGCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTGTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTAGTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCaggaag                                                       <  2:374642/244‑1 (MQ=255)                                                                                                                                                                                                   gACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCAGGAAGGACGGGTGAAGaaa                                         <  1:464110/250‑1 (MQ=255)                                                                                                                                                                                                   gACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCAGGAAGGACGGGTGAAGaaa                                         <  2:306746/250‑1 (MQ=255)                                                                                                                                                                                                        cccAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGGTTGTTATTTGAGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCAGGAAGGACGGGTGAAGAAATTTAt                                    <  1:261219/250‑1 (MQ=255)                                                                                                                                                                                                               tggtttttATCTCGGTGGTG‑GTCTGTATGTTTGTGATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGGCGGCGTGCATAAATGCAATGGTAAAATCATGGGCACTTGTGGATGTAGTGAGAGCAGTGCCACTTCGAAGAAAATCATTTTCGGCGGCCCATTAACTGCGGGCAGTTTAAAAGCAGAAATTTCCGACGGCAAAGTAAATATCGGCCAGGAAGGCCGGGTGAAGAAATTTATTCggtaa                             <  2:36808/247‑4 (MQ=255)                                                                                                                                                                                                                 aTTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCAGGAAGGACGGGTGAAGAAATTTATTCGGGAACTa                          <  2:801312/251‑1 (MQ=255)                                                                                                                                                                                                                  tttttATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCAGGAAGGACGGGTGAAGAAATTTATTCGGGAACTAc                         >  1:13561/1‑251 (MQ=255)                                                                                                                                                                                                                    tttATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCAGGAAGGACGGGTGAAGAAATTTATTCGGGAACTACCg                       >  1:779571/1‑251 (MQ=255)                                                                                                                                                                                                                         cACGGTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCAGGAAGGACGGGTGAAGAAATTTATTCGGGAACTACCGGaaa                   <  2:63250/250‑1 (MQ=255)                                                                                                                                                                                                                            ggTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCAGGAAGGACGGGTGAAGAAATTTATTCGGGAACTACCGGAAATTAc               >  1:617535/1‑251 (MQ=255)                                                                                                                                                                                                                            ggTGGCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCAGGAAGGACGGGTGAAGAAATTTATTCGGGAACTACCGGAAATTAc               <  2:897494/251‑1 (MQ=255)                                                                                                                                                                                                                                gCG‑GTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCAGGAAGGACGGGTGAAGAAATTTATTCGGGAACTACCGGAAATTACTTTc           <  1:222050/251‑1 (MQ=255)                                                                                                                                                                                                                                        ggATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCAGGAAGGACGGGTGAAGAAATTTATTCGGGAACTACCGGAAATTACTTTCAGCGGa     >  1:764205/1‑250 (MQ=255)                                                                                                                                                                                                                                          aTGTTTGTTATGTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATGCAATGGTAAAATCATGGGCACCGTTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCAGGAAGGACGGGTGAAGAAATTTATTCGGGAACTACCGGAAATTACTTTCAGCGGaaaa  <  1:286514/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             GCATTATTCGAAATGCGTAAAGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCG‑GTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCAGGAAGGACGGGTGAAGAAATTTATTCGGGAACTACCGGAAATTACTTTCAGCGGAAAA  >  NC_000913/1776469‑1776951

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: ‑