Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,380,208	G→A	13.3%	A6A (GCG→GCA)	yhcB →	conserved protein YhcB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,380,208	0	G	A	13.3%	65.7 / 4.2	30	A6A (GCG→GCA)	yhcB	conserved protein YhcB
Reads supporting (aligned to +/- strand):  ref base G (15/11);  new base A (2/2);  total (17/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

GCGATATGCAGTCGGGATTCCACTTCTGTGGAATTAACGGTTATAGTGGCATAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCGCTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGCAACAGGCGTTGCAGTACGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTA  >  NC_000913/3380073‑3380350                                                                                                                                        |                                                                                                                                               gcaatatGCAGTCGGGATTCGACGTCTGTGGAATTAACGGGTATAGTGGCATAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCGACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCGCTAATTg                                                                                                                                         <  1:105187/140‑1 (MQ=255)    atatGCAGTCGGGATTCCGCTTCTGTGGAATTAACGGTTATAGGGGCATAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCGCTAATTGGGt                                                                                                                                      <  1:38957/143‑1 (MQ=255)        gCAGTCGGGATTCCACTTCTGTGGAATTAACGGTTATAGTGGCATAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCGCTAATTGGGTTAgt                                                                                                                                  <  1:184969/143‑1 (MQ=255)          aGTCGGGATTCCACTTCTGTGGAATTAACGGTTATAGTGGCATAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCGCTAATCGGGTTAgtcg                                                                                                                                >  2:41732/1‑143 (MQ=255)           gTCGGGATTCCACTTCTGTGGAATTAACGGTTATAGTGGCATAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgtcgt                                                                                                                               >  1:590518/1‑143 (MQ=255)           gTCGGGATTCCACTTCTGTGGAATTAACGGTTATAGTGGCATAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgtcgt                                                                                                                               >  1:446519/1‑143 (MQ=255)             cGGGATTCCACTTCTGTGGAATTAACGGTTATAGTGGCATAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGtcgtcg                                                                                                                             <  2:590518/143‑1 (MQ=255)             cGGGATTCCACTTCTGTGGAATTAACGGTTATAGTGGCATAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGtcgtcg                                                                                                                             <  2:446519/143‑1 (MQ=255)                      aCTTCTGTGGAATTAACGGTTATAGTGGCATAATCAGCCGCAGGAATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCAGGACCTGGGAATATGCGCTAATCGGGTTAGTCGTCGGcaccagca                                                                                                                    >  2:442481/1‑140 (MQ=255)                       cttcttTGGAATTAACGGTTAGAGTGCCATAATCAGCCGCAGGCATGAAGCCCGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCGCTAATTGGGTTAGTCGTCGGCATCattat                                                                                                                   <  1:42178/143‑1 (MQ=255)                         tCTGTGGAATTAACGGTTATAGTGGCATAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCGCTAATTGGGTTAGTCGTCGGCATCATTATTg                                                                                                                 <  1:646420/143‑1 (MQ=255)                             tgGAATTAACGGTTATAGTGGCATAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCGCTAATGGGGTCAGTCGTCGGCACCAGTATAGGAGc                                                                                                             >  1:465988/1‑143 (MQ=255)                                     aCGGTTATAGTGGCATAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCGGGGAATATGCGCTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCa                                                                                                     >  2:372122/1‑143 (MQ=255)                                         ttATAGTGGCATAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCGCTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCg                                                                                                 >  1:26542/1‑143 (MQ=255)                                             aGTGGCATAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCGCTAATTGGGTTAGTCGTGGGCATCATTATTGGGGCTGTGGGCATGCGattt                                                                                             >  2:209786/1‑143 (MQ=255)                                                 gCATAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAAATGTTCATGACCTGGGAATATGCGCTAATTGGGTTAGTCGTCGGCATCATTAGGGGTGCTGTGGCCATGCGGTTTGGTa                                                                                         >  1:592408/1‑143 (MQ=255)                                                    tAAACAGCCGAAGGCATGGAGCCAGCAGCCAACGCCCGACGTAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCGCTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATc                                                                                      <  1:436495/143‑1 (MQ=255)                                                                       aGCCTGAAGCGTACACCATACGGGAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCGCTAATTGGGTTAGCCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGCAACa                                                                   <  2:547018/143‑1 (MQ=255)                                                                            gAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCGCTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGCAACAGGCGt                                                              >  1:80144/1‑143 (MQ=255)                                                                                ccAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCGCTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGCAACAGGCGTTGCa                                                          <  1:355708/143‑1 (MQ=255)                                                                                  aaCACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCGCTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGCAACAGGCGTTGCAGt                                                        <  1:299535/143‑1 (MQ=255)                                                                                     acCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCGCTAATGGGGTTAGTCGTCGGCATAATTATTGGTGCTGTGGCCATGCGTTTTGGGAATCGTAAACGACGCCGGCACCAGGGGGTGCAGTagg                                                     >  2:561279/1‑141 (MQ=255)                                                                                              aaaCAAAAGACAAAGGGAGATGTTCATGACCTGGGAATATGAGCTAATGGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATACGTTTTGGGAATCGTAAACTACGCCAGCAACAGGCGTTGCAGTACGAGCTGGaaa                                            >  2:368054/1‑143 (MQ=255)                                                                                                  aaaaGACAACGGGAGATGTTCATGACCTGGGAATATGCGCTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGCAACAGGCGTTGCAGTACGAACTGGAAAAGaa                                        >  2:277540/1‑143 (MQ=255)                                                                                                  aaaaGACAACGGGAGATGTTCATGACCTGGGAATATGCGCTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGCAACAGGCGTTGCAGTAAGAACTGGACAAGaa                                        >  2:388394/1‑143 (MQ=255)                                                                                                   aaaGACAACGGGAGATGTTCATGACCTGGGAATATGCGCTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGCAACAGGCGTTGCAGTACGAACTGGAAAAGAAt                                       >  1:20892/1‑143 (MQ=255)                                                                                                         aaggggAGATGTTCATGACCTGGGAATATGCGCTAATTGGGTTAGTCGGCCGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGCAACAGGCGTTGCAGTACGAACTGGAAAAGAATAAAGCt                                 <  2:327226/140‑1 (MQ=255)                                                                                                                   gTTCATGACCTGGGAATATGCGCTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGCAACAGGCGTTGCAGTACGAACTGGAAAAGAATAAAGCTGAACTGGACg                       >  1:106119/1‑143 (MQ=255)                                                                                                                    ttCGTGACCGGGGAATATGCGCGAATGGGGTTAGACGGCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGCAACAGGCGTTGCAGTACGAACTGGAAAAGAATAAAGCTGAACTGGACGa                      <  1:255062/143‑1 (MQ=255)                                                                                                                            cTGGGAATATGCGCTAATTGGGTTAGTCGTCGGCATCATTATTGCTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGCAACAGGCGTTGCAGTACGAACTGGAAAAGAATAAAGCTAAAGGGGACGAGGATcgcg              >  2:517138/1‑143 (MQ=255)                                                                                                                                        gcTAATTGGGTTAGTCGTCGGCATCATTGTTGGTGCTCGAACCATGCGTTTTGGTAAGCGTAAACGACGCCAGCAACAGGCGTTGCAGTACGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTa  <  2:453847/143‑1 (MQ=255)                                                                                                                                        |                                                                                                                                               GCGATATGCAGTCGGGATTCCACTTCTGTGGAATTAACGGTTATAGTGGCATAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCGCTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGCAACAGGCGTTGCAGTACGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTA  >  NC_000913/3380073‑3380350

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