Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,428,954	G→A	intergenic (‑192/‑282)	rrsD ← / → yrdA	16S ribosomal RNA of rrnD operon/bacterial transferase hexapeptide domain protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,428,954	0	G	A	100.0%	115.7 / NA	34	intergenic (‑192/‑282)	rrsD/yrdA	16S ribosomal RNA of rrnD operon/bacterial transferase hexapeptide domain protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/18);  total (16/18)

TATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAATTTTTTGCAC  >  NC_000913/3428860‑3429081                                                                                               |                                                                                                                                tATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTtct                                                                                     >  1:105371/1‑139 (MQ=18)      tCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCg                                                                                <  2:320478/139‑1 (MQ=37)      tCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCg                                                                                <  1:953665/139‑1 (MQ=37)             cccACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTcc                                                                         >  1:612398/1‑139 (MQ=255)               cacacaGATTGTCTGATAAATAGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGa                                                                       >  1:673718/1‑139 (MQ=255)                  acaGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCa                                                                    >  1:925535/1‑139 (MQ=255)                      aTTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCt                                                                >  1:1090341/1‑139 (MQ=255)                       ttGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCgttgt                                                                            >  2:931254/1‑126 (MQ=255)                       ttGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCgttgt                                                                            <  1:931254/126‑1 (MQ=255)                            tGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAg                                                          >  1:1083422/1‑139 (MQ=255)                              aTAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGaa                                                           <  2:94939/136‑1 (MQ=255)                              aTAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGaa                                                           >  1:94939/1‑136 (MQ=255)                                    tGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCAc                                                  >  1:73060/1‑139 (MQ=255)                                    tGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCAc                                                  >  1:143332/1‑139 (MQ=255)                                      ttAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCg                                                                                >  1:943322/1‑107 (MQ=255)                                      ttAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCg                                                                                <  2:943322/107‑1 (MQ=255)                                            agCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTgtcgt                                          <  2:143332/139‑1 (MQ=255)                                              cAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATc                                                                  <  1:442687/113‑1 (MQ=255)                                              cAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATc                                                                  >  2:442687/1‑113 (MQ=255)                                                gTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATc                                                                  >  2:1085420/1‑111 (MQ=255)                                                gTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATc                                                                  <  1:1085420/111‑1 (MQ=255)                                                 ttGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCaa                                     <  2:1090341/139‑1 (MQ=255)                                                      aCGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAg                                <  2:73060/139‑1 (MQ=255)                                                      aCGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAg                                <  1:721922/139‑1 (MQ=255)                                                       cgcgGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAgg                               <  2:105371/139‑1 (MQ=255)                                                       cgcgGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAgg                               <  1:242903/139‑1 (MQ=255)                                                       cgcgGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAgg                               <  2:983522/139‑1 (MQ=255)                                                        gcgGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGgc                              <  2:525098/139‑1 (MQ=255)                                                         cgGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGgcg                             <  2:787565/139‑1 (MQ=255)                                                         cgGCTTTCAGCTCACTGTCGCGAGGTGGCATATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGgcg                             <  2:188165/139‑1 (MQ=255)                                                           gCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCa                           >  2:71283/1‑139 (MQ=255)                                                                        tGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATccc              >  2:380107/1‑139 (MQ=255)                                                                                 gTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAATTTTTTg     <  2:612398/139‑1 (MQ=255)                                                                                    gCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAATTTTTTGcac  >  2:315501/1‑139 (MQ=255)                                                                                               |                                                                                                                                TATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAATTTTTTGCAC  >  NC_000913/3428860‑3429081

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