Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,595,570	(A)6→5	coding (282/915 nt)	lgoR ←	putative DNA‑binding transcriptional regulator LgoR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,595,570	0	A	.	82.6%	54.5 / 12.8	23	coding (282/915 nt)	lgoR	putative DNA‑binding transcriptional regulator LgoR
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base . (10/9);  total (12/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

CTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATA‑TTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAAT  >  NC_000913/4595333‑4595818                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                         cTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATAGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATA‑TTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGCAAAA‑GCCTGCTCAAacac                                                                                                                                                                                                                                            >  1:54037/1‑251 (MQ=255)                   gAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATA‑TTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCAt                                                                                                                                                                                                                          >  1:216979/1‑251 (MQ=255)                    aGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATA‑TTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCAtt                                                                                                                                                                                                                         <  1:231501/251‑1 (MQ=255)                                    cAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATA‑TTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGGTCACTCATGGAAGAGGTAGTACAg                                                                                                                                                                                                          <  1:79998/250‑1 (MQ=255)                                      aaCAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATA‑TTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAGGCGGTAGTACAAGCg                                                                                                                                                                                                       >  1:249833/1‑251 (MQ=255)                                                             aTCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATA‑TTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGACAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGttt                                                                                                                                                                                >  1:282754/1‑251 (MQ=255)                                                                                   gggcggttcTTTTCGCGTTGAATAAGATTATACGGGCCGAAGTTTAAAA‑ATAATTCTCTTCCTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGGAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCCCTCATTGAAGCGGTAGTACAGGCGAATCCGTCATCATGTTCAGGTTTACGCGCAATAACATAGTCGttt                                                                                                                                                          <  2:216979/243‑1 (MQ=255)                                                                                        gcgcTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATA‑TTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACGGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAAc                                                                                                                                                     <  1:362571/251‑1 (MQ=255)                                                                                            ttttCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATA‑TTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGc                                                                                                                                                 >  1:109099/1‑251 (MQ=255)                                                                                              ttCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATA‑TTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCg                                                                                                                                                <  1:309215/250‑1 (MQ=255)                                                                                                 gCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATA‑TTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCaa                                                                                                                                            <  1:125183/251‑1 (MQ=255)                                                                                                                             tttAAAAGATA‑TTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTAAGGTTTACGCGCAATAACATAGTCGttt                                                                                                                                                          >  1:85733/1‑209 (MQ=255)                                                                                                                                   aGATA‑TTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATg                                                                                                           <  1:89405/251‑1 (MQ=255)                                                                                                                                                         gggCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCGTCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCAAACCTGCGTCAAGCCGCCGCATTCGCGTAAGTGGCTGTGAATATACCGCACTGTGGGGCGGCTgtt                                                                                      >  2:139635/1‑249 (MQ=255)                                                                                                                                                          ggCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGAtgt                                                                                    >  1:55941/1‑251 (MQ=255)                                                                                                                                                             tGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGAtgttgt                                                                                 >  1:374125/1‑251 (MQ=255)                                                                                                                                                                           cGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAgc                                                                   <  1:172965/251‑1 (MQ=255)                                                                                                                                                                                              aaaaCGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTg                                                >  1:157140/1‑251 (MQ=255)                                                                                                                                                                                                                  aaCTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGAcccc                              <  1:139635/250‑1 (MQ=255)                                                                                                                                                                                                                  aaCTGGCGCTGGTTTATCAGCGTGAAAAA‑GCCTGCTCAAACACTTTGGTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCACCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCcg                             <  1:337597/250‑1 (MQ=255)                                                                                                                                                                                                                         gcTGGTTTATCATCTTGAAAAA‑GCCTGCGCAATCACTGTGTTTTGTTCACTCATTGAATCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATATCATAGTCTTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTTGTTCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGTGAAGGGATATGACCCCTCGCCataa                     <  2:132757/251‑2 (MQ=255)                                                                                                                                                                                                                               ttATCATCGTGAAAAA‑GCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCCTAtcatca                >  1:283069/1‑250 (MQ=255)                                                                                                                                                                                                                                        tGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCTTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGCGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTgg       >  1:57604/1‑251 (MQ=255)                                                                                                                                                                                                                                             aaaGCCTGCTCAAACACTTTGTTTTGTCCGCGCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGGGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGTCCCCGCGCCATATCATCAATGACCTGGTTAAt  <  1:327528/251‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                         CTGAAGTGAATGTGTTTCGAGCATTTCGCGTAGCTCAAACAGCTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATA‑TTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAAT  >  NC_000913/4595333‑4595818

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