Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	5,267,269	A→C	30.9%	D197A (GAC→GCC)	rlmG →	Ribosomal RNA large subunit methyltransferase G

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	5,267,269	0	A	C	30.9%	70.2 / 3.7	49	D197A (GAC→GCC)	rlmG	Ribosomal RNA large subunit methyltransferase G
Reads supporting (aligned to +/- strand):  ref base A (13/19);  new base C (14/3);  total (27/22)
Fisher's exact test for biased strand distribution p-value = 6.91e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCGACCTGATGGAGAAGTACATCGGCCCCGTACAGGCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCACGCAACCTCGGCCGTGCGCGGGTGGCAGACCTGGGCTGCGG  >  KT2440/5267138‑5267401                                                                                                                                    |                                                                                                                                     cggcGACCTGATGGAGAAGTACAGCGGCCCCGTACTGGCTTCGATGGCACTTAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGAcgc                                                                                                                                   <  7:344668/135‑1 (MQ=255)           aTGGAGACGTACATCGGCCCCGTACAGGCTTCGCGGGCACTCAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCactt                                                                                                                         <  4:154412/135‑1 (MQ=255)             ggAGAAGTACATCGGCCCCGTACAGGCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCCCGTCCCCCAGCCCCAGCTGCTACCGCCGCGCCGCGCCGGcccctgc                                                                                                                       >  1:165331/1‑129 (MQ=255)                   gTACATCGGCCCCGTACAGGCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGt                                                                                                                 <  6:176945/135‑1 (MQ=255)                    tACATCGGCCCCGTACAGGCTTCGCTGGCACTGAAAACGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTc                                                                                                                <  6:123601/135‑1 (MQ=255)                     aCATCGGCCCCGTACAGGCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCCAGTCCCCCGACCCCACCAGCTACCGCCTCGCCGCGCCGGCACGTGACTTGGTCa                                                                                                               >  7:356887/1‑135 (MQ=255)                         cGGCCCCGTACAGGCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCCAGTCCCCCGACCCCCGCGGCTACCGCCACGCCGCGCCCGCACTTGCCGTGGTCAcccc                                                                                                           >  5:247081/1‑134 (MQ=255)                          ggCCCCGTACAGGCTTCGCTGTCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTGGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACg                                                                                                         <  4:158499/136‑1 (MQ=255)                           gcgccGTACAGGCTGCGCGGGGACTGAAAAAGGCCCTTTTGCTGACTGCATCGGTTGCCGAGCGGCCGGGCGCCAAGTCCCCCTCCCCCAGCTGCTACCGCCTCGCCGCGCCGGCACTTGACTTGGTCAACCACg                                                                                                         <  7:362155/132‑1 (MQ=255)                            ccccGTACAGGCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCGCGCCAAGTCGCCCTACCCCCGCTGCTACCGCCTCGCCGCGCCGGCCCTTGACTTGGTCACCCCCGc                                                                                                        >  5:368953/1‑135 (MQ=255)                                 tACAGGCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACgt                                                                                                  <  3:140525/136‑1 (MQ=255)                                 tACAGGCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCGTCCGCCTCGCCGCGCCCGCACTTGACTTGGTCAACCAAGCCAAcc                                                                                                   >  8:269666/1‑134 (MQ=255)                                  aCAGGCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACgt                                                                                                  <  7:99861/135‑1 (MQ=255)                                     ggCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGtt                                                                                               <  6:89907/135‑1 (MQ=255)                                         tCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGCCGCGCCGGCACGTGCCTTGGTCAACCACGCCAACGTGTTCTGc                                                                                           >  7:362656/1‑135 (MQ=255)                                         tCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCGACCCCAGCGGCTACCGCCGCGACGCGCCGGCACTTGACTGGGTGAACCACGCCAACGGGTTCTGcc                                                                                          >  3:210334/1‑136 (MQ=255)                                         tCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAACTCCCCCTCCCCCAGCTGCTAACCCCCCGCCGCGCCGGCCCATGACGAGGTCAACCACCCCAACCTGTTCTGc                                                                                           >  1:163722/1‑135 (MQ=255)                                          cGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGcc                                                                                          <  3:133218/135‑1 (MQ=255)                                           gCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCGCGCCAAGTCCCCCTACCCCCCCTGCGACCGCCTCGCCGCGCCCGCCCGTGCCGTGGTCACCCACCCCCACGGGTTATGCCg                                                                                         >  2:327874/1‑135 (MQ=255)                                            cTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGt                                                                                        <  6:7203/135‑1 (MQ=255)                                             tGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGCGCGGCCGCTCGCCAAGTCCCCCTAACCCCGCGGCTACCGCCTCGACGCGCCGGCCCATGCCGTGGTGAACCACGCCCACGGGTTCGGcccg                                                                                        >  1:119591/1‑132 (MQ=255)                                              ggCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCCGCTGCTACCGCCTCGCCGCGCCGGCCCGTGACTTGGTCAACCCCGCCAACGTGTTCTGCCGt                                                                                        >  2:296134/1‑133 (MQ=255)                                              ggCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGCCGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGt                                                                                        >  2:164603/1‑133 (MQ=255)                                                  cTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGc                                                                                  <  2:34794/135‑1 (MQ=255)                                                   tGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGCGCGGCCGCTCGCCAAGTCCCCCTACCCCCCCTGCTACCGCCGCGACGCGCCGGCCCATTACTTGGTCAAACACGCCAACCGGTTCTGCCGGGAAGGcc                                                                                 >  7:41580/1‑135 (MQ=255)                                                   tGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGCGCGGCCGCACGCCAAGTCCCCCCACCCCACCAGCTACCCCCACGCCCCGCCCGCACATGACTTGGTCAACCACCCCCACCGGTTCTGCCGTGAAGGcc                                                                                 >  7:47639/1‑135 (MQ=255)                                                   tGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGcc                                                                                 >  6:105418/1‑135 (MQ=255)                                                        aaGGCCCGTTTGCTGACAGCAACGGTTGCCGCGCGGCCGCTCGCCCCGTCCCCCTCCCCCCGCTGCTACCGCCGCGCCGCGCCCGCACTTGACATGGTCAACCCCGCCACCGTGGTCTGCCGGGCAGGCCggggt                                                                            >  1:271602/1‑133 (MQ=255)                                                        aaGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCCACGTGTTCTGCCGTGAAGGCCTGGat                                                                            >  2:283094/1‑135 (MQ=255)                                                              cGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCCCTTGCCTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCggg                                                                      >  6:264632/1‑135 (MQ=255)                                                               gTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGAc                                                                    <  4:330543/136‑1 (MQ=255)                                                               gTTTGCGGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGAc                                                                    <  4:336181/136‑1 (MQ=255)                                                                  tGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGGGTTCTGCCGTGGAGGCCTGGATATCGGGAcgc                                                                  >  2:268814/1‑135 (MQ=255)                                                                            aaCGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCCGCCCTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCAGGATATCGGGACGCGGGCGTtcct                                                        >  7:278201/1‑135 (MQ=255)                                                                             aCGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTtcctc                                                       <  7:56025/135‑1 (MQ=255)                                                                                      gAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCGCGCCGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCGGGCTATCGGGACGCGGGCGGTCCTCCCGCATCTAc                                             >  4:233445/1‑136 (MQ=255)                                                                                         cGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCa                                           <  5:105418/135‑1 (MQ=255)                                                                                            ccGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCACGCCGCGCCGGCCCTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCACGc                                        >  7:192630/1‑135 (MQ=255)                                                                                                         cccccTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCACGCAACCTCGGCCGTgc                          <  3:337854/136‑1 (MQ=255)                                                                                                               aCCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCACGCAACCTCGGCCGTgc                          <  5:262798/130‑1 (MQ=255)                                                                                                               aCCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCACGCAACCTCGGCCGTGCGCggg                     >  1:234519/1‑135 (MQ=255)                                                                                                               aCCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCACGCAACCACGGCCGTgc                          >  6:262798/1‑130 (MQ=255)                                                                                                                     gctgctACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCACGCAACCTCGGCCGGGCGCGGGTGGCAg               >  1:79071/1‑135 (MQ=255)                                                                                                                         ctACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGtt                                                                                               >  7:369789/1‑51 (MQ=255)                                                                                                                         ctACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGtt                                                                                               <  8:369789/51‑1 (MQ=255)                                                                                                                          tcccGCCTCTCCGCGCCGGCACTTGACTGGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCACGCAACCTCGGCCGTGCGCGGGTGGCAGACCTg          <  2:169230/133‑1 (MQ=255)                                                                                                                          tACCGCCTCGCCGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCACGCAACCTCGGCCGTGCGCGGGTGGCAGACCTg          <  1:296134/135‑1 (MQ=255)                                                                                                                          tACCGACTCGACGCGCCGGCACTTGCCGTGGTCCACCACGCCAACGTGTGCGGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCACGCAACCTCGGCCGTGCGCGGGTGGCAGACCTg          <  8:259750/135‑1 (MQ=255)                                                                                                                                  cGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCACGCAACCTCGGCCGTGCGCGGGTGGCAGCCCTGGGCTGCgg  >  1:236514/1‑135 (MQ=255)                                                                                                                                    |                                                                                                                                     CGGCGACCTGATGGAGAAGTACATCGGCCCCGTACAGGCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCACGCAACCTCGGCCGTGCGCGGGTGGCAGACCTGGGCTGCGG  >  KT2440/5267138‑5267401

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