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

ATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCG  >  NC_000913/4441881‑4442141
                                                                                                                              |                                                                                                                                      
aTGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGa                                                                                                                            >  1:8990/1‑139 (MQ=255)
 tGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGaa                                                                                                                           >  2:473619/1‑139 (MQ=255)
  gCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGaaa                                                                                                                          >  1:543137/1‑139 (MQ=255)
         gCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCcac                                                                                                                                    <  1:555018/122‑1 (MQ=255)
         gCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCcac                                                                                                                                    >  2:555018/1‑122 (MQ=255)
           ttCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATc                                                                                                                 >  1:242703/1‑139 (MQ=255)
             cGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGc                                                                                                               >  2:626357/1‑139 (MQ=255)
                cATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCaaa                                                                                                            <  2:242703/139‑1 (MQ=255)
                              ccGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCa                                                                                              >  2:685859/1‑139 (MQ=255)
                                ggAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCAtt                                                                                            <  1:473619/139‑1 (MQ=255)
                                  aGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTcc                                                                                          <  2:309709/139‑1 (MQ=255)
                                         ttCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGt                                                                                   <  2:91165/139‑1 (MQ=255)
                                                   gTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCa                                                                         >  1:343946/1‑139 (MQ=255)
                                                   gTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCa                                                                         >  1:1040925/1‑139 (MQ=255)
                                                    tCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAACCGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAAAAACCACCCTCCGCAACAATGCCAATCTCCATTCCGTCAGGTACATTGGTCAt                                                                        <  2:172023/139‑1 (MQ=255)
                                                        aTTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAg                                                                    >  2:880017/1‑139 (MQ=255)
                                                        aTTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAg                                                                    <  2:1040925/139‑1 (MQ=255)
                                                               gTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGt                                                             <  1:685859/139‑1 (MQ=255)
                                                                 ataGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTg                                                           <  2:343946/139‑1 (MQ=255)
                                                                  taGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAg                                                                                     >  2:423902/1‑112 (MQ=255)
                                                                  taGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAg                                                                                     <  1:423902/112‑1 (MQ=255)
                                                                  taGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGa                                                          >  2:1402366/1‑139 (MQ=255)
                                                                         cgGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACCGCCAGAACAGACGCCCCACACCCCAGCAACACCCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCTTg                                                   <  1:1268254/139‑1 (MQ=255)
                                                                             ataGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGATCGCATGCAGc                                               <  1:469585/139‑1 (MQ=255)
                                                                             ataGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGc                                               <  1:880017/139‑1 (MQ=255)
                                                                                     cggcggTGCTGTAAACGCCGGTTAACTGCCAGAACAGACGCCCCACAcc                                                                                                                                 <  2:916305/49‑1 (MQ=38)
                                                                                     cggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCcaca                                                                                                                                   >  1:1245339/1‑47 (MQ=255)
                                                                                     cggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCcaca                                                                                                                                   <  2:1245339/47‑1 (MQ=255)
                                                                                     cggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCca                                                                                                                                     <  2:626565/45‑1 (MQ=255)
                                                                                     cggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCca                                                                                                                                     <  2:1435503/45‑1 (MQ=255)
                                                                                     cggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCca                                                                                                                                     >  1:626565/1‑45 (MQ=255)
                                                                                     cggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCca                                                                                                                                     >  1:1435503/1‑45 (MQ=255)
                                                                                     cggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACAcc                                                                                                                                 >  1:916305/1‑49 (MQ=255)
                                                                                     cggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTAc                                       >  1:975174/1‑139 (MQ=255)
                                                                                      ggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACAcc                                                                                                                                 >  1:472250/1‑48 (MQ=255)
                                                                                      ggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACAcc                                                                                                                                 <  2:472250/48‑1 (MQ=255)
                                                                                      ggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTc                                                                                       <  2:422969/90‑1 (MQ=255)
                                                                                      ggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTc                                                                                       >  1:422969/1‑90 (MQ=255)
                                                                                      ggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGt                                                             >  1:624744/1‑116 (MQ=255)
                                                                                      ggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGt                                                             <  2:624744/116‑1 (MQ=255)
                                                                                      ggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCg                                              >  1:557236/1‑131 (MQ=255)
                                                                                      ggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCg                                              <  2:557236/131‑1 (MQ=255)
                                                                                          gTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTc                                                                                       <  2:23783/86‑1 (MQ=255)
                                                                                          gTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTc                                                                                       >  1:23783/1‑86 (MQ=255)
                                                                                          gTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCa                                  >  1:575861/1‑139 (MQ=255)
                                                                                           tGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCAt                                 >  1:1222419/1‑139 (MQ=255)
                                                                                           tGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCAt                                 <  2:575861/139‑1 (MQ=255)
                                                                                              tGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCgg                              <  1:1125432/139‑1 (MQ=255)
                                                                                                   aCGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGgtgt                         <  2:543137/139‑1 (MQ=255)
                                                                                                    cGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCgg                              <  2:837025/133‑1 (MQ=255)
                                                                                                    cGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCgg                              >  1:837025/1‑133 (MQ=255)
                                                                                                     gCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATg                                                   <  2:602486/111‑1 (MQ=255)
                                                                                                     gCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATg                                                   >  1:602486/1‑111 (MQ=255)
                                                                                                     gCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTa                       >  1:178840/1‑139 (MQ=255)
                                                                                                       cGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACAt                                                                               <  2:268180/81‑1 (MQ=255)
                                                                                                       cGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACAt                                                                               >  1:268180/1‑81 (MQ=255)
                                                                                                       cGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACg                     >  1:1055975/1‑139 (MQ=255)
                                                                                                                   aGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGc         >  1:405486/1‑139 (MQ=255)
                                                                                                                     aaCAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCAt       <  2:8990/139‑1 (MQ=255)
                                                                                                                        aGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATcgg    <  2:1222419/139‑1 (MQ=255)
                                                                                                                          aCGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATcggcg  <  1:626357/139‑1 (MQ=255)
                                                                                                                              |                                                                                                                                      
ATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCG  >  NC_000913/4441881‑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: