Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,944,244	C→A	intergenic (+13/‑80)	rrlC → / → rrfC	23S ribosomal RNA/5S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,944,244	0	C	A	100.0%	179.4 / NA	74	intergenic (+13/‑80)	rrlC/rrfC	23S ribosomal RNA/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (36/38);  total (36/38)

GGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGG  >  NC_000913/3944221‑3944333                        |                                                                                          ggCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCtgatg                      <  1:1050279/93‑1 (MQ=14) ggCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCtgatg                      <  1:3435518/93‑1 (MQ=14) ggCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCtgatg                      <  2:3992509/93‑1 (MQ=14)  gCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCtgatga                     >  2:2982216/1‑93 (MQ=17)  gCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCtgatga                     <  2:1853666/93‑1 (MQ=17)  gCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCtgatga                     >  1:2699505/1‑93 (MQ=17)   cTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGaa                    >  1:4119137/1‑93 (MQ=18)    ttAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGaaa                   >  1:1815215/1‑93 (MQ=18)    ttAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGaaa                   <  1:2042638/93‑1 (MQ=18)    ttAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGaaa                   >  2:385494/1‑93 (MQ=18)    ttAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGaaa                   >  2:937338/1‑93 (MQ=18)     tAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAAc                  <  1:2733946/93‑1 (MQ=18)     tAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAAc                  <  2:4234084/93‑1 (MQ=18)     tAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAAc                  >  2:1664248/1‑93 (MQ=18)     tAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAAc                  <  2:2288963/93‑1 (MQ=18)      aaCCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACa                 >  2:2067242/1‑93 (MQ=18)      aaCCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACa                 <  2:2078087/93‑1 (MQ=21)      aaCCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACa                 <  2:136860/93‑1 (MQ=18)      aaCCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACa                 >  2:4191974/1‑93 (MQ=18)      aaCCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACa                 <  2:2406335/93‑1 (MQ=18)      aaCCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACa                 <  2:3283775/93‑1 (MQ=18)      aaCCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACa                 >  1:4335575/1‑93 (MQ=18)      aaCCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACa                 <  2:810504/93‑1 (MQ=18)      aaCCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACa                 <  1:3242376/93‑1 (MQ=18)       aCCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAg                <  2:955174/93‑1 (MQ=18)        ccTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGa               <  1:3648192/93‑1 (MQ=18)        ccTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGa               <  2:2459683/93‑1 (MQ=21)        ccTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGa               >  1:2317758/1‑93 (MQ=18)        ccTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGa               >  2:1616688/1‑93 (MQ=18)        ccTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGa               >  2:2983303/1‑93 (MQ=18)        ccTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGa               <  1:498307/93‑1 (MQ=21)        ccTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGa               >  2:4630886/1‑93 (MQ=18)        ccTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGa               <  1:435136/93‑1 (MQ=21)         cTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGaa              >  2:3408636/1‑93 (MQ=18)         cTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGaa              >  1:2676521/1‑93 (MQ=18)         cTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGaa              >  2:1822031/1‑93 (MQ=18)         cTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGaa              >  2:1499810/1‑93 (MQ=18)         cTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGaa              >  1:4135164/1‑93 (MQ=18)         cTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGaa              >  1:60878/1‑93 (MQ=18)         cTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGaa              <  1:426670/93‑1 (MQ=21)         cTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGaa              >  1:4619479/1‑93 (MQ=18)         cTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGaa              >  1:4607543/1‑93 (MQ=18)          ttACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAAt             >  1:1932375/1‑93 (MQ=18)          ttACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAAt             <  1:3347364/93‑1 (MQ=21)          ttACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAAt             >  1:2924603/1‑93 (MQ=18)          ttACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAAt             >  1:2778749/1‑93 (MQ=18)           tACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAAtt            <  2:11645/93‑1 (MQ=21)            aCAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAAttt           <  1:2253695/93‑1 (MQ=21)            aCAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAAttt           <  1:181538/93‑1 (MQ=21)             cAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTg          <  2:1728921/93‑1 (MQ=18)              aaCGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGc         >  1:4184596/1‑93 (MQ=18)              aaCGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGc         >  1:3943799/1‑93 (MQ=18)              aaCGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGc         <  1:2598852/93‑1 (MQ=21)               aCGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGcc        <  1:3056697/93‑1 (MQ=21)               aCGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGcc        <  2:670221/93‑1 (MQ=18)               aCGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGcc        <  1:41893/93‑1 (MQ=21)                cGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCt       <  1:3389919/93‑1 (MQ=18)                cGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCt       <  1:839202/93‑1 (MQ=18)                cGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCt       <  2:3379717/93‑1 (MQ=21)                cGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCt       <  1:3696253/93‑1 (MQ=21)                cGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCt       <  2:3868621/93‑1 (MQ=21)                 gCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTg      <  1:1154739/93‑1 (MQ=18)                 gCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTg      <  2:3994569/93‑1 (MQ=21)                 gCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTg      <  1:860359/93‑1 (MQ=21)                  ccGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTgg     >  2:491838/1‑93 (MQ=18)                   cGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTggc    >  2:3555539/1‑93 (MQ=18)                   cGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTggc    <  1:1887922/93‑1 (MQ=21)                    gAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTggcg   >  1:4540938/1‑93 (MQ=18)                    gAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTggcg   >  1:2922516/1‑93 (MQ=18)                     aaGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTggcgg  >  2:341818/1‑93 (MQ=18)                     aaGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTggcgg  >  1:558757/1‑93 (MQ=18)                     aaGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTggcgg  >  1:3104526/1‑93 (MQ=18)                     aaGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTggcgg  >  2:82723/1‑93 (MQ=18)                     aaGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTggcgg  >  1:286363/1‑93 (MQ=18)                        |                                                                                          GGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGG  >  NC_000913/3944221‑3944333

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