Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,325,054	G→A	W103* (TGG→TGA)	yciQ →	DUF2207 domain‑containing protein YciQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,325,054	0	G	A	96.5%	88.7 / NA	29	W103* (TGG→TGA)	yciQ	DUF2207 domain‑containing protein YciQ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (12/16);  minor base T (0/1);  total (12/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGGCATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGTTTGAAATTTATAAGACCCGTTTTTCTCTCCAGTTCTCTAATATTGCGGGTAATCCATTTAGCGAAATCGATCTT  >  NC_000913/1324815‑1325297                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    gAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTg                                                                                                                                                                                                                                          >  1:220837/1‑251 (MQ=255)          gagaTCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATTGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCgg                                                                                                                                                                                                                                 <  1:468413/251‑1 (MQ=255)                           agaagGGCTCTCGCTCTTTGATGTTGAAGCCACTCTGCAACAAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAACTCAGCAGCTTAAGAACGGCTTTTATCGTGATTTACCACGACTATGGCTGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCTTCTTGTTGGCGTCACCCGGGATGGTATTCCTGAACCCTGACCTCTTGACTGGCCTATCGGGTTAAGGAgaaatttcg                                                                                                                                                                                                                >  2:419960/1‑243 (MQ=255)                            gaagGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCAGTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGt                                                                                                                                                                                                               <  2:2771/251‑1 (MQ=255)                                      cGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCg                                                                                                                                                                                                          <  1:198717/246‑1 (MQ=255)                                      cGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATaaa                                                                                                                                                                                                     <  2:102190/251‑1 (MQ=255)                                           tttGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGGGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACg                                                                                                                                                                                                >  2:387916/1‑251 (MQ=255)                                                           ctctGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGTCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATTGATTCATCCTGATGGGGACGCTGCACTGCTGAACTATCATATTTTTGGCGGCACCCGTGATGTTATTCCTGAACCCTGACATCTTGACTGTCATAGCGGGTTAATGAGTATTGTCGTTTGCGATAAACAACGTTTCTTGCCTCAAg                                                                                                                                                                                  <  1:437951/249‑1 (MQ=255)                                                                       aTGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTAt                                                                                                                                                                    <  1:254137/251‑1 (MQ=255)                                                                                                          cAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCa                                                                                                                                                                        >  1:309129/1‑212 (MQ=255)                                                                                                          cAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCa                                                                                                                                                                        <  2:309129/212‑1 (MQ=255)                                                                                                                cgAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAggg                                                                                                                           >  1:81193/1‑251 (MQ=255)                                                                                                                   aaTCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAggggg                                                                                                                         >  1:412947/1‑250 (MQ=255)                                                                                                                    aTCAGCAGATTAAGGACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTAGTGGGCg                                                                                                                                                                                                          >  2:213157/1‑168 (MQ=255)                                                                                                                    aTCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCg                                                                                                                                                                                                          <  1:213157/168‑1 (MQ=255)                                                                                                                                                            ccTATGGATGCAGCCTGATGGAGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGc                                                                               <  1:248642/250‑1 (MQ=255)                                                                                                                                                               aTGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGGGAGGGGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGt                                                                            <  1:50251/251‑1 (MQ=255)                                                                                                                                                                      cAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCGGCTAATCTGGAACGTGACCGGTAACCACTGGCCGTTTGAAAt                                                                     <  1:419960/251‑1 (MQ=255)                                                                                                                                                                        gCCTGATGGGGACGCTTCACGGCGGAACTATCATATTGGTGGCGTCTCCCGTGATGGTGTTCCTGAATGCTGACATCTTGACTGGCATAACGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGGTTCTTGCCGCAAGGCGGCTATTATTATCAAATTCATTATAAGGTTAAAAATGCTTGCCTGCGTGAGGGGGATGCTGATCTGCTAATCTGTAACGTGCCCGGTAGACAGAGGCCGTTTGAAAttt                                                                   <  2:49153/251‑1 (MQ=255)                                                                                                                                                                             aTGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGTTTGAAATTtata                                                                >  1:413289/1‑249 (MQ=255)                                                                                                                                                                                 ggACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGTTTGAAATTTATAAGAccc                                                          <  1:277476/251‑1 (MQ=255)                                                                                                                                                                                  gACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGTCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGTTTGAAATTTATAAGACCCg                                                         >  1:404343/1‑251 (MQ=255)                                                                                                                                                                                   aCGCTGCACTGCTGAACTATCATATTGTTGGCGTCCCCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGGGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGTTTGAAATTTATAAGTCCCGt                                                        <  2:378577/251‑1 (MQ=255)                                                                                                                                                                                             gctgAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGTTTGAAATTTATAAGACCCGATTTTCTCTcc                                              >  1:341263/1‑251 (MQ=255)                                                                                                                                                                                                taaCTATCATGTTGTTGGCGTCACCCGGGATTGTATTCCTGAACCCTGTCAACTTGACTTGCATATCGGGTTTATGAGTATTTTCGTGGGCGATAAACAACGTTTCTTGCCTAAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGTGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGTTTGAAATTTATAAGACCCGTTTTTCTCTcct                                             <  2:361/248‑2 (MQ=255)                                                                                                                                                                                                     aTCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGTTTGAAATTTATAAGACCCGTTTTTCTCTCCAGTTctct                                      <  1:168375/251‑1 (MQ=255)                                                                                                                                                                                                                                gTATTCCTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGTTTGAAATTTATAAGACCCGTTTTTCTCTCCAGTTCTCTAATATTGCGGGTAATCCATTTAGCGa            >  2:311502/1‑250 (MQ=255)                                                                                                                                                                                                                                      cTGAACCCTGACATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGTTTGAAATTTATAAGACCCGTTTTTCTCTCCAGTTCTCTAATATTGCGGGTAATCCATTTAGCGAAATCGAt     >  1:226159/1‑251 (MQ=255)                                                                                                                                                                                                                                         aaCCCTGACTTCCTGACTGGCATATCGTGTTAATGAGTATTGTCTGGGGCTATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGCTTAAAAATGCTTTCCTGCGTTAGGTGGATTCTGATCTGCTAATCTGGAACGTGACCTGTATCCACTGGCCGTTTGAAATTTATAAGACCCGTTTTTCTCTCCAGTTCTCTAATATGGCGGGTAATCCATTTAGCGAAATCGATCtt  <  2:313253/251‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    GAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGGCATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGTGGGCGATAAACAACGTTTCTTGCCTCAAGGCGACTATCATTATCAAATTCATTATCAGGTTAAAAATGCTTTCCTGCGTGAGGGGGATTCTGATCTGCTAATCTGGAACGTGACCGGTAACCACTGGCCGTTTGAAATTTATAAGACCCGTTTTTCTCTCCAGTTCTCTAATATTGCGGGTAATCCATTTAGCGAAATCGATCTT  >  NC_000913/1324815‑1325297

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