Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,442,007	T→C	E57G (GAG→GGG)	msrA ←	methionine sulfoxide reductase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,442,007	0	T	C	100.0%	173.4 / NA	49	E57G (GAG→GGG)	msrA	methionine sulfoxide reductase A
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/28);  total (21/28)

TCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGC  >  NC_000913/4441872‑4442140                                                                                                                                        |                                                                                                                                      tCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCca                                                                                                                                    <  2:1091212/139‑1 (MQ=255) tCGTAAACAATGCGTACCGCTTCGACATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCca                                                                                                                                    <  1:538361/139‑1 (MQ=255)      aaCAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACAccc                                                                                                                               <  1:390686/139‑1 (MQ=255)       aCAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACAcccc                                                                                                                              <  2:1094776/139‑1 (MQ=255)            gCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGaaa                                                                                                                         >  1:1217129/1‑139 (MQ=255)                aCCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTGGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCcac                                                                                                                                   >  1:662298/1‑125 (MQ=255)                aCCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTGGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCcac                                                                                                                                   <  2:662298/125‑1 (MQ=255)                   gCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCa                                                                                                                  >  2:297263/1‑139 (MQ=255)                                ccGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGc                                                                                                     <  2:335906/139‑1 (MQ=255)                                 cGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCa                                                                                                    <  1:1002780/139‑1 (MQ=255)                                   tATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAAt                                                                                                  >  2:37917/1‑139 (MQ=255)                                        ccGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCa                                                                                             <  2:1217129/139‑1 (MQ=255)                                            aGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTcc                                                                                         <  2:179828/139‑1 (MQ=255)                                            aGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTcc                                                                                         <  1:399775/139‑1 (MQ=255)                                             gCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCg                                                                                        <  1:297263/139‑1 (MQ=255)                                             gCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCg                                                                                        >  2:471866/1‑139 (MQ=255)                                               acacTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTc                                                                                      <  1:597896/139‑1 (MQ=255)                                               acacTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCAGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTc                                                                                      <  1:113671/139‑1 (MQ=255)                                                 acTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAg                                                                                    <  2:294187/139‑1 (MQ=255)                                                      ccGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACa                                                                               >  2:302232/1‑139 (MQ=255)                                                           aaGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGt                                                                          >  1:223418/1‑139 (MQ=255)                                                              tCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCAt                                                                       >  1:41239/1‑139 (MQ=255)                                                                      ggCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACAt                                                                              >  2:563437/1‑124 (MQ=255)                                                                      ggCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACAt                                                                              >  2:215710/1‑124 (MQ=255)                                                                      ggCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACAt                                                                              <  1:215710/124‑1 (MQ=255)                                                                      ggCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACAt                                                                              <  1:563437/124‑1 (MQ=255)                                                                        cGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACAcccc                                                                                                                              <  2:1066921/74‑1 (MQ=255)                                                                        cGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACAcccc                                                                                                                              >  1:1066921/1‑74 (MQ=255)                                                                        cGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACAccc                                                                                                                               >  2:1102518/1‑73 (MQ=255)                                                                        cGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACAccc                                                                                                                               <  1:1102518/73‑1 (MQ=255)                                                                         gTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGt                                                            <  1:367985/139‑1 (MQ=255)                                                                         gTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGt                                                            <  1:37917/139‑1 (MQ=255)                                                                         gTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGt                                                            <  2:223418/139‑1 (MQ=255)                                                                                       ataGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGc                                              <  1:302232/139‑1 (MQ=255)                                                                                             tGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATgcc                                                                                                     >  1:873333/1‑76 (MQ=255)                                                                                             tGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGc                                                                                                     <  2:873333/78‑1 (MQ=255)                                                                                               cggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTAc                                      >  1:962724/1‑139 (MQ=255)                                                                                               cggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTAc                                      >  1:249120/1‑139 (MQ=255)                                                                                                     tGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCACATGCAGCGTGGCTACGGGCAt                                >  1:459723/1‑139 (MQ=255)                                                                                                      gCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATc                               >  2:291785/1‑139 (MQ=255)                                                                                                         gTAAACGCCGGGTAACTGCCAGAACAGACGCCCcac                                                                                                                                   <  2:159919/36‑1 (MQ=255)                                                                                                         gTAAACGCCGGGTAACTGCCAGAACAGACGCCCcac                                                                                                                                   >  1:159919/1‑36 (MQ=255)                                                                                                             aCGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGgtgt                        <  2:962724/139‑1 (MQ=255)                                                                                                              cGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAgg                                                                                   <  1:992716/79‑1 (MQ=255)                                                                                                              cGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAgg                                                                                   >  2:992716/1‑79 (MQ=255)                                                                                                                       aCTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCcagg              <  2:41239/139‑1 (MQ=255)                                                                                                                       aCTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCcagg              <  1:471866/139‑1 (MQ=255)                                                                                                                                   gACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATcggc  >  1:603251/1‑139 (MQ=255)                                                                                                                                   gACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATcggc  >  2:941290/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                      TCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGC  >  NC_000913/4441872‑4442140

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: ‑