Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,735,897 T→C Y174H (TAT→CAT)  sodB → superoxide dismutase, Fe

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,735,8970TC100.0% 148.8 / NA 40Y174H (TAT→CAT) sodBsuperoxide dismutase, Fe
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/23);  total (17/23)

AAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTA  >  NC_000913/1735761‑1736027
                                                                                                                                        |                                                                                                                                  
aaaCAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCCAt                                                                                                                                  <  1:133908/139‑1 (MQ=255)
      cGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCACTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCCATCTGGAg                                                                                                                            >  1:863548/1‑139 (MQ=255)
       gATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCCATCTGGAGc                                                                                                                           <  2:1009034/139‑1 (MQ=255)
             aaaCTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCCATCTGGAGCACTTCt                                                                                                                     <  1:363250/139‑1 (MQ=255)
                          tttCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCCATCTGGAGCACTTCTggg                                                                                                                  >  2:902231/1‑129 (MQ=255)
                          tttCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCCATCTGGAGCACTTCTggg                                                                                                                  <  1:902231/129‑1 (MQ=255)
                           ttCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAAc                                                                                                       <  1:924034/139‑1 (MQ=255)
                                           ggTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGa                                                                                       >  1:230472/1‑139 (MQ=255)
                                                    cTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCg                                                                              <  2:393861/139‑1 (MQ=255)
                                                    cTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCg                                                                              <  1:62416/139‑1 (MQ=255)
                                                     ttgctaccGGTCCGACCCCGGTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTACTGGCCATCTGGAGCATTTCTGGGCTGTGGTGAACTGGAAATTCGTAACGAAAACTCTCGc                                                                             >  2:1049451/6‑139 (MQ=255)
                                                     tGACCCCCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATTCACGTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGc                                                                             <  1:688707/139‑1 (MQ=255)
                                                     tGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGc                                                                             <  2:508439/139‑1 (MQ=255)
                                                       accaccGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTg                                                                           <  1:771759/139‑1 (MQ=255)
                                                       accaccGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTg                                                                           <  1:501832/139‑1 (MQ=255)
                                                                       cGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGc                                                           <  2:899901/139‑1 (MQ=255)
                                                                         ctgctgACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGc                                                                             <  2:50621/119‑1 (MQ=255)
                                                                         ctgctgACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGc                                                                             >  1:50621/1‑119 (MQ=255)
                                                                                 cGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCg                                                                              <  1:879933/110‑1 (MQ=255)
                                                                                 cGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCg                                                                              >  2:879933/1‑110 (MQ=255)
                                                                                                       tattaCATCGACTATCGCAATGCACGTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTg                                                                           <  1:906043/91‑1 (MQ=255)
                                                                                                       tattaCATCGACTATCGCAATGCACGTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTg                                                                           >  2:906043/1‑91 (MQ=255)
                                                                                                       tattaCATCGACTATCGCAATGCACGTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTAt                           <  1:965790/139‑1 (MQ=255)
                                                                                                       tattaCATCGACTATCGCAATGCACGTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTAt                           <  2:83755/139‑1 (MQ=255)
                                                                                                         ttaCATCGACTATCGCAATGCACGTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCa                         >  1:632817/1‑139 (MQ=255)
                                                                                                              tCGACTATCGCAATGCACGTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGa                                                                                                         <  1:795575/54‑1 (MQ=255)
                                                                                                              tCGACTATCGCAATGCACGTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGa                                                                                                         >  2:795575/1‑54 (MQ=255)
                                                                                                              tCGACTATCGCAATGCACGTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCt                    >  2:493215/1‑139 (MQ=255)
                                                                                                                gACTATCGCAATGCACGTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTa                   >  2:179066/1‑138 (MQ=255)
                                                                                                                gACTATCGCAATGCACGTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTa                   <  1:179066/138‑1 (MQ=255)
                                                                                                                gACTATCGCAATGCACGTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTAc                  >  2:20565/1‑139 (MQ=255)
                                                                                                                 aCTATCGCAATGCACGTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTAca                 >  1:329162/1‑138 (MQ=255)
                                                                                                                      cGCAATGCACGTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTg                                                                           >  1:716844/1‑76 (MQ=255)
                                                                                                                      cGCAATGCACGTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTg                                                                           <  2:716844/76‑1 (MQ=255)
                                                                                                                        cAATGCACGTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGatca          <  2:230472/139‑1 (MQ=255)
                                                                                                                           tGCACGTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGc                                                                             >  2:599554/1‑69 (MQ=255)
                                                                                                                           tGCACGTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGc                                                                             <  1:599554/69‑1 (MQ=255)
                                                                                                                                gTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAAtctc                                                                               <  2:674163/62‑1 (MQ=255)
                                                                                                                                gTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAAtctc                                                                               >  1:674163/1‑62 (MQ=255)
                                                                                                                                gTCCTGGCCATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTa  >  2:546416/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                  
AAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTA  >  NC_000913/1735761‑1736027

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