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

CCGATCAGGCCTACGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑GTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGTG  >  NC_000913/4441408‑4441672                                                                                                                                   |                                                                                                                                        ccGATCAGGCCTACGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTc                                                                                                                                 <  2:167412/139‑1 (MQ=255)    aTCAGGCCTACGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCcgg                                                                                                                              <  2:284109/139‑1 (MQ=255)           cTACGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCAGcaga                                                                                                                       >  2:607361/1‑139 (MQ=255)                   tCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgcc                                                                                                               >  2:928772/1‑139 (MQ=255)                   tCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgcc                                                                                                               >  2:509595/1‑139 (MQ=255)                   tCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgcc                                                                                                               >  1:369358/1‑139 (MQ=255)                   tCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgcc                                                                                                               >  2:789401/1‑139 (MQ=255)                         cAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCAGCAGACAGAcgccaattcc                                                                                                         <  1:607361/139‑1 (MQ=255)                             atatTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgcca                                                                                                     >  2:859736/1‑139 (MQ=255)                                 tGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaattc                                                                                                 >  2:383089/1‑139 (MQ=255)                                       tGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑Gt                                                                                          <  2:173227/139‑1 (MQ=255)                                        ggACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑GTa                                                                                         >  2:705395/1‑139 (MQ=255)                                        ggACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑GTa                                                                                         >  1:579300/1‑139 (MQ=255)                                        ggACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑GTa                                                                                         >  2:804773/1‑139 (MQ=255)                                                ttGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑GTAACCATACg                                                                                 <  2:369358/139‑1 (MQ=255)                                                         tGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑GTAACCATACGGGTTTTTAt                                                                        <  1:509595/139‑1 (MQ=255)                                                           gTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgc                                                                                                       <  2:115576/107‑1 (MQ=255)                                                           gTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgc                                                                                                       >  1:115576/1‑107 (MQ=255)                                                            tAAGTAGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACTTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑GTAACCATACGGGTTTTTATGCa                                                                     <  2:628036/139‑1 (MQ=255)                                                               gTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑GTAACCATACGGGTTTTTATGCAGat                                                                  <  2:739489/139‑1 (MQ=255)                                                                    ttACGCCGCATCAAGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑GTAACCATACGGGTTTTTATGCAGATATTGc                                                             <  1:383089/139‑1 (MQ=255)                                                                     tACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑GTAACCATACGGGTTTTTATGCAGATATTGCt                                                            <  1:859736/139‑1 (MQ=255)                                                                            cATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑GTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTc                                                     >  1:719138/1‑139 (MQ=255)                                                                                 ggCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑GTAACCATACGGGTTTTTATGCa                                                                     >  2:762102/1‑118 (MQ=255)                                                                                 ggCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑GTAACCATACGGGTTTTTATGCa                                                                     <  1:762102/118‑1 (MQ=255)                                                                                    caggACAAAAAGGCCATTAACAACAAACTACATTTGTTACCCTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑GTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCgg                                             <  1:732278/136‑1 (MQ=255)                                                                                      aGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑Gt                                                                                          >  2:1006946/1‑92 (MQ=255)                                                                                      aGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑Gt                                                                                          <  1:1006946/92‑1 (MQ=255)                                                                                      aGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑GTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCa                                           <  2:579300/139‑1 (MQ=255)                                                                                      aGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑GTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCa                                           <  1:928772/139‑1 (MQ=255)                                                                                        aCAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑GTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCata                                         <  1:980318/139‑1 (MQ=255)                                                                                            agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑GTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCAtaataa                                     >  2:65625/1‑139 (MQ=255)                                                                                             gagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑GTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATaaa                                    <  1:804773/139‑1 (MQ=255)                                                                                                cACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGcaga                                                                                                                       >  1:276050/1‑54 (MQ=255)                                                                                                cACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGcaga                                                                                                                       <  2:276050/54‑1 (MQ=255)                                                                                                  cTTTATCAACAAACTACATTTGTACCCGTACTGTTATGCTTCCGGCGGAAGACAAAGGCCAATTCCGCC‑ATTTCACAGGTAACCATCCGGGTTTTTAGGAAGATATTGCTGGTGGCAATCTTCGCCTTAATAAAAccgc                               >  2:732278/1‑136 (MQ=255)                                                                                                  cTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAGTTCCACA‑GTAACCATACGGGTTTTTATGCAGATATTGCTggtggt                                                      >  1:687045/1‑116 (MQ=255)                                                                                                  cTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAGTTCCACA‑GTAACCATACGGGTTTTTATGCAGATATTGCTggaggt                                                      <  2:687045/116‑1 (MQ=255)                                                                                                     tATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑GTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCa                                           >  1:716437/1‑124 (MQ=255)                                                                                                     tATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑GTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCa                                           <  2:716437/124‑1 (MQ=255)                                                                                                             aaCTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaa                                                                                                    >  2:664488/1‑60 (MQ=255)                                                                                                             aaCTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaa                                                                                                    <  1:664488/60‑1 (MQ=255)                                                                                                             aaCTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACACCAATTCCGCCAATTCCACA‑GTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGc                    >  1:907025/1‑139 (MQ=255)                                                                                                                      tGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑GTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCgt           <  1:65625/139‑1 (MQ=255)                                                                                                                        tACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑GTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCa                                           <  1:326293/105‑1 (MQ=255)                                                                                                                        tACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑GTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCa                                           >  2:326293/1‑105 (MQ=255)                                                                                                                        tACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑GTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCgtgg         >  1:194984/1‑139 (MQ=255)                                                                                                                         aCCCGTAACGTTATGCTTCCGGCGGCAGACAGACACCAATTCCGCCAATTCCACA‑GTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCgtggt        <  2:907025/139‑1 (MQ=255)                                                                                                                             gTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑GTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCataata                                      >  2:593796/1‑105 (MQ=255)                                                                                                                             gTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑GTAACCATACGGGTTTTTATGCAGATATTGCTGATGGTCATCTTCGGCataata                                      <  1:593796/105‑1 (MQ=255)                                                                                                                               aaCGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑GTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGAtgtg  >  2:1001437/1‑139 (MQ=255)                                                                                                                                   |                                                                                                                                        CCGATCAGGCCTACGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACA‑GTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGTG  >  NC_000913/4441408‑4441672

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