Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	193,327	C→T	D152D (GAC→GAT)	frr →	ribosome‑recycling factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	193,327	0	C	T	96.9%	92.6 / NA	32	D152D (GAC→GAT)	frr	ribosome‑recycling factor
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (19/12);  total (20/12)
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.30e-01

GTTCAATGTCTCCGGCCGTTGAAAAAGCGATTATGGCGTCCGATCTTGGCCTGAACCCGAACTCTGCGGGTAGCGACATCCGTGTTCCGCTGCCGCCGCTGACGGAAGAACGTCGTAAAGATCTGACCAAAATCGTTCGTGGTGAAGCAGAACAAGCGCGTGTTGCAGTACGTAACGTGCGTCGTGACGCGAACGACAAAGTGAAAGCACTGTTGAAAGATAAAGAGATCAGCGAAGACGACGATCGCCGTTCTCAGGACGATGTACAGAAACTGACTGATGCTGCAATCAAGAAAATTGAAGCGGCGCTGGCAGACAAAGAAGCAGAACTGATGCAGTTCTGATTTCTTGAACGACAAAAACGCCGCTCAGTAGATCCTTGCGGATCGGCTGGCGGCGTTTTGCTTTTTATTCTGTCTCAACTCTGGATGTTTCATGAAGCAACTCACCATTCTGGGCTCGACCGGCTCGATTGGTTGCAGC  >  NC_000913/193086‑193568                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  gTTCAATGTCTCCGGCCGTTGAAAAAGCGATTATGGCGTCCGATCTTGGCCTGAACCCGAACTCTGCGGGTAGCGACATCCGTGTTCCGCTGCCGCCGCTGACGGAAGAACGTCGTAAAGATCTGACCAAAATCGTTCGTGGTGAAGCAGAACAAGCGCGTGTTGCAGTACGTAACGTGCGTCGTGACGCGAACGACAAAGTGAAAGCACTGTTGAAAGATAAAGAGATCAGCGAAGACGATGATCGCCGt                                                                                                                                                                                                                                          <  1:245940/251‑1 (MQ=255)              ggCCGTTGACCAAGCGATTATTGCATCAGATCTTAGCCTGAACCCGAACTCTGCGGGTAGCGACATCCGTGTTCCGCTGCCGCCGCTGACGGCAGAACGTCGTAAAGATCTGACCAACATCGTTCGTGGTGAAGCAGAACAAGCGCGTGTTGCCGTAAGTAACGTGCGGCGTGACGCGAACGCCAAAGGGTACGCACTGTTGAAAGATAAAGAGATCAGCGAAGACGATGATCGCCGTTCTCAGGACGATg                                                                                                                                                                                                                             >  2:255503/1‑251 (MQ=255)                         aaGCGATTATGGCGTCCGATCTTGGCCTGAACCCGAACTCTGCGGGTAGCGACATCCGTGTTCCGCTGACGCCGCTGACGGAAGAACGTCGTAAAGATCTGACCAAAATCGTTCGTGGTGAAGCAGAACAAGCGCGTGTTGCAGTACGTAACATGCGTCGTGACGCGAACGACAAAGTGAAAGCACTGTTGAAAGATAAAGAGATCAGAGAAGACGATGATCGCCGTTCTCATGACGATGTACAGAAactg                                                                                                                                                                                                                  >  2:199180/1‑251 (MQ=255)                                   ggCGTCCGATCTTGGCCTGAACCCGAACTCTGCGGGTAGCGACATCCGTGTTCCGCTGCCGCCGCTGACGGAAGAACGTCGTAAAGATCTGACCAAAATCGTTCGTGGTGAAGCAGAACAAGCGCGTGTTGCAGTACGTAACGTGCGTCGTGACGCGAACGACAAAGTGAAAGCACTGTTGAAAGATAAAGAGATCAGCGAAGACGATGATCGCCGTTCTCAGGACGATGTACAGAAACTGACTGAtgctg                                                                                                                                                                                                        >  1:98694/1‑251 (MQ=255)                                     cGTCCGATATTGGCCTGAACCCGAACTCTGCGGGTAGCGACATCCGTGTTCCGCTGCCGCCGCTGACGGAAGAACGTCGTAAAGATCTGACCAAAATCGTTCGTGGTGAAGCAGAACAAGCGCGTGTTGCAGTACGTAACGTGCGTCGTGACGCGAACGACAAAGTGAAAGCACTGTTGAAAGATAAAGAGATCAGCGAAGACGATGATCGACGTTCTCAGGACGATGTACAGAAACTGACTGAtgctgc                                                                                                                                                                                                       >  2:169042/1‑250 (MQ=255)                                         cGATCTTGGCCTGAACCCGAACTCTGCGGGTAGCGACATCCGTGTTCCGCTGCCGCCGCTGACGGAAGAACGTCGTCAAGATCTGACCACAATAGTTCGTGGTGAAGCACAACAAGCGCGTGTTGCAGTACGTAACGTGCGTCGTGACGCGAACGACAAAGTGAAAGCACTGTTGAAACATAAAGAGATCAGCGAAGACGATGATCGCCGTTCTCAGGACGATGTACAGAAACTGACTGATGCTGCAATc                                                                                                                                                                                                   >  1:258777/1‑250 (MQ=255)                                         cGATCTTGGCCTGAACCCGAACTCTGCGGGTAGCGACATCCGTGTTCCGCTGCCGCCGCTGACGGAAGAACGTCGTAAAGATCTGACCAAAATCGTTCGTGGTGAAGCAGAACAAGCGCGTGTTGCAGTACGTAACGTGCGTCGTGACGCGAACGACAAAGTGAAAGCACTGTTGAAAGATAAAGAGATCAGCGAAGACGATGATCGCCGTTCTCAGGACGATg                                                                                                                                                                                                                             >  1:82727/1‑224 (MQ=255)                                         cGATCTTGGCCGGAACCCGAACTCTGCGGGTAGCGACATCCGTGTTCCGCTGCCGCCGCTGACGGAAGAACGTCGTAAAGATCTGACCAAAATCGTTCGTGGTGAAGCAGAACAAGCGCGTGTTGCAGTACGTAACGTGCGTCGTGACGCGAACGACAAAGTGAAAGCACTGTTGAAAGATAAAGAGATCAGCGAAGACGATGATCGCCGTTCTCAGGACGATg                                                                                                                                                                                                                             <  2:82727/224‑1 (MQ=255)                                          gctCTTGGCCTGAACCCGCACTCTGCGGGTCGCGACATCCGTGTTCCGCTGCCGCCGCTGACGGAAGAACGTCGTAAAGATCTGACCAAAATCGTTCGTGGTGAAGCAGAACAAGCGCGTGTTGCAGTACGTAACGTGCGTCGTGACGCGAACGACAAAGTGAAAGCACTGTTGAAAGATAAAGAGATCAGCGAAGACGATGATCGCCGTTCTCAGGACGATGTACAGAAACTGACTGATGCTGCAATCa                                                                                                                                                                                                  >  2:26227/3‑250 (MQ=255)                                                    tGAACCCGAACTCTGCGGGTAGCGACATCCGTGTTCCGCTGCCGCCGCTGACGGAAGAACGTCGTAAAGATCTGACCAAAATCGTTCGTGGTGAAGCAGAACAAGCGCGTGTTGCAGTACGTAACGTGCGTCGTGACGCGAACGACAAAGTGAAAGCACTGTTGAAAGATAAAGAGATCAGCGAAGACGATGATCGCCGTTCTCAGGACGATGTACAGAAACTGACTGATGCTGCAATCAAGAAAATTGa                                                                                                                                                                                        >  1:197546/1‑250 (MQ=255)                                                       aCCCGAACTCTGCGGGTAGCGACATCCGTGTTCCGCTGCCGCCGCTGACGGAAGAACGTCGTAAAGATCTGCCCAAAATCGTTCGTGGTGAAGCAGAACACGCGCGTGTTGCAGTACGTAACGTGCGTCGTGCCGCGAACGACAAAGTGAAAGCACTGTTGAAAGCTAAAGAGATCAGCGAAGACGATGATCGCCGTTCTCCGGACGATGTACAGAAACTGACTGATGCTGCAATCAAGAAAATTGAAgcg                                                                                                                                                                                    >  2:24073/1‑251 (MQ=255)                                                                        aGCGACATCCGTGTTCCGCTGCCGCCGCTGACGGAAGAACGTCGTAAAGATCTGACCAAAATCGTTCGTGGTGAAGCAGAACAAGCGCGTGTTGCAGTACGTAACGTGCGTCGTGACGCGAACGACAAAGTGAAAGCACTGTTGAAAGATAAAGAGATCAGCGAAGACGATGATCGCCGTTCTCAGGACGATGTACAGAAACTGACTGATGCTGCAATCAAGAAAATTGAAGCGGCGCTGGCAGACaa                                                                                                                                                                      >  2:99917/1‑248 (MQ=255)                                                                                    gtTCCGCTGCCGCCGCTGACGGAAGAACGTCGTAAAGATCTGACCAAAATCGTTCGTGGTGAAGCAGAACAAGCGCGTGTTGCAGTACGTAACGTGCGTCGTGACGCGAACGACAAAGTGAAAGCACTGTTGAAAGATAAAGAGATCAGCGAAGACGATGATCGCCGTTCTCAGGACGATGTACAGAAACTGACTGATGCTGCAATCAAGAAAATTGAAGCGGCGCTGGCATACAAAGAAGCAGAACTGa                                                                                                                                                        <  1:75731/250‑1 (MQ=255)                                                                                            gccgccGCGGACGGAAGAACGTCGTAAAGATCTGATCAAAATCGTTCGTGGTGAAGCAGAACAAGCGCGTGTTGCAGTACGTAACGTGCGTCGTGATGCTAACGAAAAAGTGAAAGCACTGTTGAAAGAGAAAGAGATCAGTGAAGACGATGATCGCCGTCCTCAGGACGATGTACAGAAACTGACTGATGCTGCAATCAAGAAAATTGAAGCGGCGCTGGCAGACAAAGAAGCAGAACTGATGCAGtt                                                                                                                                                 <  1:131012/249‑1 (MQ=255)                                                                                              cgccgcTGACGGAAGAACGTCGTAAAGATCTGCCCAAAATCGTTCGTGGTGAAGCAGAACAAGCGCGTGTTGCAGTACGTAACGTGCGTCGTGACGCGAACGACAAAGTGAAAGCACTGTTGAAAGATAAAGAGATCAGCGAAGACGATGATCGCCGTTCTCAGGACGATGTACAGAAACTGACTGATGCTGCAATCAAGAAAATTGAAGCGGCGCTGGCAGACAAAGAAGCAGAACTGATGCAGTTCTGa                                                                                                                                             <  1:169042/251‑1 (MQ=255)                                                                                                           agaaCGTCGTAAAGATCTGACCAAAATCGTTCGTGGTGAAGCAGAACAAGCGCGTGTTGCAGTACGTAACGTGCGTCGTGACGCGAACGACAAAGTGAAAGCACTGTTGAAAGATAAAGAGATCAGCGAAGACGATGATCGCCGTTCTCAGGACGATGTACAGAAACTGACTGATGCTGCAATCAAGAAAATTGAAGCGGCGCTGGCAGACAAAGAAGCAGAACTGATGCAGTTCTGATTTCTTGAACGAc                                                                                                                                <  1:177237/251‑1 (MQ=255)                                                                                                                gtcgtAACGATCTGACCAAAATCGTTCGTGGTGAAGCAGAACAAGCGCGTGTTGCAGTACGTAACGTGCGTCGTGACGCGAACGACAAAGTGAAAGCACTGTTGAAAGATAAAGAGATCAGCGAAGACGATGATCGCCGTTCTCAGGACGATGTACAGAAACTGACTGATGCTGCAATCAAGAAAATTGAAGCGGCGCTGGCAGACAAAGAAGCAGAACTGATGCAGTTCTGATTTCTTGAACGACaaaa                                                                                                                            >  2:110738/1‑250 (MQ=255)                                                                                                                                                      gAACAAGCGCGTGTTGCAGTACGTAACGTGCGTCGTGACGCGAACGACAAAGTGAAAGCACTGTTGAAAGATAAAGAGATCAGCGAAGACGATGATCGCCGTTCTCAGGACGATGTACAGAAACTGACTGATGCTGCAATCAAGAAAATTGAAGCGGCGCTGGCAGACAAAGAAGCAGAACTGATGCAGTTCTGATTTCTTGAACGACAAAAACGCCGCTCAGTAGATCCTTGCGGATCGGCTGGCGGCGt                                                                                     <  2:14878/251‑1 (MQ=255)                                                                                                                                                        aCAAGCGCGTGTTGCAGTACGTAACGTGCGTCGTGACGCGAACGACAAAGTGAAAGCACTGTTGCAAGATAAAGAGATCAGCGAAGACGATGATCGCCGTTCTCAGGAAGCTGTACGGAAACTGACTGATGCTACAATCAAGAAAATTGAAGCGGCGCTGGCAGCCAACGGAGCAGAACTGCTGCAGTTCTGATTTCTTGAACGACAAAAACGCCGCTCAGTAGATCCTTGCGGATCGGCAGGCGGCGttt                                                                                   >  1:261007/1‑251 (MQ=255)                                                                                                                                                           aGCGCGTGTTGCAGTACGTAACGTGCGTCGTGACGCGAACGACAAAGTGAAAGCACTGTTGAAAGATAAAGAGATCAGCGAAGACGATGATCGCCGTTCTCAGGACGATGTACAGAAACTGACTGATGCTGCAATCAAGAAAATTGAAGCGGCGCTGGCAGACAAAGAAGCAGAACTGATGCAGTTCTGATTTCTTGAACGACAAAAACGCCGCTCAGTAGATCCTTGCGGATCGGCTGGCGGCGTTTTGc                                                                                <  2:20795/251‑1 (MQ=255)                                                                                                                                                             cgcgTGTTGCAGTACGTAACGTGCGTCGTGACGCGAACGACAAAGTGAAAGCACTGTTGAAAGATAAAGAGATCAGCGAAGACGATGATCGCCGTTCTCAGGACGATGTACAGAAACTGACTGATGCTGCAATCAAGAAAATTGAAGCGGCGCTGGCAGACAAAGAAGCAGAACTGATGCAGTTCTGATTTCTTGAACGACAAAAACGCCGCTCAGTAGATCCTTGCGGATCGGCTGGCGGCGTTTTGCtt                                                                              >  2:51005/1‑251 (MQ=255)                                                                                                                                                              gcgTGTTGCAGTACGTAACGTGCGTCGTGACGCGAACGACAAAGTGAAAGCACTGTTGAAAGATAAAGAGATCAGCGAAGACGATGATCGCCGTTCTCAGGACGATGTACAGAAACTGACTGATGCTGCAATCAAGAAAATTGAAGCGGCGCTGGCAGACAAAGAAGCAGAACTGATGCAGTTCTGATTTCTTGAACGACAAAAACGCCGCTCAGTAGATCCTTGCGGATCGGCTGGCGGCGTTTTGCttt                                                                             <  1:10277/251‑1 (MQ=255)                                                                                                                                                              gcgTGTTGCAGTACGTAACGTGCGTCGTGACGCGAACGACAAAGTGAAAGCACTGTTGAAAGATAAAGAGATCAGCGAAGACGATGATCGCCGTTCTCAGGACGATGTACAGAAACTGACTGATGCTGCAATCAAGAAAATTGAAGCGGCGCTGGCAGACAAAGAAGCAGAACTGATGCAGTTCTGATTTCTTGAACGACAAAAACGCCGCTCAGTAGATCCTTGCGGATCGGCTGGCGGCGTTTTGCttt                                                                             <  1:110738/251‑1 (MQ=255)                                                                                                                                                                           cGTAAAGTGCGTCGTGACGCGAACGACAAAGTGAAAGCACTGTTGAAAGATAAAGAGATCAGCGAAGACGATGATCGCCGTTCTCAGGACGATGTACAGAAACTGACTGATGCTGCAATCAAGAAAATTGAAGCGGCGCTGGCAGACAAAGAAGCAGAACTGATGCAGTTCTGATTTCTTGAACGACAAAAACGCCGCTCAGTAGATCCTTGCGGATCGGCTGGCGGCGTTTTGCTTTTTATTTTGtctc                                                                 >  2:246697/1‑250 (MQ=255)                                                                                                                                                                                   gCGTCGTGACGCGAACGACAACGTGAAAGAACTGTTGAAAGATACAGAGATCAACGAAGACGACGCTCGCCGCTCTCAGGACGATGTACAGACACTGACTGCTGCTGCAATCAAGAAAATTGAAGCGGCGCTGGCAGACAAAGAAGCAGAACTGATGCAGTTCTGATTTCTTGAACGACAAAAACGCCGCTCAGTAGATCCTTGCGGCTCGGATGGCGGCGTTTTGCTTTTTATTTTGTCTCAAATCTgg                                                         >  2:265375/1‑250 (MQ=255)                                                                                                                                                                                              cgAACGACAAAGTGAAAGCACTGGTGAAAGATAAAGAGATCAGCGAAGACGATGATCGCCGTTCTTAGGACGATGTACGGAAACTGACTGATGCTTCAATCAAGAAAATTGAAGCGGCGCTGGCAGACAAAGAAGCAGAACGGATGCAGTTCTTAGTTCTTGAACGATAAAAACGCCGCTCAGTAGATCCTTGCTGATCGGCTGGCGGCGTTTTGCTTTTTATGCCGTCTCAACTCTGGATGTTTCATGaa                                             <  2:261007/251‑1 (MQ=255)                                                                                                                                                                                                    aCAAAGTGAAAGCACTGTTGAAAGATAAAGAGATCAGCGAAGACGATGATCGCCGTTCTCAGGACGATGTACAGAAACTGACTGATGCTGCAATCAAGAAAATTGAAGCGGCGCTGGCAGACAAAGAAGCAGAACTGATGCAGTTCTGATTTCTTGAACGACAAAAACGCCGCTCAGTAGATCCTTGCGGATCGGCTGGCGGCGTTTTGCTTTTTATTCTGTCTCAACTCTGGATGTTTCATGAAGCAACt                                       <  1:25847/251‑1 (MQ=255)                                                                                                                                                                                                           gAAAGCACTGTTGAAAGATAAAGAGATCAGCGAAGACGATGATCGCCGTTCTCAGGACGATGTACAGAAACTGACTGATGCTGCAATCAAGAAAATTGAAGCGGCGCTGGCAGACAAAGAAGCAGAACTGATGCAGTTCTGATTTCTTGAACGACAAAAACGCCGCTCAGTAGATCCTTGCGGATCGGCTGGCGGCGTTTTGCTTTTTATTCTGTCTCAACTCTGGATGTTTCATGAAGCAACTCACCAtt                                >  1:44733/1‑251 (MQ=255)                                                                                                                                                                                                                      tGAAAGATAACGAGATCAGCGAAGACGATGATCGCCGTTCTCAGGACGATGTACAGAAACTGACTGATGCTGCAATCAAGAAAATTGAAGCGGCGCTGGCAGACAAAGAAGCAGAACTGATGCAGTTCTGATTTATTGAACGACAAAAACGCCGCTCAGTAGATCCTTGCGGATCGGCTGGCGGCGTTTTGCTTTTTTTTCTGTCTCAACTCTGGATGTTTCATGAAGCAACTCACCATTCTGGGCACGAc                     >  1:236400/1‑251 (MQ=255)                                                                                                                                                                                                                                 gagaTCAGAGAAGACGATGATCGCCGTTCTCAGGACAATGTACAGAAACTGACTGATGCTGCAATCACGAAAATTGAAGCGGCGCTGGCAGACAAAGAAGCAGAACTGATGCAGTTCTGATTTCTTGAACGACAAAAACGCCGCTCAGTAGATCCTTGCGGATCGGCTGGCGGCGTTTTGCTTTTTATTCTGTCTCAACTCTGGATGTTTCATGAAGCAACTCACCATTCTGGGCTCGACCGGCTCGATTg          >  2:27748/1‑251 (MQ=255)                                                                                                                                                                                                                                   gaTCAGCGAAGACGATGATCGCCGTTCTCAGGACGATGTACAGAAAATGACTGATGCTGCAATCACGAAACTTGAAGCGGCGCTGGCAGACAAAGAAGCAGAACTGATGCAGTTCTGATTTCTTGAACGACAAAAACGCCGCTCAGTAGATCCTTGCGGATCGGCTGGCGGCGTTTTGCTTTTTATTCTGTCTCAACTCTGGCTGTTTCATGACGCAACTCACCATTCTGGGCTCGACCGGCTCGATTGGt        >  1:242757/1‑251 (MQ=255)                                                                                                                                                                                                                                          gAAGACGATGATCGCCGTTCTCAGGACGATGTACAGAAACTGACTGATGCTGCAATCAAGAAAATTGAAGCGGCGCTGGCAGACAAAGAAGCAGAACTGATGCAGTTCTGATTTCTTGAACGACAAAAACGCCGCTCAGTAGATCCTTGCGGATCGGCTGGCGGCGTTTTGCTTTTTATTCTGTCTCAACTCTGGATGTTTCATGAAGCAACTCACCATTCTGGGCTCGACCGGCTCGATTGGTTgcagc  >  1:244971/1‑250 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  GTTCAATGTCTCCGGCCGTTGAAAAAGCGATTATGGCGTCCGATCTTGGCCTGAACCCGAACTCTGCGGGTAGCGACATCCGTGTTCCGCTGCCGCCGCTGACGGAAGAACGTCGTAAAGATCTGACCAAAATCGTTCGTGGTGAAGCAGAACAAGCGCGTGTTGCAGTACGTAACGTGCGTCGTGACGCGAACGACAAAGTGAAAGCACTGTTGAAAGATAAAGAGATCAGCGAAGACGACGATCGCCGTTCTCAGGACGATGTACAGAAACTGACTGATGCTGCAATCAAGAAAATTGAAGCGGCGCTGGCAGACAAAGAAGCAGAACTGATGCAGTTCTGATTTCTTGAACGACAAAAACGCCGCTCAGTAGATCCTTGCGGATCGGCTGGCGGCGTTTTGCTTTTTATTCTGTCTCAACTCTGGATGTTTCATGAAGCAACTCACCATTCTGGGCTCGACCGGCTCGATTGGTTGCAGC  >  NC_000913/193086‑193568

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