Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,713,568 T→C W225R (TGG→CGG)  srmB → ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,713,5680TC100.0% 137.8 / NA 40W225R (TGG→CGG) srmBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/22);  total (18/22)

AGACCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCT  >  NC_000913/2713443‑2713686
                                                                                                                             |                                                                                                                      
aGACCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCgcgc                                                                                                           <  2:1147401/139‑1 (MQ=255)
  aCCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCg                                                                                                         >  2:124017/1‑139 (MQ=255)
  aCCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCg                                                                                                         >  1:590024/1‑139 (MQ=255)
  aCCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTAGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCg                                                                                                         >  2:242613/1‑139 (MQ=255)
            tttCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGa                                                                                               >  2:687355/1‑139 (MQ=255)
                     cGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCTCCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAAc                                                                                      >  2:841962/1‑139 (MQ=255)
                     cGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAAc                                                                                      <  1:687355/139‑1 (MQ=255)
                          gAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGt                                                                                 <  2:258902/139‑1 (MQ=255)
                                  tgcgatTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTg                                                                         <  1:124017/139‑1 (MQ=255)
                                  tgcgatCCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCAGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGATGGTg                                                                         >  1:461592/1‑139 (MQ=255)
                                    cgatTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCa                                                                       <  1:328585/139‑1 (MQ=255)
                                    cgatTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCa                                                                       <  1:1225525/139‑1 (MQ=255)
                                    cgatTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCa                                                                       <  1:841962/139‑1 (MQ=255)
                                    cgacTCAGGACTTTGCCGCGCGTCTGCTGTAAGATCCGGTGGAAGTTTCTCCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCa                                                                       <  1:755335/139‑1 (MQ=255)
                                                      aGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGa                                                     >  1:436461/1‑139 (MQ=255)
                                                         gTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCgcgc                                                                                                           >  2:754627/1‑82 (MQ=255)
                                                         gTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCgcgc                                                                                                           <  1:754627/82‑1 (MQ=255)
                                                            tgctgGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGAc                                               >  1:503358/1‑139 (MQ=255)
                                                               tgGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCg                                            <  2:590024/139‑1 (MQ=255)
                                                                      tCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAAtt                                     <  2:461592/139‑1 (MQ=255)
                                                                                  ttCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGt                         <  2:503358/139‑1 (MQ=255)
                                                                                        cAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAagcgt                   <  1:1165862/139‑1 (MQ=255)
                                                                                         aaTCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAagcgtg                  <  2:973301/139‑1 (MQ=255)
                                                                                           tCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCg                                                                                                         <  2:1173823/50‑1 (MQ=255)
                                                                                           tCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCg                                                                                                         >  1:1173823/1‑50 (MQ=255)
                                                                                           tCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAagcgtgag                >  1:768814/1‑139 (MQ=255)
                                                                                                  aCCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGc         <  2:768814/139‑1 (MQ=255)
                                                                                                      gTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCg                                                 >  1:58656/1‑95 (MQ=255)
                                                                                                      gTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCg                                                 <  2:58656/95‑1 (MQ=255)
                                                                                                       tGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAagcgtga                 <  2:1359531/126‑1 (MQ=255)
                                                                                                       tGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAagcgtga                 >  1:1359531/1‑126 (MQ=255)
                                                                                                         aGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCga                                                                                                        <  2:918671/37‑1 (MQ=255)
                                                                                                         aGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCga                                                                                                        >  1:918671/1‑37 (MQ=255)
                                                                                                         aGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCt  >  2:1351092/1‑139 (MQ=255)
                                                                                                           cgcAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATa                                                                                          >  1:1007715/1‑49 (MQ=255)
                                                                                                           cgcAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATa                                                                                          <  2:1007715/49‑1 (MQ=255)
                                                                                                           cgcAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACcg                                                                                    <  1:585826/55‑1 (MQ=255)
                                                                                                           cgcAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACcg                                                                                    >  2:585826/1‑55 (MQ=255)
                                                                                                             cAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGAcc                                              <  1:176651/91‑1 (MQ=255)
                                                                                                             cAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGAcc                                              >  2:176651/1‑91 (MQ=255)
                                                                                                                             |                                                                                                                      
AGACCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCT  >  NC_000913/2713443‑2713686

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