Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,657,207	C→T	N71N (AAC→AAT)	hdeD →	acid‑resistance membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,657,207	0	C	T	100.0%	92.9 / NA	29	N71N (AAC→AAT)	hdeD	acid‑resistance membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/12);  total (17/12)

GGATTTTCTTCTTAATTTTATAGGGTGGTTCTATGTTATATATAGATAAGGCAACAATTTTGAAGTTTGATCTGGAGATGCTTAAAAAACATCGCAGAGCAATCCAGTTTATTGCCGTGCTGCTGTTTATCGTCGGGTTGCTGTGTATCAGTTTCCCGTTCGTCTCTGGCGATATTTTAAGCACAGTAGTGGGTGCATTATTAATCTGCTCGGGTATTGCGCTTATTGTCGGGTTATTCAGCAACCGCAGTCATAATTTCTGGCCGGTATTATCCGGTTTCCTCGTCGCAGTCGCCTATTTATTGATCGGCTATTTCTTCATCCGCGCACCGGAGCTGGGCATTTTTGCCATTGCGGCATTTATTGCCGGTTTGTTCTGTGTTGCAGGGGTTATTCGCCTGATGAGTTGGTATCGTCAGCGTTCAATGAAAGGCAGCTGGCTACAGCTCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGA  >  NC_000913/3656963‑3657445                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               ggATTTTCTTCTTAATTTTATAGGGAGGTTCTATGTTATATATAGATAAGGCAACAATTTTGAAGTTTGATCTGGAGATGCTTAAAAAACATCGCAGAGCAATCCAGTTTATTGCCGTGCTGCTGTTTATCGTCGGGTTGCTGTGAATCAGTTTCCCGTTCGTCTCTGGCGATATTTTAAGCACAGTAGTGGGTGCATTATTAATCTGCTCGGGTATTGCGCTTATTGTCGGGTTATTCAGCAATCGCAGt                                                                                                                                                                                                                                          >  1:223508/1‑251 (MQ=255)             tAATTTTATAGGGTGGTTCTATGTTATATATAGATAAGGCAACAATTTTGAAGTTTGATCTGGAGATGCTTAAAAAACATCGCAGAGCAATCCAGTTTATTGCCGTGCTGCTGTTTATCGTCGGGTTGCTGTGTATCAGTTTCCCGTTCGTCTCTGGCGATATTTTAAGCACAGTAGTGGGTGCATTATTAATCTGCTCGGGTATTGCGCTTATTGTCGGGTTATTCAGCAATCGCAGTCATAATTTCTg                                                                                                                                                                                                                               <  1:163860/250‑1 (MQ=255)                             ttCTATGTTATATATAGATAAGGCAACAATTTTGAAGTTTGATCTGGAGATGCTTAAAAAACATCGCAGAGCAATCCAGTTTATTGCCGTGCTGCTGTTTATCGTCGGGTTGCTGTGTATCAGTTTCCCGTTCGTCTCTGGCGATATTTTAAGCACAGTAGTGGGTGCATTATTAATCTGCTCGGGTATTGCGCTTATTGTCGGGTTATTCAGCAATCGCAGTCATAATTTCTGGCCGGTATTATc                                                                                                                                                                                                                   >  2:36869/1‑246 (MQ=255)                             ttCTATGTTATATATAGATAAGGCAACAATTTTGAAGTTTGATCTGGAGATGCTTAAAAAACATCGCAGAGCAATCCAGTTTATTGCCGTGCTGCTGTTTATCGTCGGGTTGCTGTGTATCAGTTTCCCGTTCGTCTCTGGCGATATTTTAAGCACAGTAGTGGGTGCATTATTAATCTGCTCGGGTATTGCGCTTATTGTCGGGTTATTCAGCAATCGCAGTCATAATTTCTGGCCGGTATTATc                                                                                                                                                                                                                   <  1:36869/246‑1 (MQ=255)                                                         aTTTTGAAGATTGATCTGGAGATGCTTCAAAAACATCGCAGAACAATCCAGTTTATTGCCGTGCTGCTGTTTATCGTCGGGTTGCTGTGTATCAGTTTCCCGTTCGTCTCTGGCGATATTTTAAGCACAGTAGTGGGTGCATTATTAATCTGCTCGGGTATTGCGCTTATTGTCGGGTTATTCAGCAATCGCAGTCATAATTTCTGGCCGGTATTATCCGGTTTCCTCATCGCAGTCGCCTATTTATTGAt                                                                                                                                                                                  >  1:240569/1‑251 (MQ=255)                                                                    tGATCTGGAGATGCTTAAAAAACATCGCAGAGCAATCCAGTTTATTGCCGTGCTGCTGTTTATCGTCGGGTTGCTGTGTATCAGTTTCCCGTTCGTCTCTGGCGATATTTTAAGCACAGTAGTGGGTGCATTATTAATCTGCTCGGGTATTGCGCTTATTGTCGGGTTATTCAGCAATCGCAGTCATAATTTCTGGCCGGTATTATCCGGTTTCCTCGTCGCAGTCGCCTATTTATTGATCGGCTATttct                                                                                                                                                                       >  2:204707/1‑251 (MQ=255)                                                                    tGATCTGGAGATGCTTAAAAAACATCGCAGAGCAATCCAGTTTATTGCCGAGCTGCTGTTTATCGTCGGGCTGCTGGGTATCAGATTCCCGTTCGTCTCTGGCGCTATTTTAAGCACAGTAGTGGGTGCATTATTAATCTGTTCGGGCATTGCGCTTATTGTCGGGTTATTTAGCACTCGCAGTCATAATTTCTGGCCGGTATTATCCGGTTTCCTCTTCGCAGCCGCCTATTTATTGATCGGTTATttct                                                                                                                                                                       >  2:332202/1‑251 (MQ=255)                                                                                             cGCAGAGCAATCCAGTTTATTGCCGTGCTGCTGTTTATCGTCGGGTTGCTGTGTATCAGTTTCCCGTTCGTCTCTGGCGATATTTTAAGCACAGTAGTGGGTGCATTATTAATCTGCTCGGGTATTGCGCTTATTGTCGGGTTATTCAGCAATCGCAGTCATAATTTCTGGCCGGTATTATCCGGTTTCCTCGTCGCAGTCGCCTATTTCTTGATCGGCTATTTCTTCATCCGCGCACCGGAGCTGGNCAt                                                                                                                                              >  2:364280/1‑251 (MQ=255)                                                                                               cAGAGCAATCCAGTTTATTGCCGTGCTGCTGTTTATCGTCGGGTTGCTGTGTATCAGTTTCCCGTTCGTCTCTGGCGATATTTTAAGCACAGTAGTGGGTGCATTATTAATCTGCTCGGGTATTGCGCTTATTGTCGGGTTATTCAGCAATCGCAGTCATAATTTCTGGCCGGTATTATCCGGTTTCCTCGTCGCAGTCTCCTATTTATTGATCGGCTATTTCTTCATCCGCGCACCGGAGCTGGGCAttt                                                                                                                                            >  2:107417/1‑251 (MQ=255)                                                                                                agagCAATCCAGTTTCTTGCCGTGCTGCTGTTTATCGTCGGGTTGCTGTGTATCAGTTTCCCGTTCGTCTCTGGCGATATTTTAAGCACAGTAGTGGGTGCATTATTAATCTGCTCGGGTATTGCGCTTATTGTCGGGTTATTCAGCAATCGCAGTCATAATTTCTGGCCGGTATTATCCGGTTTCCTCGTCGCAGTCGCCTATTTATTGATCGGCTATTTCTTCATCCGCGCACCGGAGCTGGGCAtttt                                                                                                                                           >  1:77192/1‑251 (MQ=255)                                                                                                          ggTTTATTGCCGTGCTGCTGTTTATCGTCGGGTTGGTGTGTATCAGTTTCCCGTTCGTCTCTGGCGATATTTTAAGCACAGTAGTGGGTGCATTATTAATCTGCTCGGGTATTGCGCTTATTGTCGGGTTATTCAGCAATCGCAGTCATAATTTCTGGCCGGTATTATCCGGTTTCCTCGTCGCAGTCGCCTATTTATTGATCGGCTATTTCTTCATCCGCGCACCGGAGCTGGGCATTTTTGCCATTGCg                                                                                                                                 >  2:24267/2‑251 (MQ=255)                                                                                                           gTTTATTGCCGTGCTGCTGTTTATCGTCGGGTTGCTGTGTATCAGTTTCCCGTTCGTCTCTGGCGATATTTTAAGCACAGTAGTGGGTGCATTATTAATCTGCTCGGGTATTGCGCTTATTGTCGGGTTATTCAGCAATCGCAGTCATAATTTCTGGCCGGTATTATCCGGTTTCCTCGTCGCAGTCGCCTATTTATTGATCGGCTATTTCTTCATCCGCGCACCGGCGCTGGGCATTTTTGCCATTGCgg                                                                                                                                >  1:22746/1‑251 (MQ=255)                                                                                                                                 gtcgtcgGGTTGCTGTGTATCAGTTTCCCGTTCGTCTCTGGCGGTATTTTAAGCACAGTAGTGGGTGCATTAGTAATCTGCTCGGGTATTGCGCTTATTGTCGGGTTATTCAGCAATCGCAGTCATAATTTCTGGCCGGTATTATCCGGTTTCCTCGTCGCAGTCGCCTATTTATTGATCGGCTATTTCTTCGTCCGCGCACCGGAGCTGGGCAGTTTTGCCATTGCGGCATTGATTGCCGGTTCGTTCtg                                                                                                          <  2:378295/250‑1 (MQ=255)                                                                                                                                 aTCGTCGGGTTGCTGTGTATCAGTTTCCCGTTCGTCTCTGGCGATATTTTAAGCACAGTAGTGGGTGCATTATTAATCTGCTCGGGTATTGCGCTTATTGTCGGGTTATTCAGCAATCGCAGTCATAATTTCTGGCCGGTATTATCCGGTTTCCTCGTCGCAGTCGCCTATTTATTGATCGGCTATTTCTTCATCCGCGCACCGGAGCTGGGCATTTTTGCCATTGCGGCATTTATTGCCGGTTTGTTCtg                                                                                                          <  2:4216/251‑1 (MQ=255)                                                                                                                                    gtcgGGTTGCTGTGTATCAGTTTCCCGTTCGTCTCTGGCGATATTTTAAGCACAGTAGTGGGTGCATTATTAATCTGCTCGGGTATTGCGCTTATTGTCGGGTTATTCAGCAATCGCAGTCATAATTTCTGGCCGGTATTATCCGGTTTCCTCGTCGCAGTCGCCTATTTATTGATCGGCTATTTCTTCATCCGCGCACCGGAGCTGGGCATTTTTGCCATTGCGGCATTTATTGCCGGTTTGTTCTgtgt                                                                                                       >  1:374157/1‑251 (MQ=255)                                                                                                                                    gtcgGGATGCTGTGTATCAGTTTCCCGTTCGTCTCTGGCGATCTTTTAAGCACAGTAGTGGGTGCATTATTAATCTGCTCGGGTATTGCGCTTATTGTCGGGTTATTAAGCAATCGCAGTCATAATTTCTGGCCGGTATTATCCGGTTTCCTCGTCGCAGTCGCCTATTTATTGATCGGCTATTTCTTCATCCGCGCACCGGAGCTGGGCATTTTTGCCATTGCGGCATTTATTGCCGGTTTGTTCTgtgt                                                                                                       >  2:161348/1‑251 (MQ=255)                                                                                                                                            gCTGTGTATCAGTTTCCCGTTCGTCTCTGGCGATATTTTAAGCACAGTAGTGGGTGCATTATTAATCTGCTCGGGTATTGCGCTTATTGTCGGGTTATTCAGCAATCGCAGTCATAATTTCTGGCCGGTATTATCCGGTTTCCTCGTCGCAGTCGCCTATTTATTGATCGGCTATTTCTTCATCCGCGCACCGGAGCTGGGCATTTTTGCCATTGCGGCATTTATTGCCGGTTTGTTCTGTGTTGCAgtg                                                                                                <  2:402425/250‑3 (MQ=255)                                                                                                                                            gCTGTGTATCAGTTTCCCGTTCGTCTCTGGCGATATTTTAAGCACAGTAGTGGGTGCATTATTAATCTGCTCGGGTATTGCGCTTATTGTCGGGTTATTCAGCAATCGCAGTCATAATTTCTGGCCGGTATTATCCGGTTTCCTCGTCGCAGTCGCCTATTTATTGATCGGCTATTTCTTCATCCGCGCACCGGAGCTGGGCATTTTTGCCATTGCGGCATTTATTGCCGGTTTGTTCTGTGTTGCAgggg                                                                                               <  1:261005/251‑1 (MQ=255)                                                                                                                                                 gtATCAGTTTCCCGTTCGTCTCTGGCGATATTTTAAGCACAGTAGTGGGTGCATTATTAATCTGCTCGGGTATTGCGCTTATTGTCGGGTTATTCAGCAATCGCAGTCATAATTTCTGGCCGGTATTATCCGGTTTCCTCGTCGCAGTCGCCTATTTATTGATCGGCTATTTCTTCATCCGCGCACCGGAGCTGGGCATTTTTGCCATTGCGGCATTTATTGCCGGTTTGTTCTGTGTTGCAGGGGTTAt                                                                                           <  1:364280/250‑1 (MQ=255)                                                                                                                                                          tCCCGTTCGTCTCTGGCGATATTTTAAGCACAGTAGTGTGTGCATTATTAATCTGCTCGGGTATTGCGCTTATTGTCGGGTTATTCAGCAATCGCAGTCATAATTTCTGGCCGGTATTATCCGGTTTCCTCGTCGCAGTCGCCTATTTATTGATCGGCTATTTCTTCATCCGCGCACCGGAGCTGGGCATTTTTGCCATTGCGGCATTTATTGCCGGTTTGTTCTGTGTTGCAGGGGTTATTCGCCtgat                                                                                  <  2:195579/250‑1 (MQ=255)                                                                                                                                                          tCCCGTTCGTCTCTGGCGATATTTTAAGCACAGTAGTGGGTGCATTATTAATCTGCTCGGGTATTGCGCTTATTGTCGGGTTATTCAGCAATCGCAGTCATAATTTCTGGCCGGTATTATCCGGTTTCCTCGTCGCAGTCGCCTATTTATTGATCGGCTATTTCTTCATCCGCGCAACGGAGCTGGGCATTTTTGCCATTGCGGCATTTATTGCCGGTTTGTTCTGTGTTGCAGGGGTTATTCGCCtgatg                                                                                 >  1:44395/1‑251 (MQ=255)                                                                                                                                                                        ggCGATATTTAAAGCACAGTAGTGGGTGCATTATTAATCGTCGCGGGTATTGCGCTTATTGTCGTGTTATTCAGCAATCGCAGTCATAATTTCTGGCCGGTATTAACCGGTTTCCTCGTCGCAGTCGCCTATTTATTGATCGGCTATTTCAGCATCCGCGCACCGGAGCTGTGCATTTTTGCCATTGCGTCATTTATTGCCGGTTTGTTCTGTGTTGCAGGGGTTATTCGCCTTATTAGTTGGTATCGTCa                                                                   <  2:204877/251‑1 (MQ=255)                                                                                                                                                                                   aaGCACCGTAGTGGGTGCATTATTAATCTGCTCGGGTATTGCGCTTATTGTCGGGTTATTCAGCAATCGCAGTCATAATTTCTGGCCGGTATTATCCGGTTTCCTCGTCGCAGTCGCCTATTTATTGATCGGCTATTTCTTCATCCGCGCACCGGAGCTGGGCATTTTTGCCATTGCGGCATTTATTGCCGGTTTGTTCTGTGTTGCAGGGGTTATTCGCCTGATGAGTTGGTATCGTCAGCGTTCAATg                                                         >  1:301907/1‑250 (MQ=255)                                                                                                                                                                                                                           gcgcTTATTGTCGGGTTATTCAGCAATCGCAGTCATAATTTCTGGCCGGTATTATCCGGTTTCCTCGTCGCAGTCGCCTATTTATTGATCGGCTATTTCTTCATCCGCGCACCGGAGCTGGGCATTTTTGCCATTGCGGCATTTATTGCCGGTTTGTTCTGTGTTGCAGGGGTTATTCGCCTGATGAGTTGGTATCGTCAGCGTTCAATGAAAGGCAGCTGGCTACAGCTCGTTATTGGTGTGCTGGatat                >  1:178094/1‑251 (MQ=255)                                                                                                                                                                                                                                    gTCGGGTTATTCAGCAATCGCAGTCATAATTTCTGGCCGGTATTATCCGGTTTCCTCGTCGCAGTCGCCTATTTATTGATCGGCTATTTCTTCATCCGCGCACCGGAGCTGGGCATTTTTGCCATTGCGGCATTTATTGCCGGTTTGTTCTGTGTTGCAGGGGTTATTCGCCTGATGAGTTGGTATCGTCAGCGTTCAATGAAAGGCAGCTGGCTTCAGCTCGTTATTGGTGTGCTGGATATCGTCATTGc       <  1:10149/251‑1 (MQ=255)                                                                                                                                                                                                                                       gggTTATTCAGCAATCGCAGTCATAATTTCTGGCCGGTATTATCCGGTTTCCTCGTCGCAGTCGCCTATTTATTGATCGGCTATTTCTTCATCCGCGCACCGGAGCTGGGCATTTTTGCCATTGCGGCATTTATTGCCGGTTTGTTCTGTGTTGCAGGGGTTATTCGCCTGATGAGTTGGTATCGTCAGCGTTCAATGAAAGGCAGCTGGCTACAGCTCGTTATTGGTGTGCTGGATATCGTCATTGCCTg    >  1:187259/1‑251 (MQ=255)                                                                                                                                                                                                                                         gTTATTCAGCAATCGCAGTCATAATTTCTGGCCGGTATTATCCGGTTTCCTCGTCGCAGTCGCCTATTTATTGATCGGCTATTTCTTCATCCGCGCACCGGAGCTGGGCATTTTTGCCATTGCGGCATTTATTGCCGGTTTGTTCTGTGTTGCAGGGGTTATTCGCCTGATGAGTTGGTATCGTCAGCGTTCAATGAAAGGCAGCTGGCTACAGCTCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGa  <  2:279427/251‑1 (MQ=255)                                                                                                                                                                                                                                             ttCAGCAATCGCAGTCATAATTTCTGGCCGGTATTATCCGGTTTCCTCGTCGCAGTCGCCTATTTATTGATCGGCTATTTCTTCATCCGCGCACCGGAGCTGGGCATTTTTGCCATTGCGGCATTTATTGCCGGTTTGTTCTGTGTTGCAGGGGTTATTCGCCTGATGAGTTGGTATCGTCAGCGTTCAATGAACGGCAGCTGGCTACAGCTCGTTATTGGTGTGCTGGATATCGTCAtt         >  2:384777/1‑240 (MQ=255)                                                                                                                                                                                                                                             ttCAGCAATCGCAGTCATAATTTCTGGCCGGTATTATCCGGTTTCCTCGTCGCAGTCGCCTATTTATTGATCGGCTATTTCTTCATCCGCGCACCGGAGCTGGGCATTTTTGCCATTGCGGCATTTATTGCCGGTTTGTTCTGTGTTGCAGGGGTTATTCGCCTGATGAGTTGGTATCGTCAGCGTTCAATGAAAGGCAGCTGGCTACAGCTCGTTATTGGTGTGCTGGATATCGTCAtt         <  1:384777/240‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               GGATTTTCTTCTTAATTTTATAGGGTGGTTCTATGTTATATATAGATAAGGCAACAATTTTGAAGTTTGATCTGGAGATGCTTAAAAAACATCGCAGAGCAATCCAGTTTATTGCCGTGCTGCTGTTTATCGTCGGGTTGCTGTGTATCAGTTTCCCGTTCGTCTCTGGCGATATTTTAAGCACAGTAGTGGGTGCATTATTAATCTGCTCGGGTATTGCGCTTATTGTCGGGTTATTCAGCAACCGCAGTCATAATTTCTGGCCGGTATTATCCGGTTTCCTCGTCGCAGTCGCCTATTTATTGATCGGCTATTTCTTCATCCGCGCACCGGAGCTGGGCATTTTTGCCATTGCGGCATTTATTGCCGGTTTGTTCTGTGTTGCAGGGGTTATTCGCCTGATGAGTTGGTATCGTCAGCGTTCAATGAAAGGCAGCTGGCTACAGCTCGTTATTGGTGTGCTGGATATCGTCATTGCCTGGA  >  NC_000913/3656963‑3657445

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