Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,210,838	A→G	K20E (AAG→GAG)	rpsU →	30S ribosomal subunit protein S21

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,210,838	0	A	G	100.0%	136.2 / NA	44	K20E (AAG→GAG)	rpsU	30S ribosomal subunit protein S21
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/22);  total (22/22)

TGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGAAAACGC  >  NC_000913/3210700‑3210974                                                                                                                                           |                                                                                                                                         tGCAGTAATATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAACCGCGCCGTTCGACGTAGCTCTGCGTCGCTTCGAGCg                                                                                                                                      <  2:753362/143‑1 (MQ=255)     gTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAATCCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCGAGCGTTcc                                                                                                                                  <  2:1107725/143‑1 (MQ=255)      tAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCGTTTTTCGCCAATAAAAGTACGTGAAAACGAGACATTAGAAGTAGCAACGCGTCGCTTCGAGCGttggt                                                                                                                                 >  1:905917/1‑140 (MQ=255)      tAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCGAGCGTTCCt                                                                                                                                 >  1:1141616/1‑143 (MQ=255)              cGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCGAGCGTTCCTGAGaaaaa                                                                                                                         >  1:221021/1‑143 (MQ=255)               gCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCGAGCGTTCCTGCGAAAAAg                                                                                                                        <  1:761190/143‑1 (MQ=255)                 aCCATGTTGAAACAAGCTGGCGGCGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCGAGCGTTCCTGCGAAAAAGCa                                                                                                                      <  1:922731/143‑1 (MQ=255)                  ccATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCGAGCGTTCCTGCGAAAAAGCAg                                                                                                                     <  2:320781/143‑1 (MQ=255)                                 cTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAg                                                                                                      >  1:711114/1‑143 (MQ=255)                                            gCCGGCGGCAAACCGAATTAAGCAAAGGTGAGAGGCACAGGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTcgtcgtcgt                                                                                           <  1:831154/143‑1 (MQ=255)                                              cAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGa                                                                                         <  1:802911/143‑1 (MQ=255)                                                          ggaTTAATCAAAGGTGAGAGGCACATGCCGGTAATGAAAGTACGGGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGaaaa                                                                             <  2:46959/141‑1 (MQ=255)                                                          gAATTAATCAAAGGGGAGAGGCACATGCCGGTAATTAAAGTACGTGGAGACGAGCCGTGCGACGTAGCTCTGCGTCGCGTCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGaaaa                                                                             <  1:569698/143‑1 (MQ=255)                                                           aaTTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGaaaaa                                                                            <  1:843454/143‑1 (MQ=255)                                                            atgaaTCGAAGGTGAGGGGCACATGCCGGTAATTAAAGTACGTGAAGACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAAc                                                                           <  2:1037205/140‑1 (MQ=255)                                                             ttAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAAcc                                                                          <  1:841379/143‑1 (MQ=255)                                                               aaTCGAAGGTGAGAGGCACATGGCGGTAATTAAAGTACGTGACAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGa                                                                        <  1:885654/143‑1 (MQ=255)                                                                aTCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCgcc                                                                       >  1:181097/1‑141 (MQ=255)                                                                 tCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCAAcc                                                                      >  1:273255/1‑142 (MQ=255)                                                                   aaaGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTAc                                                                    >  1:1002136/1‑143 (MQ=255)                                                                   aaaGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTAc                                                                    >  2:907069/1‑143 (MQ=255)                                                                   aaaGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGCTCGTCGTCGCGAGTTCTATGAAAAACCGACTAc                                                                    >  2:64521/1‑143 (MQ=255)                                                                   aaaGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTGCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGCTCGTCGTCGTGAGTTCTATGAAAAAGCGAGTAc                                                                    >  2:1127462/1‑143 (MQ=255)                                                                       gTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGaa                                                                <  1:599936/143‑1 (MQ=255)                                                                       gTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGGTCGTCGTCGTGGGTTCTATGAAAAACCGACTACCGaa                                                                >  1:775498/1‑143 (MQ=255)                                                                                  aTGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAACCCCGCCTACCGATCTTAAgctctc                                                     >  1:998733/1‑139 (MQ=255)                                                                                      cgttAATTAAAGTACGTGAGAGCGAGCCGTTCGACGTAGCTCTGCGTCGCCTCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTaaa                                                 <  2:488752/140‑1 (MQ=255)                                                                                        gTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCGAGCGTTCCTGCGCAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGATAAACCGCGTACCCAACTTAACAGCGCTAAAgg                                               >  2:870741/1‑142 (MQ=255)                                                                                            ttAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTATGCGTCGCTTCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTACAGCTTCt                                           >  1:755525/1‑143 (MQ=255)                                                                                              acaGTACGTGAACACGAGCCGATCGACGTAGCTCTGCGTCGCTTCGAGCGTTCCTGCGAAAAAGCCGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGCCTACCGACCGTAAGCGCGCCAAACCtacccc                                         >  2:1089674/3‑138 (MQ=255)                                                                                                   aCGTCAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGa                                    <  1:530127/143‑1 (MQ=255)                                                                                                       gAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGGTCTATGAAAAACCGACTAGCGAACGTAAGCGCGCTAAAGCGTCTGCAGTGAGACg                                >  2:176383/1‑143 (MQ=255)                                                                                                           gcGAGCTGTTCGGCGTAGCTCCTCGTCGCTTCGAGCGTTCCTGCGAAAAAGCAGTTGTTCTGGCGGAATTTCGTCGTCGTGAGTTCTATGAAATACCGACTACCGAATGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCAc                            <  2:579701/142‑1 (MQ=255)                                                                                                             gAGCCGTGCGACGTAGCTCTGCGTCGCTTCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCAcgc                          <  1:374422/143‑1 (MQ=255)                                                                                                                 cGTTCGACGTAGCTCTGCGTCGCTTCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCgaag                      >  1:445293/1‑143 (MQ=255)                                                                                                                   ttCGACGTAGCTCTGCGTCGCTTCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCCAAAGCTTCTGCAGTGGAACGGCACGCgacgaa                    >  2:578881/1‑143 (MQ=255)                                                                                                                       aCGTAGCTCTGCGTCGCTTCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTg                <  2:653342/143‑1 (MQ=255)                                                                                                                        cGTAGCTCTGCGTCGCGTCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAGCGTCACGCGAAGAAACTgg               <  2:882582/143‑1 (MQ=255)                                                                                                                           aGCTCTGCGTCGCTTCGAGCGTTCCTGCGAAAAAGCCGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGACAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTc            <  2:958152/143‑1 (MQ=255)                                                                                                                            gCTCTGCGGCGCTTCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCGAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTcg           >  2:948888/1‑143 (MQ=255)                                                                                                                             ctctGCGTCGCTTCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAACAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCCGCAGTGAAACGTCACGCGAAGGAACTGGCTcgc          >  2:212567/1‑143 (MQ=255)                                                                                                                              tctGCGTCGCTTCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGGAAAACCGACTACCGAACGTAAGCGCGCTACAGCTTCTGCAGTGAAACGTCACGCGACGAAACTGGCTcgcg         >  1:33583/1‑143 (MQ=255)                                                                                                                               ctGCGTCGCTTCGAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGa        <  2:1230046/143‑1 (MQ=255)                                                                                                                                     cGCTTCGAGCGTTCCTGCGAAGAAGCAAATAGCCCTCCGGAAGTGCGTCGTCGTGAGTTCTATGAAAAACCGACCACGAAACATAAGCGCGCTAAATCTTCAGCGGTGAAACGTCACGCGAAGGAACTGGCTCGCGGGAacgc  >  2:942915/1‑143 (MQ=255)                                                                                                                                           |                                                                                                                                         TGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGAAAACGC  >  NC_000913/3210700‑3210974

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