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%	208.3 / NA	58	E57G (GAG→GGG)	msrA	methionine sulfoxide reductase A
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/31);  total (27/31)

TCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGG  >  NC_000913/4441872‑4442128                                                                                                                                        |                                                                                                                          tCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCca                                                                                                                        <  2:332292/139‑1 (MQ=255)      aaCAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACAccc                                                                                                                   >  2:311570/1‑139 (MQ=255)       aCAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACAcccc                                                                                                                  <  1:646617/139‑1 (MQ=255)       aCAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACAcccc                                                                                                                  <  1:639874/139‑1 (MQ=255)          aTGCGTACCGCTTCGGCTTGACCTGTATCACCGGAGCACACTTCCCGAGAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCATGCGGTGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGa                                                                                                               >  2:1278611/1‑139 (MQ=255)           tGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCAGACTTCACGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAAAAGACGCCCCACACCCCAGaa                                                                                                              >  1:957467/1‑139 (MQ=255)           tGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGaa                                                                                                              >  1:1179012/1‑139 (MQ=255)                aCCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAaacaac                                                                                                         >  2:902541/1‑139 (MQ=255)                            tGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCaaaaa                                                                                             <  1:496846/139‑1 (MQ=255)                                ccGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGc                                                                                         <  2:735169/139‑1 (MQ=255)                                 cGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCa                                                                                        <  1:902541/139‑1 (MQ=255)                                        ccGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCa                                                                                 >  2:510588/1‑139 (MQ=255)                                         cGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGTGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCAt                                                                                >  1:672616/1‑139 (MQ=255)                                         cGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACAcccc                                                                                                                  >  1:1101473/1‑105 (MQ=255)                                         cGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACAcccc                                                                                                                  <  2:1101473/105‑1 (MQ=255)                                           gAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCcac                                                                                                                       <  1:1271080/98‑1 (MQ=255)                                           gAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCcac                                                                                                                       >  2:1271080/1‑98 (MQ=255)                                            aGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTcc                                                                             <  1:311570/139‑1 (MQ=255)                                            aGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTcc                                                                             <  1:5574/139‑1 (MQ=255)                                            aGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTcc                                                                             <  2:550442/139‑1 (MQ=255)                                              cacaCTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGt                                                                           >  2:221270/1‑139 (MQ=255)                                              cacaCTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGt                                                                           <  1:1096947/139‑1 (MQ=255)                                                 acTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAg                                                                        >  1:420309/1‑139 (MQ=255)                                                             gTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCa                                                            >  2:945514/1‑139 (MQ=255)                                                                             aGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGAc                                            <  2:1179012/139‑1 (MQ=255)                                                                                 gccgGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCa                                        <  1:945514/139‑1 (MQ=255)                                                                                 gccgGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCa                                        <  1:990538/139‑1 (MQ=255)                                                                                  ccgGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCAt                                       >  1:670812/1‑139 (MQ=255)                                                                                        taGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCg                                 <  2:670812/139‑1 (MQ=255)                                                                                         aGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGc                                                                                         >  2:688966/1‑82 (MQ=255)                                                                                         aGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGc                                                                                         <  1:688966/82‑1 (MQ=255)                                                                                               cggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACAcc                                                                                                                    >  1:814720/1‑49 (MQ=255)                                                                                               cggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACAcc                                                                                                                    <  2:814720/49‑1 (MQ=255)                                                                                               cggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCAtt                                                                               >  1:825777/1‑86 (MQ=255)                                                                                               cggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCAtt                                                                               <  2:825777/86‑1 (MQ=255)                                                                                               cggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCtcc                          >  2:760010/1‑137 (MQ=255)                                                                                               cggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTAc                          >  1:1092528/1‑139 (MQ=255)                                                                                               cggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTAc                          >  2:677535/1‑139 (MQ=255)                                                                                               cggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTAc                          >  2:39022/1‑139 (MQ=255)                                                                                                ggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCca                                                                                                                        >  1:1155686/1‑44 (MQ=255)                                                                                                ggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCca                                                                                                                        <  2:1155686/44‑1 (MQ=255)                                                                                                     tGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACATCCCATCGCAAAAATGGCAAt                                                                                      <  2:761334/73‑1 (MQ=255)                                                                                                     tGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAAt                                                                                      >  1:761334/1‑73 (MQ=255)                                                                                                     tGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCAt                    <  1:510588/139‑1 (MQ=255)                                                                                                         gTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCggg                <  2:1242507/139‑1 (MQ=255)                                                                                                             aCGCCGGGTAACTGCCAGAACTGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGgtgt            <  1:760010/139‑1 (MQ=255)                                                                                                             aCGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGgtgt            <  2:672616/139‑1 (MQ=255)                                                                                                               gCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTcc                                                                             <  1:926017/72‑1 (MQ=255)                                                                                                               gCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTcc                                                                             >  2:926017/1‑72 (MQ=255)                                                                                                                 cGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCAtt                                                                               <  2:1283770/68‑1 (MQ=255)                                                                                                                 cGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCAtt                                                                               >  1:1283770/1‑68 (MQ=255)                                                                                                                 cGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTAc                          <  1:279421/121‑1 (MQ=255)                                                                                                                 cGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTAc                          >  2:279421/1‑121 (MQ=255)                                                                                                                       aCTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCcagg  <  1:677535/139‑1 (MQ=255)                                                                                                                              gAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCa                                                            >  2:670604/1‑74 (MQ=255)                                                                                                                              gAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCa                                                            <  1:670604/74‑1 (MQ=255)                                                                                                                                   gACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCa                                                                                        <  1:1287008/41‑1 (MQ=255)                                                                                                                                   gACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCa                                                                                        >  2:1287008/1‑41 (MQ=255)                                                                                                                                        |                                                                                                                          TCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGG  >  NC_000913/4441872‑4442128

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