Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	551,583	A→G	11.3%	V188V (GTA→GTG)	rplC →	50S ribosomal protein L3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	551,583	0	A	G	11.3%	110.3 / 5.5	44	V188V (GTA→GTG)	rplC	50S ribosomal protein L3
Reads supporting (aligned to +/- strand):  ref base A (18/21);  new base G (3/2);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCAGGACAACACCCACGGTAACTCCGTGTCGCACCGTGTTCCTGGTTCCATCGGCCAGTGCCAGACTCCTGGTCGTGTGTTCAAGGGCAAGAAAATGTCCGGTCACATGGGCGCCGAGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTGGCGGCGACGTGGTTGTACGTCCAGCTGTCAAGGCTCGCGGTTAAGGGGAAACTGACATGCAACTTAATGTAAATGACG  >  NC_002947/551438‑551722                                                                                                                                                  |                                                                                                                                            gTCAGGACAACACCCACGGTAACTCCGTGTCGCACCGTGTTCCTGGTTCCATCGGCCAGTGCCAGACTCCTGGTCGTGTGTTCAAGGGCAAGAAAATGTCCGGTCACATGGGCGCCGAGCGCGTGACTGTTCAGTCCCTGGAAGTAGtt                                                                                                                                          >  8:53114/1‑149 (MQ=255)       acaacaCCCACGGTAACTCCGTGTCGCACCGTGTTCCTGGTTCCATCGGCCAGTGCCAGACTCCTGGTCGTGTGTTCAAGGGCAAGAAAATGTCCGGTCACATGGGCGCCGAGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTa                                                                                                                                    <  4:11989/149‑1 (MQ=255)                 cGGTAACTCCGTGTCGCACCGTGTTCCTGGTTCCATCGGCCAGTGCCAGACTCCTGGTCGTGTGTTCTAGGGCAAGAAAATGTCCGGTCACATGGGCGCCGAGCGCGTGACTGTTCAGTCCCTTGAAGTAGTTCGCGTAGACGCTGAAc                                                                                                                          <  6:137517/149‑1 (MQ=255)                        tCCGTGTCGCACCGTGTTCCTGGTTCCATCGGCCAGTGCCAGACTCCTGGTCGTGTGTTCAAGGGCAAGAAAATGTCCGGTCACATGGGCGCCGAGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACct                                                                                                                   >  7:287040/1‑149 (MQ=255)                          cGTGTCGCACCGTGTTCCTGGTTCCATCGGCCAGTGCCAGACTCCTGGTCGTGTGTTCAAGGGCAAGAAAATGTCCGGTCACATGGGCGCCGAGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACctgc                                                                                                                 <  3:313833/149‑1 (MQ=255)                                        ttCCTGGTTCCATCGGCCAGTGCCAGACTCCTGGTCGTGTGTTCAAGGGCAAGAAAATGTCCGGTCACATGGGCGCCGAGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGt                                                                                                   >  4:297304/1‑149 (MQ=255)                                         tCCTGGTTCCATCGGCCAGTGCCAGACTCCTTGTCGTGTGTTCAAGGGCAAGAAAATGTCCGGTCACATGGGCGCCGAGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGt                                                                                                   <  8:278738/148‑1 (MQ=255)                                           cTGGTTCCATCGGCCAGTGCCAGACTCCTGGTCGTGTGTTCAAGGGCAAGAAAATGTCCGGTCACATGGGCGCCGAGCGCGTGACTGTTCAGTCCCTGGAAGTGGTTcgcg                                                                                                                                      <  7:6523/111‑1 (MQ=255)                                           cTGGTTCCATCGGCCAGTGCCAGACTCCTGGTCGTGTGTTCAAGGGCAAGAAAATGTCCGGTCACATGGGCGCCGAGCGCGTGACTGTTCAGTCCCTGGAAGTGGTTcgcg                                                                                                                                      >  8:6523/1‑111 (MQ=255)                                                 ccATCGGCCAGTGCCAGACTCCTGGTCGTGTGTTCAAGGGCAAGAAAATGTCCGGTCACATGGGCGCCGAGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTg                                                                                                  >  3:9560/1‑141 (MQ=255)                                                 ccATCGGCCAGTGCCAGACTCCTGGTCGTGTGTTCAAGGGCAAGAAAATGTCCGGTCACATGGGCGCCGAGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTg                                                                                                  <  4:9560/141‑1 (MQ=255)                                                 ccATCGGCCAGTGCCAGACTCCTGGTCGTGTGTTCAAGGGCAAGAAAATGTCCGGTCACATGGGCGCCGAGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCt                                                                                          >  7:148368/1‑149 (MQ=255)                                                   aTCGGCCAGTGCCAGACTCCTGGTCGTGTGTTCAAGGGCAAGAAAATGTCCGGTCACATGGGCGCCGAGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTgg                                                                                        <  8:287040/149‑1 (MQ=255)                                                   aTCGGCCAGTGCCAGACTCCTGGTCGTGTGTTCAAGGGCAAGAAAATGTCCGGTCACATGGGCGCCGAGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTgg                                                                                        <  4:52222/149‑1 (MQ=255)                                                            tGCCAGACTCCTGGTCGTGTGTTCAAGGGCAAGAAAATGTCCGGTCACATGGGCGCCGAGCGCGTGACTGTTCAGTCCCTGGAAGTGGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTgg                                                                               >  7:9809/1‑149 (MQ=255)                                                              ccAGACTCCTGGTCGTGTGTTCAAGGGCAAGAAAATGTCCGGTCACATGGGCGCCGAGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTggc                                                                              <  1:233497/148‑1 (MQ=255)                                                                  aCTCCTGGTCGTGTGTTCAAGGGCAAGAAAATGTCCGGTCACATGGGCGCCGAGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTGGCGGCGa                                                                         <  7:53114/149‑1 (MQ=255)                                                                       tGGTCGTGTGTTCAAGGGCAAGAAAATGTCCGGTCACATGGGCGCCGAGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTGGCGGCGACGTgg                                                                    >  1:175559/1‑149 (MQ=255)                                                                                 ttCAAGGGCAAGAAAATGTCCGGTCACATGGGCGCCGAGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTGGCGGCGACGTGGTTGTACGTc                                                           <  3:297304/148‑1 (MQ=255)                                                                                          aaGAAAATGTCCGGTCACATGGGCGCCGAGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTGGCGGCGACGTGGTTGTACGTCCAGCTGTCaa                                                 <  5:81064/149‑1 (MQ=255)                                                                                                 tGTCCGGTCACATGGGCGCCGAGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTGGCGGCGACGTGGTTGTACGTCCAGCTGTCAAGGCTcgc                                          <  8:148368/149‑1 (MQ=255)                                                                                                   tCCGGTCACATGGGCGCCGAGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTGGCGGCGACGTGGTTGTACGTCCAGCTGTCAAGGCTCGCgg                                        >  7:257173/1‑149 (MQ=255)                                                                                                     cgggctcaTGGGCGCCGAGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTGGCGGCGACGTGGTTGTACGTCCAGCTGTCAAGGCTCGCGGtt                                      <  7:302990/143‑1 (MQ=255)                                                                                                     cGGTCACATGGGCGCCGAGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTGGCGGCGACGTGGTTGTACGTCCAGCTGTCAAGGCTCGCGGtt                                      <  2:175559/149‑1 (MQ=255)                                                                                                          acaTGGGCGCCGAGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTGGCGGCGACGTGGTTGTACGTCCAGCTGTCAAGTCTCGCGGTTAAggg                                 <  1:120917/149‑1 (MQ=255)                                                                                                            aTGGGCGCCGAGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTGGCGGCGACGTGGTTGTACGTCCAGCTGt                                                    >  5:276789/1‑128 (MQ=255)                                                                                                            aTGGGCGCCGAGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTGGCGGCGACGTGGTTGTACGTCCAGCTGTCAAGGCTCGCGGTTAAgggg                                >  5:292117/1‑148 (MQ=255)                                                                                                            aTGCGCGCCGAGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTGGCGGCGACGTGGTTGTACGTCCAGCTGt                                                    <  6:276789/128‑1 (MQ=255)                                                                                                              gggCGCCGAGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTGGCGGCGACGTGGTTGTACGTCCAGCTGTCAAGGCTCGCGGTTAAGGGGaaa                             >  6:328962/1‑149 (MQ=255)                                                                                                                gcgcCGAGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTGGCGGCGACGTGGTTGTACGTCCAGCTGTCAAGGCTCGCGGTTAAGGGGAAACt                           >  7:127208/1‑149 (MQ=255)                                                                                                                gcgcCGAGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTGGCGGCGACGTGGTTGTACGTCCAGCTGTCAAGGCTCGCGGTTAAGGGGAAACt                           <  6:292117/149‑1 (MQ=255)                                                                                                                  gcCGAGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTGGCGGCGACGTGGTTGTACGTCCAGCTGTCAAGGCTCGCGGTTAAGGGGAAACTGa                         >  5:107908/1‑149 (MQ=255)                                                                                                                      aGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTGGCGGCGACGTGGTTGTACGTCCAGCTGTCAAGGCTCGCGGTTAAGGGGAAACTGACATg                     <  5:94734/149‑1 (MQ=255)                                                                                                                      aGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTGGCGGCGACGTGGTTGTACGTCCAGCTGTCAAGGCTCGCGGTTAAGGGGAAACTGACATg                     <  4:2962/149‑1 (MQ=255)                                                                                                                         gcgTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTGGCGGCGACGTGGTTGTACGTCCAGCTGTCAAGGCTCGCGGTTAAGGGGAAACTGACATGCaa                  >  8:78265/1‑149 (MQ=255)                                                                                                                          cgTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTGGCGGCGACGTGGTTGTACGTCCAGCTGTCAAGGCTCGCGGTTAAGGGGAAACTGACATGCAAc                 <  8:257173/149‑1 (MQ=255)                                                                                                                              aCTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTGGCGGCGACGTGGTTGTACGTCCAGCTGTCAAGGCTCGCGGTTAAGGGGAAACTGACATGCAACTTaa             >  5:224676/1‑149 (MQ=255)                                                                                                                                    cAGTCCCTGGAAGTGGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTGGCGGCGACGTGGTTGTACGTCCAGCTGTCAAGGCTCGCGGTTAAGGGGAAACTGACATGCAACTTAATGTaa        <  8:9809/148‑1 (MQ=255)                                                                                                                                     aGTCCCTGGAAGTGGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTGGCGGCGACGTGGTTGTACGTCCAGCTGTCAAGGCTCGCGGTTAAGGGGAAACTGACATGCAACTTAATGTAAAt      >  7:132201/1‑149 (MQ=255)                                                                                                                                     aGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTGGCGGCGACGTGGTTGTACGTCCAGCTGTCAAGGCTCGCGGTTAAGGGGAAACTGACATGCAACTTAATGTAAAt      >  4:86748/1‑149 (MQ=255)                                                                                                                                         ccTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTGGCGGCGACGTGGTTGTACGTCCAGCTGTCAAGGCTCGCGGTTAAg                                   >  7:157199/1‑116 (MQ=255)                                                                                                                                         ccTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTGGCGGCGACGTGGTTGTACGTCCAGCTGTCAAGGCTCGCGGTTAAg                                   <  8:157199/116‑1 (MQ=255)                                                                                                                                         ccTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTGGCGGCGACGTGGTTGTACGTCCAGCTGTCAAGGCTCGCGGTTAAGGGGAAACTGACATGCAACTTAATGTAAATGACg  >  5:261998/1‑149 (MQ=255)                                                                                                                                         ccTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTGGCGGCGACGTGGTTGTACGTCCAGCTGTCAAGGCTCGCGGTTAAGGGGAAACTGACATGCAACTTAATGTAAATGACg  >  4:307752/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                            GTCAGGACAACACCCACGGTAACTCCGTGTCGCACCGTGTTCCTGGTTCCATCGGCCAGTGCCAGACTCCTGGTCGTGTGTTCAAGGGCAAGAAAATGTCCGGTCACATGGGCGCCGAGCGCGTGACTGTTCAGTCCCTGGAAGTAGTTCGCGTAGACGCTGAACGCAACCTGCTGCTGATCAAGGGTGCCGTACCTGGCGCTACTGGCGGCGACGTGGTTGTACGTCCAGCTGTCAAGGCTCGCGGTTAAGGGGAAACTGACATGCAACTTAATGTAAATGACG  >  NC_002947/551438‑551722

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