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% 138.1 / NA 43*213* (TAG→TAAmsrAmethionine sulfoxide reductase A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/21);  total (22/21)

TATGCCGATCAGGCCTACGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGTGA  >  NC_000913/4441404‑4441673
                                                                                                                                      |                                                                                                                                        
tATGCCGATCAGGCCTACGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATg                                                                                                                                      >  2:83338/1‑139 (MQ=255)
            gCCTACGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGca                                                                                                                          >  1:241049/1‑139 (MQ=255)
              cTACGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGcaga                                                                                                                        <  1:640857/139‑1 (MQ=255)
                      tCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACTCGTAACGTTATGCTTCCGGCGGCAGACAGAcgcc                                                                                                                >  2:33628/1‑139 (MQ=255)
                      tCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCcgg                                                                                                                               <  2:719690/124‑1 (MQ=255)
                      tCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCcgg                                                                                                                               >  1:719690/1‑124 (MQ=255)
                      tCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgcc                                                                                                                >  1:24384/1‑139 (MQ=255)
                      tCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgcc                                                                                                                >  1:325902/1‑139 (MQ=255)
                           gCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattc                                                                                                           >  1:238408/1‑139 (MQ=255)
                              atatatTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgc                                                                                                        <  1:445448/139‑1 (MQ=255)
                                atatTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgcca                                                                                                      >  1:277214/1‑139 (MQ=255)
                                 tatTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaa                                                                                                     >  1:341795/1‑139 (MQ=255)
                                           ggACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaat                                                                                                    >  1:315212/1‑130 (MQ=255)
                                           ggACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaat                                                                                                    <  2:315212/130‑1 (MQ=255)
                                                             ggTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATg                                                                        <  2:277214/139‑1 (MQ=255)
                                                              gTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGc                                                                       <  2:325902/139‑1 (MQ=255)
                                                                  gTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGat                                                                   <  2:24384/139‑1 (MQ=255)
                                                                   tCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATAc                                                                                    >  2:753711/1‑122 (MQ=255)
                                                                   tCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATAc                                                                                    <  1:753711/122‑1 (MQ=255)
                                                                       ttACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGc                                                              <  2:502912/139‑1 (MQ=255)
                                                                       ttACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGc                                                              <  1:489405/139‑1 (MQ=255)
                                                                        tACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCt                                                             >  1:427650/1‑139 (MQ=255)
                                                                          cgccgcATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCtgg                                                           <  2:341795/139‑1 (MQ=255)
                                                                                aTCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCa                                                     >  1:10271/1‑139 (MQ=255)
                                                                                       aTAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCgg                                              <  2:427650/139‑1 (MQ=255)
                                                                                         aGACAAAGAGCACTTTATCAACAAACTACATTTGTACTCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCa                                            <  1:33628/139‑1 (MQ=255)
                                                                                               agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaattc                                                                                                  >  2:493809/1‑80 (MQ=255)
                                                                                               agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaattc                                                                                                  <  1:493809/80‑1 (MQ=255)
                                                                                               agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCa                                            <  1:693094/133‑1 (MQ=255)
                                                                                               agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCa                                            >  2:693094/1‑133 (MQ=255)
                                                                                               agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCAtaataa                                      >  2:249137/1‑139 (MQ=255)
                                                                                                           caacaaACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCg                          >  1:83084/1‑139 (MQ=255)
                                                                                                             acaaACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGc                                                                       >  2:12155/1‑92 (MQ=255)
                                                                                                             acaaACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGc                                                                       <  1:12155/92‑1 (MQ=255)
                                                                                                                aaCTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaattc                                                                                                  >  2:69127/1‑63 (MQ=255)
                                                                                                                aaCTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaattc                                                                                                  <  1:69127/63‑1 (MQ=255)
                                                                                                                        ttGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTcc              <  1:83338/139‑1 (MQ=255)
                                                                                                                        ttGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCg             <  2:238408/139‑1 (MQ=255)
                                                                                                                           tACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgc                                                                                                        <  2:243113/46‑1 (MQ=255)
                                                                                                                           tACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgc                                                                                                        >  1:243113/1‑46 (MQ=255)
                                                                                                                             cccGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCgtggtg        >  2:219266/1‑139 (MQ=255)
                                                                                                                             cccGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCgtggtg        >  1:331771/1‑139 (MQ=255)
                                                                                                                                   aCGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGTGa  <  2:241049/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                        
TATGCCGATCAGGCCTACGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGTGA  >  NC_000913/4441404‑4441673

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: