Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_vch_tpiA	2,214,152	A→T	Y165N (TAT→AAT)	yehU ←	inner membrane putative sensory kinase in two‑component system with YehT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_vch_tpiA	2,214,152	0	A	T	92.9%	129.8 / ‑5.0	42	Y165N (TAT→AAT)	yehU	inner membrane putative sensory kinase in two‑component system with YehT
Reads supporting (aligned to +/- strand):  ref base A (0/3);  new base T (16/23);  total (16/26)
Fisher's exact test for biased strand distribution p-value = 2.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGAAAAAGCCGAAGTGTATTTTTCAAACATCGCGCGTTTATCGAGCAATATACGCATAAACAGCGCCGCGCCGACGGTATTGGTGACCATCATTGGCGCAGCAATATTACTCACCAGACGCACCGCATCTTCATAAGGTCGGGCGATCGCAAGGATGATCAGCATTTGCACCATTTCAGCGACGAACGTGACGGCACCGGCGGTAATGGGGTTAAAGACTTTATCAGTGCGCCCGCGGCGGATCAGGATGCTGTGTACCAGGCCACCGAGTAATCCTTCA  >  NC_000913_3_vch_tpiA/2214017‑2214297                                                                                                                                        |                                                                                                                                                  cAGAAAAAGCCGAAGTACAAATGACTATCATCGCGCGTTTATCGAGCAATAGACGCATAAACAGCGCCGCGCCGACAGTATTGGTGACCAGCATTGGCGCAGCAATATTACTCACCAGACGCACCGCATCTTCATTAGGTCGGGCGATc                                                                                                                                      <  2:595730/149‑1 (MQ=255)     aaaaGCTGAAGGGTATTTTTCAAACATCGCGCGTTTATCGAGCAATATACGCATTAACAGCGCCGCGCCGACGGTATTGGTGACCATCATTGGCGCAGCAATATTACTCACCAGACGCACCGCATCTTCATTAGGTCGGGCGATCGCaa                                                                                                                                  <  1:726048/149‑1 (MQ=255)     aaaaGCCGAAGTGTATTTTTCAAACATCGCGCGTTTATCGAGCAATATACGCATAAACAGCGCCGCGCCGACGGTATTGGTGACCATCATTGGCGCAGCAATATTACTCACCAGACGCACCGCATCTTCATTAGGTCGGGCGATCGCaa                                                                                                                                  <  1:710432/149‑1 (MQ=255)     aaaaGCCGAAGTGTATTTTGCAAACATCGCGCGTTGATCGAGCAATATACGCATAAACAGCGCCGCGCCGACGGTATTGGTGACCATCATTGGCGCAGCAATATTACTCACCAGACGCACCGCATCTTCATTAGGTCGGGCGATCGCaa                                                                                                                                  <  1:510982/149‑1 (MQ=255)       aaGCCGAAGTGTATTTTTCAAACATCGCGCGTTTATCGAGCAATATACGCATAAACAGAGCCGCGCCGACGGTATTGGTGACCATCATTGGCGCAGCAATATTACTCACCAGACGCACCGCAGCTTCATTCGGTCGGGCGATCGCAAgg                                                                                                                                >  2:572022/1‑149 (MQ=255)        aGCCGAAGTGTAGTTTTAAAACATCGCGCGGTTATCGAGCAATATACGCATAAACACCGCCGCGACGACGGTATTGGTGACCATCATTGGCGCAGCAATATTACTCACCAGACGCACCGCATCTTCATTAGGTCGGGCGATCGCAAGga                                                                                                                               <  2:678259/149‑1 (MQ=255)          gcGAAGTGTATTTTTCAAACATCACGCGTTTATCGAGCAATATACGCATAAACAGCGCCGCGCCGACGGTACTGGTGACCATCATCGGCACAGCAATATTACTCACCAGACGCACCGCATCTTCATTAGGTCTGGAGATCGGAAgagcg                                                                                                                             >  2:778140/2‑145 (MQ=255)                tgtATTTTTCAAACAGCGCGCGTTTATCGAGCAATATACGCATAAACAGCGCCGCGCCGACGGTATTGGTGACCATCATTGGCGCAGCAATATTACTCACCAGACGCACCGCATCTTCATTAGGTCGGGCGATCGCAAGGATGATCAGc                                                                                                                       <  1:62192/149‑1 (MQ=255)                 gtATTTTTCAAACATCGCGCGTTTATCGAGCAATATACGCATAAACAGCGCCGCGCCGACGGTATTGGTGACCATCATTGGCGCAGCAATATTACTCACCAGACGTACCGCATCTTCATTAGGTCGGGCGATCGCAAGgatgat                                                                                                                           <  2:739161/144‑1 (MQ=255)                 gtATTTTTCAAACATCGCGCGTTTATCGAGCAATATACGCATAAACAGCGCCGCGCCGACGGTATTGGTGACCATCATTGGCGCAGCAATATTACTCACCAGACGCACCGCATCTTCATTAGGTCGGGCGATCGCAAGgatgat                                                                                                                           >  1:739161/1‑144 (MQ=255)                    tttttCAAACATCGCGCGTTTATCGAGCAATATACGCATAAACAGCGCCGCGCCGACGGTATTGGTGACCATCATTGGCGCAGCAATATTACTCACCAGACGCACCGCATCTTCATTAGGTCGGGCGATCGCAAGGATGATCAGCAttt                                                                                                                   <  2:610676/149‑1 (MQ=255)                        tCAAACATCGCGCGTTTATCGAGCAATATACGCATAAACAGCGCCGCGCCGACGGTATTGGTGACCATCATTGGCGCAGTAATATTACTCACCAGACGCACCGCATCTTCATTAGGTCGGGCGATCGCAAGGATGATCAGCATTTGCAc                                                                                                               <  1:538542/149‑1 (MQ=255)                          aaaCAGCGCGCGTTTATCGAGCAATATACGCATAAACAGCGCCGCGCCGACGGTATTGGTGACCATCATTGGCGCAGCAATATTACTCACCAGACGCACCGCATCTTCATTAGGTCGGGCGATCGCAAGGATGATCAGCATTTGCACCa                                                                                                             <  1:49421/149‑1 (MQ=255)                           aaCATCGCGCGTTTATCGAGCAATATACGCATAAACAGCGCCGCGCCGACGGTATTGGTGACCATCATTGGCGCAGCAATATTACTCACCAGACGCACCGCATCTTCATTAGGTCGGGCGATCGCAAGGATGATCAGCATTTGCACCAt                                                                                                            <  1:611181/149‑1 (MQ=255)                                 gcgcgTTTATCGAGCCATATACGCATAAACAGCGCCGCGCCGACGGTATTGGTGACCATCATTGGCGCAGCAATATTACTCACCAGACGCACCGCATCTTCATTAGGTCGGGCGATCGCAAGGATGATCAGCATTTGCACCATTTCAGc                                                                                                      >  2:522215/1‑149 (MQ=255)                                   gcgTTTATCGAGCAATATACGCATAAACAGCGCCGCGCCGACGGTATTGGTGACAATAATTGGCGCAGCAATATTACTCACCAGACGCACCGCATCCGCGTTAGGTCGGGCGATGGCAAGGATAGTCAGCATTTGCACCTGTTCAGcga                                                                                                    >  1:745307/1‑149 (MQ=255)                                                        cATAAATAGCGCCGCGCCGACGGTATTGGGTACCATCATGGGCGCAGCAATATTACTCACCAGACGCACCGCATCTTCATTAGGTCGGGCGATCGCAAGGATGATCAGCATTTGCACCATTTCAGCGACGAACGTGACGGCACcggcgg                                                                               <  2:525330/149‑1 (MQ=255)                                                                gcgccgcgccGACGGTATTGGTGACCATCATTGGCGCAGCAATATTACTCACCAGACGCACCGCATCTTCATTAGGTCGGGCGATCGCAAGGATGATCAGCATTTGCACCATTTCAGCGACGAACGTGACGGCACCGGCGGTAATgggg                                                                       <  1:339734/149‑1 (MQ=255)                                                                       cccgACGGAATTGGTGACAATCAATGGCGAAGCAATATTACTCACCAGACGCACCGCATCTTCATTAGGTCGGGCGAACGCAGGGATGATCAGCATTTGAACCATTTCAGCGACGAACGTGACGGCACCGGCGGTAATGGGGTTAAAGa                                                                <  2:21709/148‑1 (MQ=255)                                                                          gACGGTATTGGTGACCATCATTGGCGCAGCAATATTACTCACCAGACGCACCGCATCTTCATTAGGTCGGGCGATCGCAAGGATGATCAGCATTTGCACCATTTCAGCGACGAACGTGACGGCACCGGCGGTAATGGGGTTAAAGACtt                                                             >  1:70783/1‑149 (MQ=255)                                                                          gACGGTATTGGTGACCATCATTGGCGCAGCAATATTACTCACCAGACGCACCGACTATTCATTAGGTCGGGCGATCGCAAGGATGATCAGCATTTGCACCATTTCAGCGACGAACGTGACGGCCCCGACGGTAATGGGGTTAAAGACtt                                                             >  1:163037/1‑149 (MQ=255)                                                                             ggTATTGGTGACCATCATTGGCGCAGCAATATTACTCACCAGACGCACCGCATCTTCATTAGGTCGGGCGATCGCAAGGATGATCAGCATTTGCACCATTTCAGCGACGAACGTGACGGCACCGGCGGTAATGGGGTTAAAGACTTTAt                                                          <  1:592084/149‑1 (MQ=255)                                                                                 ttGGTGACCATCATTGGCTCAGCAATATTACTCACCAGACGCACCGCATCTTCATAATCTCCGGCGATCGCAAGGATGATCAGCATTTGCACCATTTCAGCGACGAACGTGACGGCACCGGCGGTAATGGGGTTAAAGACTTTATCAGt                                                      <  1:195109/149‑1 (MQ=255)                                                                                 ttGGTGACAATCATTGGCGCAGCAATATTACTCACCAGACGCACCGCATCTTTAAATTAATCTATGATCGCAAGGATGATCAGCATTTGCACCATTTCAGCGACGAACGTGACGGCACCGGCGGTAATGGGGTTAAAGACTTTATCAGt                                                      <  1:826715/149‑1 (MQ=255)                                                                                 ctGGTGACCATCCTTGGCGTAGCAATATTACTCACCAGACGCACCGCATCTTCATTAGGTCGGGTGAGCGCAAGGATGATCAGCATGTGCACCATTTCAGCGACGAACGTGACGGCACCGGCGGTAATGGGGTTAAAGACTTTATCAGt                                                      <  2:783680/148‑1 (MQ=255)                                                                                              ttGGCGCAGCAATATTACTCACCAGACGCACCGCATCTTCATTAGGTCGGGCGATCGCAAGGATGATCAGCATTTGCACCATTTCAGCGACGAACGTGACGGCGCCGGCGGTAATGGGGTTAGAGACTTTTTCTGTGAGCTCTgggcgg                                         >  1:45782/1‑149 (MQ=255)                                                                                                 gcgcAGCAATATTACTCACCAGACGCACCGCATCTTCATTAGGTCGGGCGATCGCAAGGATGATCAGCATTTGCACCATTTCAGCGACGAACGTGACGGCACCGGCGGTAATGGGGTTAAAGACTTTATCAGTGCGCCCGCGGCGGATc                                      <  1:753762/149‑1 (MQ=255)                                                                                                 gcgcAGCAATATTACTCACCAGACGCACCGCATCTTCATTAGGTCGGGCGATCGCAAGGATGATCAGCATTTGCACCATTTCAGCGACGAACGTGACGGCACCGGCGGTAATGGGGTTAAAGACTTTATCAGTGCGCCCGCGGCGGATc                                      <  2:70783/149‑1 (MQ=255)                                                                                                  cgcAGCAATATTACTCACCAGACGCACCGCATCTTCATTAGGTCGGGCGATCGCAAGGATGATCAGCATTTGCACCATTTCAGCGACGAACGTGACGGCACCGGCGGTAATGGGGTTAAAGACTTTATCAGTGCGCCCGCGGCGGATCa                                     <  1:769467/149‑1 (MQ=255)                                                                                                   gcagcaATATTACTCACCAGACGCACCGCATCTTCATTAGGTCGGGCGATCGCAAGGATGATCAGCATTTGCACCATTTCAGCGACGAACGTGACGGCACCGGCGGTAATGGGGTTAAAGACTTTATCAGTGCGCCCGCGGCGGATCAg                                    >  1:546716/1‑149 (MQ=255)                                                                                                     agtaaTATTCTTCACCACACGCACTGGATCTTCATTAGGTCGGGCGATCGCAAGGATGATCAGCATTTGCACCATTTCAGCGACGAACGTGACGGCACCGGCGGTAATGGGGTTAAAGACTTTATCAGTGCGCCCGCGGCGGATCAGGa                                  <  2:159433/146‑1 (MQ=255)                                                                                                                 cACCAGACGCACCGCATCTTCATTAGGTCGGGCGATCGCAAGGATGATCAGCATTTGCACCATTTCAGCGACGAACGTGACGGCACCGGCGGTAATGGGGTTAAAGACTTTATCAGTGCGCCCGCGGCGGATCAGGATGCTGTGTACCa                      <  1:522215/149‑1 (MQ=255)                                                                                                                          cACCGCATCTTCATTAGGTCGGGCGATCGCAAGGATGATCAGCATTTGCACCATTTCAGCGACGAACGTGACGGCACCGGCGGTAATGGGGTTAAAGACTTTATCAGTGCGCCCGCGGCGGATCAGGATGCTGTGTACCAGGCCACCGa             <  1:671709/149‑1 (MQ=255)                                                                                                                           aCCGCATCTTCATTAGGTCGGGCGATCGCAAGGATGATCAGCATTTGCACCATTTCAGCGACGAACGTGACGGCACCGGCGGTAATGGGGTTAAAGACTTTATCAGTGCGCCCGCGGCGGATCAGGATGCTGTGTACCAGGCCACCGAg            >  1:734873/1‑149 (MQ=255)                                                                                                                           aCCGCATCTTCATTAGGTCGGGCGATCGCAAGGATGATCAGCATTTGCACCATTTCAGCGACGAACGTGACGGCACCGGCGGTAATGGGGTTAAAGACTTTATCAGTGCGCCCGCGGCGGATCAGGATGCTGTGTACCAGGCCACCGAg            >  1:669656/1‑149 (MQ=255)                                                                                                                            ccGCATCTTCATTAGGTCGGGCGATCGCAAGGATGATCAGCATTTGCACCATTTCAGCGACGAACGTGACGGCACCGGCGGTAATGGGGTTAAAGACTTTATCAGTGCGCCCGCGGCGGATCAGGATGCTGTGTACCAGGCCACCGAGt           >  1:306394/1‑149 (MQ=255)                                                                                                                            ccGCATCTTCATTAGGTCGGGCGATCGCAAGGATGATCAGCATTTGCACCATTTCAGCGACGAACGTCACGCCACCGGCGGTAATGGGGTTAAAAACTTTATCAGTGCGCCCGCGGCGCACCAGGAGCCTGTGTACCAGGccccggact           >  2:153908/1‑142 (MQ=255)                                                                                                                               gaTCTTCATAGGGTCGAGCGATCGCAAGGATGATCAGCATTTGCACCATTGCAGCGACGAACGTGACGGCACCGGCGGCAATGGGGTTAAAGACTTGATCAGTGCGCCCGCGGCGGATCAGGATGCTGTGTACCAGGCCACCGAGTAAt        <  2:549100/148‑1 (MQ=255)                                                                                                                               cATCTTCATTAGGTCGGGCGATCGCAAGGATGATCAGCATTTGCACCATTTCAGCGACGAACGTGACGGCACCGGCGGTAATGGGGTTAAAGACTTTATCAGTGCGCCCGCGGCGGATCAGGATGCTGTGTACCAGGCCACCGAGTAAt        >  1:563481/1‑149 (MQ=255)                                                                                                                                   ttCATTAGGTCGGGCGATCGCAAGGATGATCAGCATTTg                                                                                                                  <  2:152118/39‑1 (MQ=255)                                                                                                                                   ttCATTAGGTCGGGCGATCGCAAGGATGATCAGCATTTg                                                                                                                  >  1:152118/1‑39 (MQ=255)                                                                                                                                     cATTAGGTCGGGCGATCGCAAGGATGATCAGCATTTGAACCATTTCAGCGACGAACGTGACGGCACCGGCGGTAATGGGGTTAAAGACTTTATCAGTGCGCCCGCGGCGGCACAGGATGCTGTGTACCAGGCAACCGCGTCATCCTTCa  >  1:255116/1‑149 (MQ=255)                                                                                                                                        |                                                                                                                                                  CAGAAAAAGCCGAAGTGTATTTTTCAAACATCGCGCGTTTATCGAGCAATATACGCATAAACAGCGCCGCGCCGACGGTATTGGTGACCATCATTGGCGCAGCAATATTACTCACCAGACGCACCGCATCTTCATAAGGTCGGGCGATCGCAAGGATGATCAGCATTTGCACCATTTCAGCGACGAACGTGACGGCACCGGCGGTAATGGGGTTAAAGACTTTATCAGTGCGCCCGCGGCGGATCAGGATGCTGTGTACCAGGCCACCGAGTAATCCTTCA  >  NC_000913_3_vch_tpiA/2214017‑2214297

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