Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,282,332	(A)5→6	coding (744/1293 nt)	yjcF ←	pentapeptide repeat‑containing protein YjcF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,282,332	1	.	A	93.8%	51.8 / ‑2.5	16	coding (744/1293 nt)	yjcF	pentapeptide repeat‑containing protein YjcF
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (7/8);  total (8/8)
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.67e-01

AGTCGGTACATATAACCATTTCGAGGAACGTATTTGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AAAAA‑TTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTAAAATTCGTTTTATACATATTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCG  >  NC_000913/4282097‑4282576                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     agTCGGTACATATAACCATTTCGAGGAACGTATTTGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AAAAAATTTGAATAATCa                                                                                                                                                                                                                                          >  1:6976/1‑249 (MQ=255)      gTACATATAACCATTTCGAGGAACGTATTTGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATACAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAAATAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AAAAAATTTGAATAATCaatgaagg                                                                                                                                                                                                                                   <  2:479900/251‑1 (MQ=255)         cATATAACCATTTCGAGGAACGTATTTGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AAAAAATTTGAATAATCAAGGAAGGTAt                                                                                                                                                                                                                                >  1:238056/1‑251 (MQ=255)                         ggAACGTATTTGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGctct                                                                                                                                                                                                                <  2:309475/251‑1 (MQ=255)                          gAACGTATTTGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTa                                                                                                                                                                                                               <  1:199592/251‑1 (MQ=255)                                                                                                                      cAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATc                                                                                                                   <  2:6976/251‑1 (MQ=255)                                                                                                                       aGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGCTTTACCTCAGCCATGTATGCCATAT‑AAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCg                                                                                                                  >  1:62286/1‑251 (MQ=255)                                                                                                                                              cAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCAt                                                                                                                                                                                                  <  1:19589/148‑1 (MQ=255)                                                                                                                                              cAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCAt                                                                                                                                                                                                  >  2:19589/1‑148 (MQ=255)                                                                                                                                                   ttAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTAtttttt                                                                                      >  2:278588/1‑251 (MQ=255)                                                                                                                                                      gagaAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCa                                                                                   <  2:88965/251‑1 (MQ=255)                                                                                                                                                                       tctAAATAAATTAATATATGGCGCTATTACTTTATCAAGATTTACCTCAGCCATGTATGCCATATCAAAAA‑TTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTaaa                                                                   >  2:79673/1‑250 (MQ=255)                                                                                                                                                                                    atatatGGCGCTATTACTTTATAAAGATTTACCTCAGCCCTGTATGCCATAT‑AAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGGATTTTTTCCATTATGGCGAAATTAAAATTCGTTTTATAc                                                      >  1:459681/1‑250 (MQ=255)                                                                                                                                                                                            cgcTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTAAAATTCGTTTTATACATATTACAg                                             <  1:49631/251‑1 (MQ=255)                                                                                                                                                                                                                ttACCTCAGCCATGTATGCCATAT‑AAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTAGCAAAAAGTATTTTTTCCATTATGGCGAAATTAAAATTCGTTTTATACATATTACAGTCCTGAAAAATAGCGtttt                          <  1:278588/250‑1 (MQ=255)                                                                                                                                                                                                                                tGCCATAT‑AAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAACAGTATTTTTTCCATTATGGCGAAATTAAAATTCGTTTTATACATATTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGc          >  2:175467/1‑250 (MQ=255)                                                                                                                                                                                                                                        aaaaaa‑TTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTACAATTCGTTTTATACATATTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCg  >  2:449518/2‑250 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     AGTCGGTACATATAACCATTTCGAGGAACGTATTTGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AAAAA‑TTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTAAAATTCGTTTTATACATATTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCG  >  NC_000913/4282097‑4282576

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