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_0009134,441,5380CT100.0% 191.9 / NA 59*213* (TAG→TAAmsrAmethionine sulfoxide reductase A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/27);  total (32/27)

GCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGAT  >  NC_000913/4441431‑4441657
                                                                                                           |                                                                                                                       
gCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGcaga                                                                                                       >  1:77531/1‑126 (MQ=255)
gCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGcaga                                                                                                       <  2:77531/126‑1 (MQ=255)
gCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattc                                                                                          >  1:1541951/1‑139 (MQ=255)
   atatatTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgc                                                                                       <  2:1541951/139‑1 (MQ=255)
     atatTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgc                                                                                       >  2:227202/1‑137 (MQ=255)
     atatTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgc                                                                                       <  1:227202/137‑1 (MQ=255)
       atTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaat                                                                                   >  1:1365001/1‑139 (MQ=255)
           aaTTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCca                                                                               <  2:1365001/139‑1 (MQ=255)
              ttGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAg                                                                            >  1:299323/1‑139 (MQ=255)
                      ttttGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATa                                                                    <  2:299323/139‑1 (MQ=255)
                                cTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTa                                                          >  1:1415407/1‑139 (MQ=255)
                                  ggTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGcagaca                                                                                                     >  1:1446832/1‑94 (MQ=255)
                                  ggTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGcagaca                                                                                                     <  2:1446832/94‑1 (MQ=255)
                                   gTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGc                                                       >  1:904376/1‑139 (MQ=255)
                                    tAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCa                                                      <  2:1415407/139‑1 (MQ=255)
                                          gTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATg                                                        <  1:1146329/131‑1 (MQ=255)
                                          gTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATg                                                        >  2:1146329/1‑131 (MQ=255)
                                             tACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGAGAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCt                                             >  1:819355/1‑139 (MQ=255)
                                             tACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATa                                                                    >  2:385634/1‑116 (MQ=255)
                                             tACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATa                                                                    <  1:385634/116‑1 (MQ=255)
                                             tACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCa                                                      >  1:1356293/1‑130 (MQ=255)
                                             tACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCa                                                      <  2:1356293/130‑1 (MQ=255)
                                                                aCAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCata                          >  1:389239/1‑139 (MQ=255)
                                                                    agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGTTTCCGGCGGCAGACAGAcgccaattccgccaattc                                                                                 >  2:1196320/1‑80 (MQ=255)
                                                                    agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaattc                                                                                 <  1:1196320/80‑1 (MQ=255)
                                                                    agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCca                                                                               <  2:1461044/82‑1 (MQ=255)
                                                                    agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCca                                                                               >  1:1461044/1‑82 (MQ=255)
                                                                    agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATAc                                                                   >  1:1324901/1‑94 (MQ=255)
                                                                    agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATAc                                                                   <  2:1324901/94‑1 (MQ=255)
                                                                    agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTAt                                                         >  2:585413/1‑104 (MQ=255)
                                                                    agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTAt                                                         <  1:585413/104‑1 (MQ=255)
                                                                    agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCAtaataa                      >  2:1376852/1‑139 (MQ=255)
                                                                    agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCACATATTGCTGGTTGTCATCTTCGGCAtaataa                      >  1:1268023/1‑139 (MQ=255)
                                                                    agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGACAGACAGACGCCAATTCCGCCAATTCCATAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCAtaataa                      >  2:1128925/1‑139 (MQ=255)
                                                                          cTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattcc                                                                                         >  1:1430477/1‑66 (MQ=255)
                                                                          cTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattcc                                                                                         <  2:1430477/66‑1 (MQ=255)
                                                                                caacaaACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCg          <  2:904376/139‑1 (MQ=255)
                                                                                  acaaACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattcc                                                                                         >  1:778319/1‑58 (MQ=255)
                                                                                  acaaACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattcc                                                                                         <  2:778319/58‑1 (MQ=255)
                                                                                     aaCTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCca                                                                               <  1:975958/65‑1 (MQ=255)
                                                                                     aaCTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGc     >  1:312352/1‑139 (MQ=255)
                                                                                     aaCTACATTTGTACCCGTAACGTTATGCTACCGGCGGCAGACAGATGCCAATTCCGCCAATTCca                                                                               >  2:975958/1‑65 (MQ=255)
                                                                                        tACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGAt  <  2:1268023/139‑1 (MQ=255)
                                                                                           aTTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGt                                                                           <  2:1205952/63‑1 (MQ=255)
                                                                                           aTTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGt                                                                           >  1:1205952/1‑63 (MQ=255)
                                                                                                tACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCcaca                                                                             >  2:1300347/1‑56 (MQ=255)
                                                                                                tACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCcaca                                                                             <  1:1300347/56‑1 (MQ=255)
                                                                                                tACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCcaca                                                                             <  2:1378894/56‑1 (MQ=255)
                                                                                                tACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCcaca                                                                             >  1:1378894/1‑56 (MQ=255)
                                                                                                     gTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTAt                                                         <  1:668984/71‑1 (MQ=255)
                                                                                                     gTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTAt                                                         >  2:668984/1‑71 (MQ=255)
                                                                                                      tAACGTTATGCTTCCGGCGGCAGACAGAcgccaatt                                                                                           <  2:1364488/36‑1 (MQ=255)
                                                                                                      tAACGTTATGCTTCCGGCGGCAGACAGAcgccaatt                                                                                           >  1:1364488/1‑36 (MQ=255)
                                                                                                       aaCGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaattc                                                                                 <  2:560275/45‑1 (MQ=255)
                                                                                                       aaCGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaattc                                                                                 >  1:560275/1‑45 (MQ=255)
                                                                                                       aaCGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGt                                                                           <  2:717753/51‑1 (MQ=255)
                                                                                                       aaCGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGt                                                                           >  1:717753/1‑51 (MQ=255)
                                                                                                       aaCGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAAc                                                                        >  1:1438593/1‑54 (MQ=255)
                                                                                                       aaCGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAAc                                                                        <  2:1438593/54‑1 (MQ=255)
                                                                                                           |                                                                                                                       
GCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGAT  >  NC_000913/4441431‑4441657

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: