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%	52.4 / ‑2.3	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 = 1.95e-01

AGAGTCGGTACATATAACCATTTCGAGGAACGTATTTGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AAAAA‑TTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTAAAATTCGTTTTATACATATTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAAC  >  NC_000913/4282095‑4282564                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                         agagTCGGTACATATAACCATTTCGAGGAACGTATTTGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AAAAAATTTGAATAATc                                                                                                                                                                                                                               >  2:346744/1‑250 (MQ=255)    gTCGGTACATATAACCATTTCGAGGAACGTATTTGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AAAAAATTTGAATAATCaagg                                                                                                                                                                                                                           <  1:264956/251‑1 (MQ=255)                        cGAGGAACGTATTTGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTgg                                                                                                                                                                                                        <  1:315169/250‑1 (MQ=255)                                                                    aTTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAAc                                                                                                                                                           <  1:58485/251‑1 (MQ=255)                                                                       aaGTTTATAGTCTTGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTATAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCAGAT‑AAAAAATTTGAATAATCAAGGAAGGGATTTTTAATATTGGCTCTAGTCATTGGGGCATTGTAGAATTTCACATTCATTGCTGAACATg                                                                                                                                                        <  2:287367/251‑1 (MQ=255)                                                                                               ttATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCAt                                                                                                                                  <  1:349468/249‑1 (MQ=255)                                                                                                                                                            aGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTAt                                                                   <  1:198258/251‑1 (MQ=255)                                                                                                                                                               tAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGc                                                                >  2:200528/1‑251 (MQ=255)                                                                                                                                                                 aGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATCAAAAA‑TTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCg                                                               >  1:395841/1‑250 (MQ=255)                                                                                                                                                                     cGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAAtt                                                          <  1:391864/251‑1 (MQ=255)                                                                                                                                                                           taaataaaTTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTAAAAtt                                                    >  1:154044/1‑251 (MQ=255)                                                                                                                                                                             aataaatTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTAAAATTCg                                                  >  1:436052/1‑251 (MQ=255)                                                                                                                                                                              ataaatTAATATATGGCTCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AAAAAATGTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATTCATAGCGGAACATGAATTTAGAGTTCCGGCAGGTATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATAAAAAAGTATGTTTTCCAGTATGGCGAAATGAAAATTCGt                                                 <  2:158682/251‑1 (MQ=255)                                                                                                                                                                                  atTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTAAAATTCGtttt                                              >  2:82422/1‑250 (MQ=255)                                                                                                                                                                                           tGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATTCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTAAAATTCGTTTTATACATAttt                                    >  2:322880/1‑250 (MQ=255)                                                                                                                                                                                                                             aTGTATGCCATAT‑AAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTAAAATTCGTTTTATACATATTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAAc  >  2:296199/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                         AGAGTCGGTACATATAACCATTTCGAGGAACGTATTTGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AAAAA‑TTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTAAAATTCGTTTTATACATATTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAAC  >  NC_000913/4282095‑4282564

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