Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,645,285	A→G	intergenic (‑103/‑100)	prlC ← / → rlmJ	oligopeptidase A/23S rRNA m(6)A2030 methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,645,285	0	A	G	100.0%	130.2 / NA	39	intergenic (‑103/‑100)	prlC/rlmJ	oligopeptidase A/23S rRNA m(6)A2030 methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/27);  total (12/27)

AGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATGAGAAACGTTACGCTTTCGCGGCGGCGGCCTCTTTTCGGTATACTGTCCTGATACGCTTTTGTGCGCCCCGAATACGGGCCGATTTTT  >  NC_000913/3645209‑3645362                                                                             |                                                                              aGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATGAGAAACGTTgcgc                                                                            <  1:640954/80‑1 (MQ=255) aGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATGAGAAACGTTGCGCt                                                                           <  2:640952/81‑1 (MQ=255)    aaGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATGAGAAACGTTGCGCTTTc                                                                        <  2:640955/81‑1 (MQ=255)    aaGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATGAGAAACGTTGCGCTTTc                                                                        <  2:640954/81‑1 (MQ=255)    aaGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATGAGAAACGTTGCGCTTTc                                                                        >  1:640999/1‑81 (MQ=255)     aGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATGAGAAACGTTGCGCTTTcg                                                                       >  1:641003/1‑81 (MQ=255)        gcTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTCAAGTGTAGTGAATTTCAATGAGAAACGTTGCGCTTTCgcgg                                                                    <  2:640956/81‑1 (MQ=255)                 gcgcTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATGAGAAACGTTGCGCTTTCGCGGCGGCGGCct                                                           <  1:640958/81‑1 (MQ=255)                   gcTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATGAGAAACGTTGCGCTTTCGCGGCGGCGGCctct                                                         <  1:640959/81‑1 (MQ=255)                   gcTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATGAGAAACGTTGCGCTTTCGCGGCGGCGGCctct                                                         >  2:641003/1‑81 (MQ=255)                   gcTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATGAGAAACGTTGCGCTTTCGCGGCGGCGGCctct                                                         >  1:641006/1‑81 (MQ=255)                   gcTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATGAGAAACGTTGCGCTTTCGCGGCGGCGGCctct                                                         <  1:640960/81‑1 (MQ=255)                    cTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATGAGAAACGTTGCGCTTTCGCGGCGGCGGCCTCtt                                                        >  1:641032/1‑81 (MQ=255)                       tACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATGAGAAACGTTGCGCTTTCGCGGCGGCGGCCTCTTTTc                                                     >  2:641013/1‑81 (MQ=255)                       tACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATGAGAAACGTTGCGCTTTCGCGGCGGCGGCCTCTTTTc                                                     >  1:641018/1‑81 (MQ=255)                        aCATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATGAGAAACGTTGCGCTTTCGCGGCGGCGGCCTCTTTTc                                                     <  1:640961/80‑1 (MQ=255)                         cATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATGAGAAACGTTGCGCTTTCGCGGCGGCGGCCTCTTTTCgg                                                   <  1:640962/81‑1 (MQ=255)                          aTAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATGAGAAACGTTGCGCTTTCGCGGCGGCGGCCTCTTTTCGGt                                                  <  2:640962/81‑1 (MQ=255)                             aTGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATGAGAAACGTTGCGCTTTCGCGGCGGCGGCCTCTTTTCGGtat                                                <  1:640964/80‑1 (MQ=255)                                tgtgtAGCATGGGGTTAAGTGTAGTGAATTTCAATGAGAAACGTTGCGCTTTCGCGGCGGCGGCCTCTTTTCGGTATACt                                             <  1:640965/80‑1 (MQ=255)                                      gCATGGGGTTAAGTGTAGTGAATTTCAATGAGAAACGTTGCGCTTTCGCGGCGGCGGCCTCTTTTCGGTATACTGTCCTGa                                      <  1:640966/81‑1 (MQ=255)                                        aTGGGGTTAAGTGTAGTGAATTTCAATGAGAAACGTTGCGCTTTCGCGGCGGCGGCCTCTTTTCGGTATACTGTCCTGGTa                                    >  2:641027/1‑81 (MQ=255)                                                aaGTGTAGTGAATTTCAATGAGAAACGTTGCGCTTTCGCGGCGGCGGCCTCTTTTCGGTATACTGTCCTGATACGCtttt                             <  2:640966/80‑1 (MQ=255)                                                   tgtAGTGAATTTCAATGAGAAACGTTGCGCTTTCGCGGCGGCGGCCTCTTTTCGGTATACTGTCCTGATACGCTTTTGTgc                         <  2:640967/81‑1 (MQ=255)                                                   tgtAGTGAATTTCAATGAGAAACGTTGCGCTTTCGCGGCGGCGGCCTCTTTTCGGTATACTGTCCTGATACGCTTTTGTgc                         <  2:640968/81‑1 (MQ=255)                                                     tAGTGAATTTCAATGAGAAACGTTGCGCTTTCGCGGCGGCGGCCTCTTTTCGGTATACTGTCCTGATACGCTTTTGTgcg                        <  2:640969/80‑1 (MQ=255)                                                          aaTTTCAATGAGAAACGTTGCGCTTTCGCGGCGGCGGCCTCTTTTCGGTATACTGTCCTGATACGCTTTTGTGCGCCCCg                   <  2:640970/80‑1 (MQ=255)                                                           aTTTCAATGAGAAACGTTGCGCTTTCGCGGCGGCGGCCTCTTTTCGGTATACTGTCCTGATACGCTTTTGTGCGCCCCGaa                 <  1:640972/81‑1 (MQ=255)                                                            tttCAATGAGAAACGTTGCGCTTTCGCGGCGGCGGCCTCTTTTCGGTATACTGTCCTGATACGCTTTTGTGCGCCCCGAAt                <  2:640972/81‑1 (MQ=255)                                                             ttCAATGAGAAACGTTGCGCTTTCGCGGCGGCGGCCTCTTTTCGGTATACTGTCCTGATACGCTTTTGTGCGCCCCGAATa               >  2:641021/1‑81 (MQ=255)                                                               cAATGAGAAACGTTGCGCTTTCGCGGCGGCGGCCTCTTTTCGGTATACTGTCCTGATACGCTTTTGTGCGCCCCGAATAc              <  2:640973/80‑1 (MQ=255)                                                                aaTGAGAAACGTTGCGCTTTCGCGGCGGCGGCCTCTTTTCGGTATACTGTCCTGATACGCTTTTGTGCGCCCCGAATACgg            >  1:641027/1‑81 (MQ=255)                                                                aaTGAGAAACGTTGCGCTTTCGCGGCGGCGGCCTCTTTTCGGTATACTGTCCTGATACGCTTTTGTGCGCCCCGAATACgg            <  2:640974/81‑1 (MQ=255)                                                                  tGAGAAACGTTGCGCTTTCGCGGCGGCGGCCTCTTTTCGGTATACTGTCCTGATACGCTTTTGTGCGCCCCGAATACGGGc          <  1:640976/81‑1 (MQ=255)                                                                    agaAACGTTGCGCTTTCGCGGCGGCGGCCTCTTTTCGGTATACTGTCCTGATACGCTTTTGTGCGCCCCGAATACGGGcc         >  2:641033/1‑80 (MQ=255)                                                                    agaAACGTTGCGCTTTCGCGGCGGCGGCCTCTTTTCGGTATACTGTCCTGATACGCTTTTGTGCGCCCCGAATACGGGCCg        <  2:640976/81‑1 (MQ=255)                                                                      aaaCGTTGCGCTTTCGCGGCGGCGGCCTCTTTTCGGTATACTGTCCTGATACGCTTTTGTGCGCCCCGAATACGGGCCGAt      <  1:640978/81‑1 (MQ=255)                                                                        aCGTTGCGCTTTCGCGGCGGCGGCCTCTTTTCGGTATACTGTCCTGATACGCTTTTGTGCGCCCCGAATACGGGCCGAttt    <  1:640979/81‑1 (MQ=255)                                                                          gTTGCGCTTTCGCGGCGGCGGCCTCTTTTCGGTATACTGTCCTGATACGCTTTTGTGCGCCCCGAATACGGGCCGAttttt  <  2:640979/81‑1 (MQ=255)                                                                          gTTGCGCTTTCGCGGCGGCGGCCTCTTTTCGGTATACTGTCCTGATACGCTTTTGTGCGCCCCGAATACGGGCCGAtttt   >  2:641020/1‑80 (MQ=255)                                                                             |                                                                              AGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATGAGAAACGTTACGCTTTCGCGGCGGCGGCCTCTTTTCGGTATACTGTCCTGATACGCTTTTGTGCGCCCCGAATACGGGCCGATTTTT  >  NC_000913/3645209‑3645362

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