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%	159.4 / NA	66	intergenic (+13/‑80)	rrlC/rrfC	23S ribosomal RNA/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (33/33);  total (33/33)

GGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGG  >  NC_000913/3944221‑3944333                        |                                                                                          ggCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCtgatg                      >  2:179311/1‑93 (MQ=14) ggCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCtgatg                      >  1:5672532/1‑93 (MQ=14)  gCTTCACCTTACAACGCCAAAGATGTTTTGGCGGCTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGCAGCGGTCtgatga                     >  1:3754669/1‑93 (MQ=11)  gCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCtgatga                     <  1:5279080/93‑1 (MQ=17)   cTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGaa                    >  1:2002565/1‑93 (MQ=18)   cTTAACCTTACAACGCCGAAGATGCTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGTGGTCTGATGaa                    >  1:3186733/1‑93 (MQ=255)    ttAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGaaa                   <  2:2282018/93‑1 (MQ=18)    ttAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGaaa                   <  1:4805217/93‑1 (MQ=18)     tAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAAc                  <  2:6119208/93‑1 (MQ=18)     tAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAAc                  >  2:5349952/1‑93 (MQ=18)     tAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAAc                  >  2:235985/1‑93 (MQ=18)      aaCCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACa                 <  2:4947261/93‑1 (MQ=21)      aaCCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACa                 <  1:4779849/93‑1 (MQ=21)      aaCCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACa                 >  1:144877/1‑93 (MQ=18)      aaCCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACa                 <  1:4606073/93‑1 (MQ=21)       aCCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAg                >  1:1410333/1‑93 (MQ=18)       aCCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAg                <  1:6403621/93‑1 (MQ=18)       aCCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAg                >  1:1032543/1‑93 (MQ=18)       aCCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAg                >  1:6331088/1‑93 (MQ=18)       aCCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAg                <  1:6670843/93‑1 (MQ=18)       aCCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAg                >  2:1754631/1‑93 (MQ=18)        ccTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGa               <  2:2689594/93‑1 (MQ=21)        ccTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGa               >  2:1962477/1‑93 (MQ=18)         cTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGaa              <  1:3636654/93‑1 (MQ=21)         cTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGaa              >  1:3436174/1‑93 (MQ=18)         cTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGaa              <  2:1260263/93‑1 (MQ=21)         cTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGaa              >  1:171956/1‑93 (MQ=18)         cTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGaa              >  2:4973660/1‑93 (MQ=18)         cTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGaa              >  2:2915266/1‑93 (MQ=18)          ttACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAAt             >  1:4319030/1‑93 (MQ=18)          ttACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAAt             <  2:6364762/93‑1 (MQ=18)           tACAACGCCGAAGATGTGTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAAtt            <  2:3475243/93‑1 (MQ=16)            aCAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAAttt           <  1:3671095/93‑1 (MQ=21)            aCAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAAttt           <  2:342703/93‑1 (MQ=21)            aCAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAAttt           >  1:3705861/1‑93 (MQ=18)             cAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTg          >  1:2163933/1‑93 (MQ=18)             cAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTg          >  1:2034307/1‑93 (MQ=18)             cAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTg          <  1:1683985/93‑1 (MQ=18)             cAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTg          >  2:6676664/1‑93 (MQ=18)              aaCGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGc         >  1:6520671/1‑93 (MQ=18)              aaCGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGc         >  2:5300880/1‑93 (MQ=18)              aaCGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGc         >  2:2191854/1‑93 (MQ=18)               aCGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGcc        <  2:814007/93‑1 (MQ=21)               aCGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGcc        <  2:766634/93‑1 (MQ=21)               aCGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGcc        >  2:6363972/1‑93 (MQ=18)               aCGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGcc        >  2:6023962/1‑93 (MQ=18)               aCGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGcc        <  1:5536832/93‑1 (MQ=21)               aCGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGcc        <  2:4484063/93‑1 (MQ=18)               aCGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGcc        <  2:3987494/93‑1 (MQ=21)                cGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCt       <  1:3327473/93‑1 (MQ=21)                cGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCt       <  1:2865636/93‑1 (MQ=18)                cGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCt       <  2:4459571/93‑1 (MQ=18)                cGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCt       <  1:5673463/93‑1 (MQ=21)                cGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCt       <  2:2303491/93‑1 (MQ=18)                 gCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTg      <  1:3279466/93‑1 (MQ=21)                 gCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTg      <  2:5678446/93‑1 (MQ=18)                 gCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTg      <  2:2479069/93‑1 (MQ=21)                 gCCGAAGATGTTGTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATGTGCCTg      <  1:1501603/93‑1 (MQ=255)                  ccGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTgg     <  1:5455522/93‑1 (MQ=21)                  ccGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTgg     <  2:2018969/93‑1 (MQ=18)                  ccGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGACTGATGAAACAGAATTTGCCTgg     >  2:1132483/1‑93 (MQ=16)                   cGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTggc    >  2:805966/1‑93 (MQ=18)                    gAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTggcg   <  1:5525089/93‑1 (MQ=21)                    gAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTggcg   >  2:5613086/1‑93 (MQ=18)                     aaGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTggcgg  >  1:6170706/1‑93 (MQ=18)                     aaGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTggcgg  >  1:2853460/1‑93 (MQ=18)                     aaGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTggcgg  >  1:2059674/1‑93 (MQ=18)                        |                                                                                          GGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGG  >  NC_000913/3944221‑3944333

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