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_0009134,442,0070TC100.0% 125.8 / NA 36E57G (GAG→GGG) msrAmethionine sulfoxide reductase A
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/18);  total (18/18)

TGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCG  >  NC_000913/4441882‑4442141
                                                                                                                             |                                                                                                                                      
tGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCca                                                                                                                                     <  2:904704/129‑1 (MQ=255)
tGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCca                                                                                                                                     >  1:904704/1‑129 (MQ=255)
              gCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCaa                                                                                                             >  2:571448/1‑139 (MQ=255)
                            aCCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTcc                                                                                               <  1:571448/139‑1 (MQ=255)
                               ggAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCAtt                                                                                            <  1:1068100/139‑1 (MQ=255)
                                  gCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACAcccc                                                                                                                               >  1:189399/1‑101 (MQ=255)
                                  gCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACAcccc                                                                                                                               <  2:189399/101‑1 (MQ=255)
                                        ttCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCCCCATTCCGTCAGGt                                                                                   >  1:434694/1‑139 (MQ=255)
                                           ccGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCca                                                                                                                                     <  2:940547/86‑1 (MQ=255)
                                           ccGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCca                                                                                                                                     >  1:940547/1‑86 (MQ=255)
                                             gATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACAccc                                                                                                                                >  2:873249/1‑89 (MQ=255)
                                             gATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACAccc                                                                                                                                <  1:873249/89‑1 (MQ=255)
                                                  gTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGt                                                                                   <  1:15834/129‑1 (MQ=255)
                                                  gTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGt                                                                                   >  2:15834/1‑129 (MQ=255)
                                                                      gccgGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCAtt                                                                                            <  2:333577/100‑1 (MQ=255)
                                                                      gccgGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCAtt                                                                                            >  1:333577/1‑100 (MQ=255)
                                                                           tataGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAtacaac                                                                                                                      >  2:737434/1‑69 (MQ=255)
                                                                           tataGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAaacaac                                                                                                                      <  1:737434/69‑1 (MQ=255)
                                                                                  tGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACAccc                                                                                                                                >  1:312158/1‑52 (MQ=255)
                                                                                  tGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACAccc                                                                                                                                <  2:312158/52‑1 (MQ=255)
                                                                                    cggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCaaaaa                                                                                                          <  2:185989/72‑1 (MQ=255)
                                                                                    cggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCaaaaa                                                                                                          >  1:185989/1‑72 (MQ=255)
                                                                                     ggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTATgg                                     >  2:153182/1‑139 (MQ=255)
                                                                                            ctctAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACTCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATg                                                   >  2:658550/4‑119 (MQ=255)
                                                                                            cTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATg                                                   <  1:658550/119‑1 (MQ=255)
                                                                                                 aaCGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTATGGGCATCGGGgtgt                         <  1:153182/139‑1 (MQ=255)
                                                                                                      cGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAg                                                                    <  2:277444/92‑1 (MQ=255)
                                                                                                      cGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAg                                                                    >  1:277444/1‑92 (MQ=255)
                                                                                                         gTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCa                                                                                              >  1:494862/1‑63 (MQ=255)
                                                                                                         gTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGgtgt                         >  2:585543/1‑132 (MQ=255)
                                                                                                         gTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGgtgt                         <  1:585543/132‑1 (MQ=255)
                                                                                                         gTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAAGCGCCa                                                                                              <  2:494862/63‑1 (MQ=255)
                                                                                                                        gACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCa                                                                                              <  1:863375/48‑1 (MQ=255)
                                                                                                                        gACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCa                                                                                              >  2:863375/1‑48 (MQ=255)
                                                                                                                         aCGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATcggcg  >  1:19523/1‑139 (MQ=255)
                                                                                                                          cGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATg                                                   <  2:968476/89‑1 (MQ=255)
                                                                                                                             |                                                                                                                                      
TGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCG  >  NC_000913/4441882‑4442141

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: