Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,441,538	C→T	*213* (TAG→TAA)	msrA ←	methionine sulfoxide reductase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,441,538	0	C	T	100.0%	190.8 / NA	58	*213* (TAG→TAA)	msrA	methionine sulfoxide reductase A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/29);  total (29/29)

CCGATCAGGCCTACGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGA  >  NC_000913/4441408‑4441668                                                                                                                                   |                                                                                                                                   ccGATCAGGCCTACGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAATTCTTTTACGCCGCATCAGGCATAGTCAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTc                                                                                                                            <  1:924097/139‑1 (MQ=255)              cGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCggcggc                                                                                                                     >  2:1381983/1‑133 (MQ=255)              cGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCggcggc                                                                                                                     <  1:1381983/133‑1 (MQ=255)                   tCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgcc                                                                                                          >  1:784033/1‑139 (MQ=255)                         cAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattcc                                                                                                    >  2:428866/1‑139 (MQ=255)                                 tGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGcagaca                                                                                                                <  1:367777/119‑1 (MQ=255)                                 tGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGcagaca                                                                                                                >  2:367777/1‑119 (MQ=255)                                 tGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaattc                                                                                            <  2:784033/139‑1 (MQ=255)                                    aTTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcg                                                                                                            <  2:1319171/120‑1 (MQ=255)                                    aTTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcg                                                                                                            >  1:1319171/1‑120 (MQ=255)                                        ggACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTa                                                                                     >  2:476864/1‑139 (MQ=255)                                        ggACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTa                                                                                     >  1:633742/1‑139 (MQ=255)                                                   tAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGt                                                                          >  2:1202936/1‑139 (MQ=255)                                                           gTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaa                                                                                               >  2:532818/1‑110 (MQ=255)                                                           gTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaa                                                                                               <  1:532818/110‑1 (MQ=255)                                                           gTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgcc                                                                                                 >  1:847290/1‑108 (MQ=255)                                                           gTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgcc                                                                                                 <  2:847290/108‑1 (MQ=255)                                                            tAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCa                                                                 <  2:527989/139‑1 (MQ=255)                                                                          cgcATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCtggtgg                                                   >  2:715285/1‑139 (MQ=255)                                                                               cAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATc                                              <  1:428866/139‑1 (MQ=255)                                                                                 ggCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCtgg                                                      >  2:1332575/1‑129 (MQ=255)                                                                                 ggCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCtgg                                                      <  1:1332575/129‑1 (MQ=255)                                                                                      aGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgc                                                                                                  >  1:1297829/1‑80 (MQ=255)                                                                                      aGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgc                                                                                                  <  2:1297829/80‑1 (MQ=255)                                                                                            agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGcaga                                                                                                                  <  1:639119/58‑1 (MQ=255)                                                                                            agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGcaga                                                                                                                  >  2:639119/1‑58 (MQ=255)                                                                                            agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGTAGACAGAcgccaattccgc                                                                                                  >  2:992635/1‑74 (MQ=255)                                                                                            agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCTGACAGAcgccaattccgc                                                                                                  <  1:992635/74‑1 (MQ=255)                                                                                            agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTc                                                 >  1:1400567/1‑123 (MQ=255)                                                                                            agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTc                                                 <  2:1400567/123‑1 (MQ=255)                                                                                            agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCAtaataa                                 >  2:1114036/1‑139 (MQ=255)                                                                                             gagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATaaa                                >  1:278495/1‑139 (MQ=255)                                                                                                  cTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCca                                                                                          >  1:1288711/1‑76 (MQ=255)                                                                                                  cTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCca                                                                                          <  2:1288711/76‑1 (MQ=255)                                                                                                  cTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGgt                           >  1:758646/1‑139 (MQ=255)                                                                                                   tttATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGgtg                          <  2:278495/139‑1 (MQ=255)                                                                                                      aTCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGgtg                          <  2:1048757/136‑1 (MQ=255)                                                                                                      aTCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGgtg                          >  1:1048757/1‑136 (MQ=255)                                                                                                             aaCTACATTTGTACCCGTGACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCca                                                                                          <  1:367123/65‑1 (MQ=255)                                                                                                             aaCTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCca                                                                                          >  2:367123/1‑65 (MQ=255)                                                                                                             aaCTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTAt                                                                    >  2:846901/1‑87 (MQ=255)                                                                                                             aaCTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTAt                                                                    <  1:846901/87‑1 (MQ=255)                                                                                                               cTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCca                                                                                          <  1:526640/63‑1 (MQ=255)                                                                                                               cTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCca                                                                                          >  2:526640/1‑63 (MQ=255)                                                                                                                     ttGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCg        <  2:633742/139‑1 (MQ=255)                                                                                                                       gTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCgtg      <  1:715285/139‑1 (MQ=255)                                                                                                                        tACCCGTAACGTTATGCTTCCGGCGGCAGTCAGAcgccaattccgcc                                                                                                 <  2:863588/47‑1 (MQ=38)                                                                                                                        tACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgcc                                                                                                 >  1:863588/1‑47 (MQ=255)                                                                                                                        tACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACggg                                                                           <  1:726174/69‑1 (MQ=255)                                                                                                                        tACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACggg                                                                           >  2:726174/1‑69 (MQ=255)                                                                                                                        tACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTc                                                 >  2:1467922/1‑95 (MQ=255)                                                                                                                        tACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTc                                                 <  1:1467922/95‑1 (MQ=255)                                                                                                                         aCCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCgtggt    <  1:1114036/139‑1 (MQ=255)                                                                                                                          cccGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCgg                                         >  2:1121643/1‑101 (MQ=255)                                                                                                                          cccGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCgg                                         <  1:1121643/101‑1 (MQ=255)                                                                                                                           ccGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGa  <  1:476864/139‑1 (MQ=255)                                                                                                                               aaCGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGt                                                                                      <  1:1446561/51‑1 (MQ=255)                                                                                                                               aaCGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGt                                                                                      >  2:1446561/1‑51 (MQ=255)                                                                                                                                   |                                                                                                                                   CCGATCAGGCCTACGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGA  >  NC_000913/4441408‑4441668

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