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%	224.6 / NA	70	*213* (TAG→TAA)	msrA	methionine sulfoxide reductase A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (35/35);  total (35/35)

TATGCCGATCAGGCCTACGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTG  >  NC_000913/4441404‑4441664                                                                                                                                       |                                                                                                                               tATGCCGATCAGGCCTACGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATg                                                                                                                            >  1:255450/1‑139 (MQ=255)                       tCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgcc                                                                                                      >  1:741370/1‑139 (MQ=255)                                     tGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCcgg                                                                                                                     >  1:325448/1‑110 (MQ=255)                                     tGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCcgg                                                                                                                     <  2:325448/110‑1 (MQ=255)                                          ttGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAg                                                                                   >  2:211236/1‑139 (MQ=255)                                           tGGACGATTTTGTAGGCCCGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGt                                                                                  >  1:1209914/1‑139 (MQ=255)                                            ggACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGcagaca                                                                                                            <  2:221482/112‑1 (MQ=255)                                            ggACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGcagaca                                                                                                            >  1:221482/1‑112 (MQ=255)                                            ggACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTa                                                                                 >  2:595634/1‑139 (MQ=255)                                              aCGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAAc                                                                               >  2:219536/1‑139 (MQ=255)                                                 aTTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCAt                                                                            >  2:912147/1‑139 (MQ=255)                                                            cTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATTAACAAACTACATTTGGACCCGTAACGTTATGCTTCCGGCGGcaga                                                                                                              <  1:964021/94‑1 (MQ=255)                                                            cTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGcaga                                                                                                              >  2:964021/1‑94 (MQ=255)                                                            cTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCca                                                                                      <  1:1423622/118‑1 (MQ=255)                                                            cTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCca                                                                                      >  2:1423622/1‑118 (MQ=255)                                                            cTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTa                                                                 <  1:595634/139‑1 (MQ=255)                                                             tGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTAt                                                                >  2:1595655/1‑139 (MQ=255)                                                             tGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTAt                                                                <  1:211236/139‑1 (MQ=255)                                                              ggTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATg                                                               <  2:1209914/139‑1 (MQ=255)                                                              ggTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATg                                                               <  2:255450/139‑1 (MQ=255)                                                               gTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgc                                                                                              >  2:1420199/1‑107 (MQ=255)                                                               gTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgc                                                                                              <  1:1420199/107‑1 (MQ=255)                                                               gTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGc                                                              >  2:971012/1‑139 (MQ=255)                                                                         tACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGt                                                                                  >  1:413833/1‑109 (MQ=255)                                                                         tACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGt                                                                                  <  2:413833/109‑1 (MQ=255)                                                                         tACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCt                                                    >  1:306680/1‑139 (MQ=255)                                                                           cgccgcATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCtgg                                                  <  1:219536/139‑1 (MQ=255)                                                                                cATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaattc                                                                                        <  1:281466/96‑1 (MQ=255)                                                                                cATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaattc                                                                                        >  2:281466/1‑96 (MQ=255)                                                                                cATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTc                                             <  2:741370/139‑1 (MQ=255)                                                                                     ggCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCtt                                        <  2:306680/139‑1 (MQ=255)                                                                                        aTAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAGCGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCca                                                                                      >  1:1181414/1‑90 (MQ=255)                                                                                        aTAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCca                                                                                      <  2:1181414/90‑1 (MQ=255)                                                                                          aGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCa                                   <  1:912147/139‑1 (MQ=255)                                                                                                agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCcg                                                                                                                      >  2:1167659/1‑50 (MQ=255)                                                                                                agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCcg                                                                                                                      <  1:1167659/50‑1 (MQ=255)                                                                                                agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattcc                                                                                                <  2:1588947/72‑1 (MQ=255)                                                                                                agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattcc                                                                                                >  1:1588947/1‑72 (MQ=255)                                                                                                agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGt                                                                                  >  1:888095/1‑86 (MQ=255)                                                                                                agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGt                                                                                  <  2:888095/86‑1 (MQ=255)                                                                                                agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTAt                                                                >  1:314179/1‑104 (MQ=255)                                                                                                agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTAt                                                                <  2:314179/104‑1 (MQ=255)                                                                                                  agCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATaaaa                           <  1:971012/139‑1 (MQ=255)                                                                                                       tttATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaatt                                                                                         <  2:1444008/72‑1 (MQ=255)                                                                                                       tttATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaatt                                                                                         >  1:1444008/1‑72 (MQ=255)                                                                                                        ttATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCa                                   >  1:95007/1‑125 (MQ=255)                                                                                                        ttATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCa                                   <  2:95007/125‑1 (MQ=255)                                                                                                                 aaCTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaattc                                                                                        >  1:1652282/1‑63 (MQ=255)                                                                                                                 aaCTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaattc                                                                                        <  2:1652282/63‑1 (MQ=255)                                                                                                                 aaCTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGc            >  1:691838/1‑139 (MQ=255)                                                                                                                   cTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCgg                                     >  1:676494/1‑112 (MQ=255)                                                                                                                   cTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCgg                                     <  2:676494/112‑1 (MQ=255)                                                                                                                         ttGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCg    <  1:1595655/139‑1 (MQ=255)                                                                                                                          tGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgc                                                                                              <  1:799401/48‑1 (MQ=255)                                                                                                                          tGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgc                                                                                              >  2:799401/1‑48 (MQ=255)                                                                                                                           gTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCgtg  <  2:691838/139‑1 (MQ=255)                                                                                                                              cccGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgc                                                                                              >  1:1166667/1‑44 (MQ=255)                                                                                                                              cccGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgc                                                                                              <  2:1166667/44‑1 (MQ=255)                                                                                                                                cGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaat                                                                                                   >  2:1433073/1‑37 (MQ=255)                                                                                                                                cGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaat                                                                                                   <  1:1433073/37‑1 (MQ=255)                                                                                                                                cGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATa                                                                           >  2:293805/1‑61 (MQ=255)                                                                                                                                cGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATa                                                                           <  1:293805/61‑1 (MQ=255)                                                                                                                                 gTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCtgg                                                  <  2:1205549/85‑1 (MQ=255)                                                                                                                                 gTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCtgg                                                  >  1:1205549/1‑85 (MQ=255)                                                                                                                                  tAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACg                                                                         <  2:174945/61‑1 (MQ=255)                                                                                                                                  tAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACg                                                                         >  1:174945/1‑61 (MQ=255)                                                                                                                                   aaCGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGc                                                              <  2:1518282/71‑1 (MQ=255)                                                                                                                                   aaCGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGc                                                              >  1:1518282/1‑71 (MQ=255)                                                                                                                                   aaCGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCgg                                     >  2:1100367/1‑96 (MQ=255)                                                                                                                                   aaCGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCgg                                     <  1:1100367/96‑1 (MQ=255)                                                                                                                                       |                                                                                                                               TATGCCGATCAGGCCTACGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTG  >  NC_000913/4441404‑4441664

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