Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,269,449	C→T	C12C (TGC→TGT)	aphA →	acid phosphatase/phosphotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,269,449	0	C	T	96.4%	92.9 / NA	28	C12C (TGC→TGT)	aphA	acid phosphatase/phosphotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (16/11);  minor base G (0/1);  total (16/12)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

TATCCTTACATCGAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACA‑CGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCGCCAGAAAGC  >  NC_000913/4269211‑4269689                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                  tATCCTTACATCGAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATCATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGCTGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGct                                                                                                                                                                                                                                       >  1:291593/1‑251 (MQ=255) tATCCTTACATCGAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGct                                                                                                                                                                                                                                       >  2:137932/1‑251 (MQ=255)                                                   gTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGGTAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCCCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAAccc                                                                                                                                                                                     >  1:125434/1‑251 (MQ=255)                                                                      taATCATCCTGCTTTTCCTCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAgg                                                                                                                                                                  >  1:365629/1‑251 (MQ=255)                                                                                   tttCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACa                                                                                                                                                       <  1:73202/249‑1 (MQ=255)                                                                                       atcaCAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCAcc                                                                                                                                                 >  1:73549/1‑251 (MQ=255)                                                                                            aacaaCCGCAGATAATACTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCACCAAATCTCACAATAAACAGTTTCAACCTACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCATCATCTCCTTCACCGCTTAACACTGGGACTAACGTTGCCCGGCTTACTGAAACGGCACCccttc                                                                                                                                            >  2:58037/4‑251 (MQ=255)                                                                                                 ccGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTgg                                                                                                                                       >  1:257838/1‑251 (MQ=255)                                                                                                         aaTCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATTTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGt                                                                                                                               >  1:329650/1‑251 (MQ=255)                                                                                                            ccttccttTCCCCGGCAGCTGGCTTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGc                                                                                                                            <  2:73549/251‑1 (MQ=255)                                                                                                                    tCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATTTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAAtt                                                                                                                     <  2:329650/250‑1 (MQ=255)                                                                                                                            aGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTCGGGAAAAATATGCGCAAGATCAAACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCACTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACCAATTGAAACTAGc                                                                                                            >  1:217794/1‑251 (MQ=255)                                                                                                                                    ttATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAgg                                                                                                    >  2:249675/1‑251 (MQ=255)                                                                                                                                                ggTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATGGAAAATAGCCTCGCAGGGGGTCCGCCAAt                                                                                        <  1:359665/251‑1 (MQ=255)                                                                                                                                                 gTTTTTGCAACAAATCCCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCGGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATg                                                                                       <  2:9292/251‑1 (MQ=255)                                                                                                                                                   ttttGCAACAAAGCGCACAATAAAAAGTTTCAACATACTGCCTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAAACAGTGCCGCTTGGTTATTGTCCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGtt                                                                                                                                    <  1:453868/204‑1 (MQ=255)                                                                                                                                                   ttatgaaacaaATCTCACACTAAACAGATTCAACATACTGACTATTTCGGGAAAAATATGCGCCAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCCTCTCCTTCACCGCTTACCCCTGGGACTACCGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGtt                                                                                                                                    >  2:453868/7‑204 (MQ=255)                                                                                                                                                    tttGCAACAACTATCACAGTTAAAAGTTTCTACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCAAGGCTTGCTGAACAGACACCCATTCATTCGGTTTCGGTCGCACAAATTGCAAATAGCATCGCAGTGCGTCCGCCAATGGCg                                                                                    >  2:372682/1‑251 (MQ=255)                                                                                                                                                                     aaTAAAACGTTAAAACATACTGACTATTTAGGAAAAAATATGCGCAAGATCACACAGGCATTCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGATCCGCTGTTGCCATGGCCTCATCTCCTTCACCGCTTAACACTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCACTGGCGGCGGGGTTTtataacg                                                                    >  2:30945/1‑247 (MQ=255)                                                                                                                                                                      agaaaaaGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGAt                                                                  <  2:16894/249‑1 (MQ=255)                                                                                                                                                                                 cAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACA‑CGGTACt                                                      <  1:221834/251‑1 (MQ=255)                                                                                                                                                                                             tATTTAGGGAAAAATATCCGCAAGCTAAAACAGGCACTCAGTGCCGTTTGTTTAATGTTCGCTCTAAAAAGTTCCGCTGTTGCCCTTGCCTCATCTACTTAACCGCTTAACCCTGGGACTAACGTTGCCAGGCGAGCTGCACAGGCACACATTCATTGGGTTTCGGTAGCAAAAATTGAAAATACCCTCGCAGGGCGTCCGACAATGGCGGCGGGGGTTGGTGTCGA‑GTCAACGGGACTTTTTTCCactac                                          >  2:25284/1‑247 (MQ=255)                                                                                                                                                                                              aTTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAGTAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACA‑CGGTACTTTTTTCCAGGCCg                                         >  1:232980/1‑251 (MQ=255)                                                                                                                                                                                                  aGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACA‑CGGTACTTTTTTCCAGTCCGGGc                                      <  2:365629/250‑1 (MQ=255)                                                                                                                                                                                                                       cacacaGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACA‑CGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAcc                 <  2:257838/250‑1 (MQ=255)                                                                                                                                                                                                                                aaTCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGACAATAGCGTCGCAGGGCGTCCGCCAATGGCGGGGGGGGTTTAGATCGGTGAAA‑CGGTACTTTTTTCCAGTCCAGGGTTATGGCGCGGcccaaaaaacatctggccgg       >  2:434099/1‑232 (MQ=255)                                                                                                                                                                                                                                aaTCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACA‑CGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCGCCAg       <  2:349093/251‑1 (MQ=255)                                                                                                                                                                                                                                     gTGCCGTTTTTTTATTGTTCGCCCTAAACAGTGACGCTGGTGCCCTGTCTTCAGCCCCTTCTCCGCTTATCCCTGGGACTAACGTTGCCAGGCTTGCTGAACATGCACCCATTCATTGGGTTGCGGTCGCACAAATTGAAAATTGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGGTTGATATCGATGACA‑CGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCGTCAGAAAGc  <  2:268142/251‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                  TATCCTTACATCGAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACA‑CGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCGCCAGAAAGC  >  NC_000913/4269211‑4269689

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