Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,408,885	(A)6→5	coding (2232/2442 nt)	rnr →	exoribonuclease R, RNase R

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,408,880	0	A	.	100.0%	210.0 / NA	43	coding (2227/2442 nt)	rnr	exoribonuclease R, RNase R
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (18/25);  total (18/25)

GGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCGAAAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGCGTCGTCAGGTCGGTAAAAAGGTAAACTTTGAGCCAGACAGCGCCTTCCGCGGTGAGAAAAAAACGAAGCCGAAAGCGGCGAAGAAAGACGCGAGAAAAGCGAAAAAGCCAT  >  NC_000913/4408730‑4409024                                                                                                                                                       |                                                                                                                                                 ggggAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCGa                                                                                                                                                  <  2:310632/151‑1 (MQ=255)  gggAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCGa                                                                                                                                                  >  2:341046/1‑150 (MQ=255)    gAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCGaaaa                                                                                                                                               <  1:627167/151‑1 (MQ=255)          aGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGa                                                                                                                                        <  2:243788/151‑1 (MQ=255)              gCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCa                                                                                                                                    <  1:342585/151‑1 (MQ=255)                cAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAgg                                                                                                                                  >  2:505276/1‑151 (MQ=255)                       aTCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTaaaaa                                                                                                                            <  1:133075/150‑1 (MQ=255)                       aTCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTaaaaa                                                                                                                            >  2:601700/1‑150 (MQ=255)                        tCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTaaaaa                                                                                                                            >  1:762341/1‑149 (MQ=255)                         cGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAgg                                                                                                                         >  2:574568/1‑151 (MQ=255)                              gggCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCTCCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCa                                                                                                                    <  2:386870/151‑1 (MQ=255)                                gCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAg                                                                                                                  <  2:500567/151‑1 (MQ=255)                                 cGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGc                                                                                                                 <  1:620081/151‑1 (MQ=255)                                   aTCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGcg                                                                                                                <  1:71683/150‑1 (MQ=255)                                    tCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAggc                                                                                                                        >  1:38248/1‑141 (MQ=255)                                    tCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAggc                                                                                                                        <  2:38247/141‑1 (MQ=255)                                    tCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGcgtc                                                                                                              >  1:763054/1‑151 (MQ=255)                                      gcgTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGcgtcgt                                                                                                            <  1:289687/151‑1 (MQ=255)                                      gcgTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGcgtcgt                                                                                                            >  1:693668/1‑151 (MQ=255)                                        gTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGCGTCGTCa                                                                                                          <  1:583057/151‑1 (MQ=255)                                              gTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGCGTCGTCAGGTCgg                                                                                                    >  1:112313/1‑151 (MQ=255)                                                   cgTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGCGTCGTCAGGTCGGTaaa                                                                                                >  2:755088/1‑150 (MQ=255)                                                   cgTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGCGTCGTCAGGTCGGTaaa                                                                                                >  1:541528/1‑150 (MQ=255)                                                    gTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGCGTCGTCAGGTCGGTaaaaa                                                                                              <  1:329522/151‑1 (MQ=255)                                                      cGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGTAAGCGTCGTCAGGTCGGTAAAAAgg                                                                                            <  2:588666/151‑1 (MQ=255)                                                       gAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGCGTCGTCAGGTCGGTAAAAAgg                                                                                            >  1:169803/1‑150 (MQ=255)                                                         aGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGCGTCGTCAGGTCGGTAAAAAGGTa                                                                                          >  1:670030/1‑150 (MQ=255)                                                            ggTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGCGTCGTCAGGTCGGTAAAAAGGTAAACt                                                                                      <  1:266002/151‑1 (MQ=255)                                                                aaTATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGCGTCGTCAGGTCGGTAAAAAGGTAAACTTTGa                                                                                  <  1:482178/151‑1 (MQ=255)                                                                   atGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGCGTCGTCAGGTCGGTAAAAAGGTAAACTTTGAGcc                                                                               >  2:535493/1‑151 (MQ=255)                                                                          aGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGCGTCGTCAGGTCGGTAAAAAGGTAAACTTTGAGCCAGACAgc                                                                        <  2:462456/151‑1 (MQ=255)                                                                              cAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGCGTCGTCAGGTCGGTAAAAAGGTAAACTTTGAGCCAGACAGCGCCt                                                                    >  2:697117/1‑151 (MQ=255)                                                                              cAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGCGTCGTCAGGTCGGTAAAAAGGTAAACTTTGAGCCAGACAGCGCCt                                                                    <  1:442230/151‑1 (MQ=255)                                                                                aaaTCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGCGTCGTCAGGTCGGTAAAAAGGTAAACTTTGAGCCAGACAGCGCCTTc                                                                  <  1:35104/151‑1 (MQ=255)                                                                                   tCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGCGTCGTCAGGTCGGTAAAAAGGTAAACTTTGAGCCAGACAGCGCCTTCcgc                                                               >  1:328573/1‑151 (MQ=255)                                                                                                 aTCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGCGTCGTCAGGTCGGTAAAAAGGTAAACTTTGAGCCAGACAGCGCCTTCCGCGGTGAGAAAAAAAc                                                 >  1:706092/1‑151 (MQ=255)                                                                                                      cAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGCGTCGTCAGGTCGGTAAAAAGGTAAACTTTGAGCCAGACAGCGCCTTCCGCGGTGAGAAAAAAAc                                                 <  2:166072/146‑1 (MQ=255)                                                                                                        gCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGCGTCGTCAGGTCGGTAAAAAGGTAAACTTTGAGCCAGACAGCGCCTTCCGCGGTGAGAAAAAAACGAAGCCg                                          <  1:441593/151‑1 (MQ=255)                                                                                                            aCGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGCGTCGTCAGGTCGGTAAAAAGGTAAACTTTGAGCCAGACAGCGCCTTCCGCGGTGAGAAAAAAACGAAGCCGAAAg                                      <  1:261945/151‑1 (MQ=255)                                                                                                                     gcgcAACGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGCGTCGTCAGGTCGGTAAAAAGGTAAACTTTGAGCCAGACAGCGCCTTCCGCGGTGAGAAAAAAACGAAGCCGAAAGCGGCgaaga                             <  2:367321/151‑1 (MQ=255)                                                                                                                         aaCGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAAGCGGCAAGCGTCGTCAGGTCGGTAAAAAGGTAAACTTTGAGCCAGACAGCGCCTTCCGCGGTGAGAAAAAAACGAAGCCGAAAGCGGCgaagaa                            >  2:710017/1‑148 (MQ=255)                                                                                                                          aCGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGCGTCGTCAGGTCGGTAAAAAGGTAAACTTTGAGCCAGACAGCGCCTTCCGCGGTGAGAAAAAAACGAAGCCGAAAGCGGCGAAGAAAGAc                        <  1:227584/151‑1 (MQ=255)                                                                                                                          aCGTCGGTAAAACGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGCGTCGTCAGGTCGGTAAAAAGGTAAACTTTGAGCCAGACAGCGCCTTCCGCGGTGAGAAAAAAACGAAGCCGAAAGCGGCGAAGAAAGAc                        <  2:293093/151‑1 (MQ=255)                                                                                                                                     aCGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGCGTCGTCAGGTCGGTAAAAAGGTAAACTTTGAGCCAGACAGCGCCTTCCGCGGTGAGAAAAAAACGAAGCCGAAAGCGGCGAAGAAAGACGCGAGAAAAGc             >  2:727632/1‑151 (MQ=255)                                                                                                                                      cGGCGCGCGAGAAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGCGTCGTCAGGTCGGTAAAAAGGTAAACTTTGAGCCAGACAGCGCCTTCCGCGGTGAGAAAAAAACGAAGCCGAAAGCGGCGAAGAAAGACGCGAGAAAAGCg            <  1:666193/151‑1 (MQ=255)                                                                                                                                                gaAAGCG‑AAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGCGTCGTCAGGTCGGTAAAAAGGTAAACTTTGAGCCAGACAGCGCCTTCCGCGGTGAGAAAAAAACGAAGCCGAAAGCGGCGAAGAAAGACGCGAGAAAAGCGAAAAAGCCAt  <  1:544985/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                 GGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCGAAAAAAGGCGATGCAGGTAAAAAAGGCGGCAAGCGTCGTCAGGTCGGTAAAAAGGTAAACTTTGAGCCAGACAGCGCCTTCCGCGGTGAGAAAAAAACGAAGCCGAAAGCGGCGAAGAAAGACGCGAGAAAAGCGAAAAAGCCAT  >  NC_000913/4408730‑4409024

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