Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,180,426	T→A	intergenic (+67/+53)	cobB → / ← ycfZ	protein‑lysine deacetylase/desuccinylase/lipoamidase/putative inner membrane protein YcfZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,180,426	0	T	A	95.7%	60.6 / ‑6.7	23	intergenic (+67/+53)	cobB/ycfZ	protein‑lysine deacetylase/desuccinylase/lipoamidase/putative inner membrane protein YcfZ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (11/11);  minor base G (1/0);  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 = 1.00e+00

GGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAGGGATTAAAAGCGGGAAGCATTGCCTGATGCACTACGCTTATCAGGCCTGCAAGGGGACGCAATTAATTGAATTTGCGGGTTTTGTGGGCCGGATAAGGCGTTCGCGCCCCATCCGGAATAATACACGCGCACTTTTCCGGCCTCTTTTATCGTTTACTCGTCAATAACAATGCGGCAATAAACGTACAGACACCAATTAGCCCTGCCCAGAGTGGCGATATTCTGGTAATCGCAAAGCTATTATCGATATACAAACATAAAGCCTGATAACCCATATTGGCAATAGTCATCGCCACTGCAAATTTTAAC  >  NC_000913/1180216‑1180641                                                                                                                                                                                                                   |                                                                                                                                                                                                                        ggCGCGCACACCGTGGAACTGAATCTTGAACCGATTCATGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGGTGCTGAAGGGATTAAAAGCGGGAAGCATTGCCTGATGCACTACGCTTATCAGGCCTGCAAGGGGACGCAATTAATTGAATTTGCGGGTTTTGTGGGCCGGAAAAGGCGTTCGCGCCCCATCCGGAATAATACACGCGCACtn                                                                                                                                                                                 <  1:263578/251‑2 (MQ=255)        acacCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAGGGATTAAAAGCGGGAAGCATTGCCTGATGCACTACGCTTATCAGTCCTGCAAGGGGACGCAATTAATTTAATTTGCGGGTGTTGTGGGCCGGAAAAGGCGTGCGCGCCCCATCCGGAATAATACACGCGCACTTTTCCggt                                                                                                                                                                          <  1:67559/251‑2 (MQ=255)           ccGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTGGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAGGGAGTAAAAGCGGGAAGCATTGCCGGATTCACTACGCTTATCAGGCCTGCAAGGGGACGCAAGTAATTGAATTTGCGGGCTTTGCGGGCCAGAAAAGGCGTTCGCGCCCCATCCGGAATAATACACGCGCACTTTTCCGGCctc                                                                                                                                                                       <  1:241774/251‑1 (MQ=255)                         cTTGAACCGAGTCAGGTTGGTAATGAATGTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTTGTGAAAAGTTGCTGAAGGGATTAAAAGCTGGAAGCATTGCCTGATGCACTACGCTGATCAGGCCTGCAAGGGTACGCAATGAATTGAATGTGCGGGTTTTGTGGGCCGGAAAAGGCGTTCGCGCCCCTTCCGGAATAATACACGCGCACTTTTCCGGCCTCTTTTATcn                                                                                                                                                               <  1:48681/245‑2 (MQ=255)                                 gAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGGGGTGACTGGGTTTGTTGCAGAGTTGCTGAAGGGTTTAAAAACGGGAAGCATTGCCTGATGCACTGCGCTTAAAGGGCCTGCAAGGGGACGCAATTAATTTAATTTGCGGGTTTTGTGGGCCGGGGAAGGCGTTCGCGCCCCATCCGGAATAATACACGCGCACTTTTCCGGCCTCtt                                                                                                                                                                     >  1:333461/1‑231 (MQ=255)                                            gTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAGGGATTAAAAGCGGGAAGCATTGCCTGATGCACTACGCTTATCAGGCCTGCAAGGGGACGCAATTAATTGAATTTGCGGGTTTTGTGGGCCGGAAAAGGCGTTc                                                                                                                                                                                                                <  1:74921/177‑1 (MQ=255)                                                         cGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAGGGATTAAAAGCGGGAAGCATTGCCTGATGCACTACGCTTATCAGGCCTGCAAGGGGACGCAATTAATTGAATTTGCGGGTTTTGTGGGCCGGAAAAGGCGTTCGCGCCCCATCCGGAATAATACACGCGCACTTTTCCGGCCTCTTTTATCGTTTACTCGTCAATAACAATGCGGCAATAAACGTAc                                                                                                                            >  1:9879/1‑248 (MQ=255)                                                                       gCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGATGCTGAAGGGATTAAAAGAGGGAAGCATTGCCTGATGCACTACGCTTCTCAGGCCTGCAAGGGGACGCAATTAATTGACTTTGCGGGTTTTGTGGGCCGGAAAAGGCGTTCGCACCCCCTCCGGAATAATACACGCGCACTTTTCCGGCCTCTTTTATCGTTTACTCGTCAATAACAATGCGGCAATAAACGTACAGaca                                                                                                                       >  1:179742/1‑239 (MQ=255)                                                                        cccGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAACGTTGCTGAAGGGATTAAAAGCGGGAAGCATTGCCTGCTGCACTACGCTTATCAGGCCTGCCAGGGGACGCAATTAATTGAATTTGCGGGTTTTGTGGGCCGGAAAAGGCGTTCGCGCCCCATCCGGAATAATACACGCGCACTTTTCCGGCCTCTTTTATCGTTTACTCGTCAATAACAATGCGGCAATAAACGTACAGa                                                                                                                         >  1:372755/1‑236 (MQ=255)                                                                        cccGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAGGGATTAAAAGCGTGAAGCATTGCCTGATGCACTACGCTTATCAGGCCTGCAAGGGGACGCAATTAATTGAATTTGCGGGTTTTGTGGGCCGGAAAAGGCGTTCGCGCCCCATCCGGAATAATACACGCGCACTTTTCCGGCCTCTTTTATCGt                                                                                                                                                              <  1:249539/199‑1 (MQ=255)                                                                                  naGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAGGGATTAAAAGCGGGAAGCATTGCCTGATGCACTACGCTTATCAGGCCTGCAAGGGGACGCAATTAATTGAATTTGCGGGTTTTGTGGGCCGGAAAAGGCGTTCGCGCCCCATCCGGAATAATACAcg                                                                                                                                                                                        >  1:282639/2‑163 (MQ=255)                                                                                              naGTTTGTTGAAAAGTTGCTGAAGGGATTAAAAGCGGGAAGCATTGCCTGATGCACTACGCTTATCAGGCCTGCAAGGGGACGCAATTAATTGAATTTGCGGGTTTTGTGGGCCGGAAAAGGCGTTCGCGCCCCATCCGGAATAATACACGCGCACTTTTCCGGcc                                                                                                                                                                         >  1:110919/2‑166 (MQ=255)                                                                                                   nnttnNNNNGTTNCTGAAGGGATTAAAAGCGGGAAGCATTGCCTGATGCACTACGCTTATCAGGCCTGCAAGGGGACGCAATTAATTGAATTTGCGGGTTTTGTGGGCCGGAAAAGGCGTTCGCGCCCCATCCGGAATAATACAcg                                                                                                                                                                                        <  2:282639/144‑1 (MQ=255)                                                                                                                                    aaGCATTGCCTGATGCACTACGCTTATCAGGCCTGCAAGGGGACGCAATTAATTGAATTTGCGGGTTTTGTGGGCCGGAAAAGGCGTTCGCGCCCCATCCGGAATAATACACGCGCACTTTTCCGGCCTCTTTTATCGTTTACTCGTCAATAACAATGCGGCAATAAACGTACAGACACCAATTAGCCCTGCCCAGAGTGGCGATATTCTgg                                                                                     >  1:104789/1‑212 (MQ=255)                                                                                                                                           gCCTGATGCACTACGCTTATCAGGCATGTAAGGGGACGAAAGTAATTGAATTTGCGGGTTTTGTGGGCCGGAAAAGGCGTTCGCGCCCCATCCTGAATAATACTCGCGCACTTTTCCGGCCTCTTTTATCGTTTACTCGTCAATAACAATGCGGCAATAAACGTACATACACCAATTAGCCCTGCCCAGAGTGGCGATATACTGGTAATTGCAAAGCTATTATCGATATACAAACATAAAGCCTGATAAc                                        <  1:424266/250‑1 (MQ=255)                                                                                                                                                        cGNTTANCNNGNNTGCNAGGGGACGCAATTAATTGAATNNGCGGGTTTTNGGGGCCGGAAAAGGCGTTCGCGCCCCATCCGGAATAATACACGCGCACTTTTCNGGCCTCTTTNNTCNTTTACNCNNNANNNACNATGCGGCAATAAACGTACANACNNCAATTANCCCTGCCCAGAGTGGCGATATTCTgg                                                                                     <  2:104789/192‑1 (MQ=255)                                                                                                                                                                  nccTGCAAGGGGACGCAATTAATTGAATTTGCGGGTTTTGTGGGCCGGAAAAGGCGTTCGCGCCCCATCCGGAATAATACACGCGCACTTTTCCGGCCTCTTTTATCGTTTACTCGTCAATAACAATGCGGCAATAAACGTACAGACACCAATTAGCCCTGCCCAGAGTGGCGATATTCTGGTAATCg                                                                               >  1:264238/2‑188 (MQ=255)                                                                                                                                                                    cTGCAAGGNNACGCAATTANTTGAATTTGCGGTTTTTGTGGGCCGGAAAAGGCGTTCGCGCCCCATCCGGAATNATACACGCGNNCTNTTCCGNCNNNTNNNATNGTTTACTCGTCAATAACAANGCNNCAATAANCGTACAGACACCAATTAGCCCTGccc                                                                                                       <  2:102951/162‑1 (MQ=255)                                                                                                                                                                    cTGCAAGGGGACGCAATTAATTGAATTTGCGGGTTTTGTGGGCCGGAAAAGGCGTTCGCGCCCCATCCGGAATAATACACGCGCACTTTTCCGGCCTCTTTTATCGTTTACTCGTCAATAACAATGCGGCAATAAACGTACAGACACCAATTAGCCCTGccc                                                                                                       >  1:102951/1‑162 (MQ=255)                                                                                                                                                                             gACGCAATTAATTGAATTTGCGGGTTTTGTGGGCCGGAAAAGGCGTTCGCGCCCCATCCGGAATAATACACGCGCACTTTTCCGGCCTCTTTTATCGTTTACTCGTCCATAACAATGCgg                                                                                                                                        >  1:35094/1‑120 (MQ=255)                                                                                                                                                                                 cAATTAATTGAATTTGCGGGTTTTGTGGGCCGGAAAAGGCGTTcg                                                                                                                                                                                                               >  1:386333/1‑45 (MQ=255)                                                                                                                                                                                                    gTTTTGTGGGCCGGAAAAGGCGTTCGCGCCCCCTCCGGAATAATACACGCGCTCTTTTCCGGCCTCTTTTATCGTTTACTCGTCAATAACAATGCTGCAATAAACGTTCAAGCACCAATTAGACCTGCCCAGAGTGGCGATATTCCGGTAATCGCAATGCTATTAGCGATATACAAACAGTAAGCCTTATAACCCATATTGGCAATAGTCATCGCCACTGCAAATTTTAAc  <  1:391004/231‑1 (MQ=255)                                                                                                                                                                                                    gTTTTGTGGGCCGGAAAAGGCGTTCGCGCCCCATCCGGAATAATACACGCGCACTTTTCCGGCCTCTTTTATCGTTTACTCGTCAATAACAATGCGGCAATAAACGTACAGACACCAATTAGCCCTGCCCAGANTGGNNNNNtt                                                                                         >  2:279664/1‑144 (MQ=255)                                                                                                                                                                                                              ccNNNNNAGGNGTTCGCGCCCCATCCGGAATAATACACGCGCACTTTTCCGGCCTCTTTTATCGTTTACTCGTCAATAACAATGCGGCAATAAACGTACAGACACCAATTAGCCCTGCCCAGAGTGGCGATATTCTGGTAATcc                                                                               <  2:264238/144‑2 (MQ=255)                                                                                                                                                                                                                   |                                                                                                                                                                                                                        GGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAGGGATTAAAAGCGGGAAGCATTGCCTGATGCACTACGCTTATCAGGCCTGCAAGGGGACGCAATTAATTGAATTTGCGGGTTTTGTGGGCCGGATAAGGCGTTCGCGCCCCATCCGGAATAATACACGCGCACTTTTCCGGCCTCTTTTATCGTTTACTCGTCAATAACAATGCGGCAATAAACGTACAGACACCAATTAGCCCTGCCCAGAGTGGCGATATTCTGGTAATCGCAAAGCTATTATCGATATACAAACATAAAGCCTGATAACCCATATTGGCAATAGTCATCGCCACTGCAAATTTTAAC  >  NC_000913/1180216‑1180641

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