Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,351,197	C→T	intergenic (‑158/+210)	fabI ← / ← ycjD	enoyl‑[acyl‑carrier‑protein] reductase/DUF559 domain‑containing protein YcjD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,351,197	0	C	T	97.1%	95.6 / NA	34	intergenic (‑158/+210)	fabI/ycjD	enoyl‑[acyl‑carrier‑protein] reductase/DUF559 domain‑containing protein YcjD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (17/16);  minor base G (1/0);  total (18/16)
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

GTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTC  >  NC_000913/1351069‑1351325                                                                                                                                 |                                                                                                                                 gtgCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGTGCTTTTCCCCt                                                                                                                       <  1:237431/140‑1 (MQ=255)                    ttCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCa                                                                                                    >  1:862229/1‑140 (MQ=255)                             ttAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGa                                                                                           <  1:596177/140‑1 (MQ=255)                                        tACGAACAGATAAACGGTTATTATAATTAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTAt                                                                                >  2:509262/1‑140 (MQ=255)                                        tACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTAt                                                                                >  1:828550/1‑140 (MQ=255)                                        tACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTAt                                                                                >  2:526548/1‑140 (MQ=255)                                         aCGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATtg                                                                              >  2:884221/1‑140 (MQ=255)                                          cGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGa                                                                                           <  2:872727/127‑1 (MQ=255)                                          cGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGa                                                                                           >  1:872727/1‑127 (MQ=255)                                                                   cAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGcaaca                                                     >  2:353311/1‑140 (MQ=255)                                                                      ccTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCaacaaac                                                  >  1:425054/1‑140 (MQ=255)                                                                         ggCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAgg                                               <  1:641511/140‑1 (MQ=255)                                                                         ggCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAgg                                               <  1:65222/140‑1 (MQ=255)                                                                         ggCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAgg                                               <  2:828550/140‑1 (MQ=255)                                                                         ggCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGc                                              <  1:884221/140‑1 (MQ=255)                                                                          gCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGc                                              <  2:425054/140‑1 (MQ=255)                                                                           cTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCt                                             >  1:396683/1‑140 (MQ=255)                                                                                aGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATAc                                        >  1:576733/1‑140 (MQ=255)                                                                                aGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATAc                                        <  1:526548/140‑1 (MQ=255)                                                                                    gCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGcc                                    >  2:60764/1‑140 (MQ=255)                                                                                     cTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATAc                                        <  2:923528/135‑1 (MQ=255)                                                                                     cTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATAc                                        <  2:547171/135‑1 (MQ=255)                                                                                     cTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATAc                                        >  1:923528/1‑135 (MQ=255)                                                                                     cTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATAc                                        >  1:547171/1‑135 (MQ=255)                                                                                         aGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTaa                               <  2:862229/140‑1 (MQ=255)                                                                                              ccAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATAc                                        <  2:166730/126‑1 (MQ=255)                                                                                              ccAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATAc                                        >  1:166730/1‑126 (MQ=255)                                                                                              ccAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGAAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCtt                          >  2:707093/1‑140 (MQ=255)                                                                                                      gACCGGAGCAGGCTGCGGCAGGGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCg                  <  2:396683/140‑1 (MQ=255)                                                                                                            aGCAGGCTGCGGCAGGGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCt            >  2:442841/1‑140 (MQ=255)                                                                                                            aGCAGGCTGCGGCAGGGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCt            >  2:571031/1‑140 (MQ=255)                                                                                                                ggCTGCGGCAGGGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGa        <  2:210875/140‑1 (MQ=255)                                                                                                                ggCTGCGGCAGGGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGa        <  1:509262/140‑1 (MQ=255)                                                                                                                      ggCAGGGGGGGTGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTc  >  1:89222/1‑140 (MQ=255)                                                                                                                                 |                                                                                                                                 GTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTC  >  NC_000913/1351069‑1351325

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