Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,449,100	T→C	T24A (ACG→GCG)	rplP ←	50S ribosomal subunit protein L16

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,449,100	0	T	C	100.0%	137.3 / NA	38	T24A (ACG→GCG)	rplP	50S ribosomal subunit protein L16
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/14);  total (24/14)

GTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTGGTTGTAACATCAGCGACGCTCCTTATTTACGGCCTTTACGCTGCTGCTTTTTAGGCTGAGCAGCCGGTTTTTCCGGTTGTTCAAC  >  NC_000913/3448956‑3449244                                                                                                                                                 |                                                                                                                                                 gTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGCACCCTg                                                                                                                                            <  1:127745/151‑1 (MQ=255)           gAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACcg                                                                                                                                  <  2:436328/151‑1 (MQ=255)                        cAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCtt                                                                                                                     >  1:247346/1‑151 (MQ=255)                            tCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTtgtg                                                                                                                 >  2:359550/1‑151 (MQ=255)                            tCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTtgtg                                                                                                                 >  1:649620/1‑151 (MQ=255)                               tACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTtgt                                                                                                                  <  2:296683/147‑1 (MQ=255)                               tACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCa                                                                                                               >  1:371016/1‑150 (MQ=255)                                      aCGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACgg                                                                                                       <  2:85496/151‑1 (MQ=255)                                        gCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGa                                                                                                      >  2:144310/1‑150 (MQ=255)                                        gCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGa                                                                                                      >  1:578065/1‑150 (MQ=255)                                         cTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAAt                                                                                                    >  1:614457/1‑151 (MQ=255)                                           tAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAAt                                                                                                    <  2:523885/149‑1 (MQ=255)                                                       tCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTgg                                                                                      >  1:654966/1‑151 (MQ=255)                                                        cATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTgg                                                                                      <  1:4330/150‑1 (MQ=255)                                                         aTAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTGGtt                                                                                    >  2:337817/1‑151 (MQ=255)                                                         aTAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTGGt                                                                                     <  2:249581/150‑1 (MQ=255)                                                            gCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTGGTTGt                                                                                  <  1:104852/150‑1 (MQ=255)                                                              acgacgTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTGGTTGTAAc                                                                               <  2:479329/151‑1 (MQ=255)                                                                  cgTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTGGTTGTAACATc                                                                            >  2:37726/1‑150 (MQ=255)                                                                   gTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTGGTTGTAACATc                                                                            >  2:632251/1‑149 (MQ=255)                                                                         ctTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAAGTTTGTACGCTTTGGTTGTAACATCAGCGACg                                                                     >  1:52560/1‑150 (MQ=255)                                                                                  gACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCTCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTGGTTGTAACATCAGCGACGCTCCTTAttt                                                           <  1:571843/151‑1 (MQ=255)                                                                                  gACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTGGTTGTAACATCAGCGACGCTCCTTa                                                              <  1:318810/148‑1 (MQ=255)                                                                                  gACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTGGTTGTAACATCAGCGACGCTCCTTAt                                                             <  1:592985/149‑1 (MQ=255)                                                                                                    cACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTGGTTGTAACATCAGCGACGCTCCTTATTTACGGCCTTTACgctgctg                                         <  2:667008/151‑1 (MQ=255)                                                                                                          cAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTGGTTGTAACATCAGCGACGCTCCTTATTTACGGCCTTTACGCTGCTGCttttt                                   >  1:233211/1‑151 (MQ=255)                                                                                                           aaCAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTGGTTGTAACATCAGCGACGCTCCTTATTTACGGCCTTTACGCTGCTGCttttt                                   >  1:358130/1‑150 (MQ=255)                                                                                                                    cAGACCGAAGCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTGGTTGTAACATCAGCGACGCTCCTTATTTACGGCCTTTACGCTGCTGCTTTTTAGGCTGagc                          >  2:245652/1‑150 (MQ=255)                                                                                                                     aGACCGAAGCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTGGTTGTAACATCAGCGACGCTCCTTATTTACGGCCTTTACGCTGCTGCTTTTTAGGCTGagc                          >  1:70441/1‑149 (MQ=255)                                                                                                                          gAAGCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTGGTTGTAACATCAGCGACGCTCCTTATTTACGGCCTTTACGCTGCTACTTTTTAGGCTGAGCAGCCGGt                   >  1:596698/1‑151 (MQ=255)                                                                                                                            aGCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTGGTTGTAACATCAGCGACGCTCCTTATTTACGGCCTTTACGCTGCTGCTTTTTAGGCTGAGCAGCCGGttt                 >  2:542451/1‑151 (MQ=255)                                                                                                                             gCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTGGTTGTAACATCAGCGACGCTCCTTATTTACGGCCTTTACGCTGCTGCTTTTTAGGCTGAGCAGCCGGtttt                <  1:144306/151‑1 (MQ=255)                                                                                                                             gCTGCCGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTGGTTGTAACATCAGCGACGCTCCTTATTTACGGCCTTTACGCTGCTGCTTTTTAGGCTGAGCAGCCGGtttt                >  1:284447/1‑151 (MQ=255)                                                                                                                                 ccGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTGGTTGTAACATCAGCGACGCTCCTTATTTACGGCCTTTACGCTGCTGCTTTTTAGGCTGAGCAGCCGGTTTTTCCg            >  2:115566/1‑151 (MQ=255)                                                                                                                                  cGAAGCTAACATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTGGTTGTAACATCAGCGACGCTCCTTATTTACGGCCTTTACGCTGCTGCTTTTTAGGCTGAGCAGCCGGTTTTTCCgg           >  1:569374/1‑151 (MQ=255)                                                                                                                                         aaCATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTGGTTGTAACATCAGCGACGCTCCTTATTTACGGCCTTTACGCTGCTGCTTTTTAGGCTGAGCAGCCGGTTTTTCCGGTTGTTc     >  1:452128/1‑150 (MQ=255)                                                                                                                                          aCATCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTGGTTGTAACATCAGCGACGCTCCTTATTTACGGCCTTTACGCTGCTGCTTTTTAGGCTGAGCAGCCGGTTTTTCCGGTTGTTCa    >  1:231496/1‑150 (MQ=255)                                                                                                                                             tCCGCACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTGGTTGTAACATCAGCGACGCTCCTTATTTACGGCCTTTACGCTGCTGCTTTTTAGGCTGAGCAGCCGGTTTTTCCGGTTGTTCAAc  >  1:357901/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                 GTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTGGTTGTAACATCAGCGACGCTCCTTATTTACGGCCTTTACGCTGCTGCTTTTTAGGCTGAGCAGCCGGTTTTTCCGGTTGTTCAAC  >  NC_000913/3448956‑3449244

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