Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,735,899 T→C 7.3% Y174Y (TAT→TAC sodB → superoxide dismutase, Fe

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,735,8990TC 7.3% 163.5 / 3.5 55Y174Y (TAT→TACsodBsuperoxide dismutase, Fe
Reads supporting (aligned to +/- strand):  ref base T (25/26);  new base C (2/2);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.82e-01

GTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGG  >  NC_000913/1735756‑1736039
                                                                                                                                               |                                                                                                                                            
gTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTgg                                                                                                                                          <  3:270911/148‑1 (MQ=255)
  gAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGCCTCCGCTGCTGCCCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCCCGTCCTGGCTATCTGGAg                                                                                                                                        >  2:396036/1‑148 (MQ=255)
  gAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAg                                                                                                                                        >  1:328624/1‑148 (MQ=255)
   aaaaaCAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCAGGAGc                                                                                                                                       >  6:7562/1‑148 (MQ=255)
        cAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGc                                                                                                                                       >  8:204834/1‑143 (MQ=255)
        cAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGc                                                                                                                                       <  7:204834/143‑1 (MQ=255)
         aGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCt                                                                                                                                 <  7:352611/148‑1 (MQ=255)
          gCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTg                                                                                                                                <  5:7562/148‑1 (MQ=255)
                 cAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTg                                                                                                                         <  3:183369/148‑1 (MQ=255)
                   aaCTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGt                                                                                                                       >  8:13226/1‑148 (MQ=255)
                         cTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTg                                                                                                                 >  6:4080/1‑148 (MQ=255)
                             cGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGaa                                                                                                             >  4:155034/1‑148 (MQ=255)
                                      ctctAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGc                                                                                                     >  5:144152/1‑147 (MQ=255)
                                        ctAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGa                                                                                                   >  1:186228/1‑147 (MQ=255)
                                           aCGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGaaaa                                                                                                <  8:90043/147‑1 (MQ=255)
                                                    cTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGc                                                                                      <  8:231548/148‑1 (MQ=255)
                                                      ccGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCa                                                                                     >  7:153262/1‑147 (MQ=255)
                                                         cTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGa                                                                                                                     <  6:135375/112‑1 (MQ=255)
                                                         cTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGa                                                                                                                     >  5:135375/1‑112 (MQ=255)
                                                            accaccGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCataata                                                                                <  1:90123/146‑1 (MQ=255)
                                                            accaccGATGCGACTACGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCataata                                                                                <  5:4080/146‑1 (MQ=255)
                                                              caccGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTg                                                                            <  2:186228/148‑1 (MQ=255)
                                                               accGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGa                                                                           <  4:233355/148‑1 (MQ=255)
                                                                 cGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATg                                                                         >  7:387934/1‑148 (MQ=255)
                                                                    tGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCa                                                                                     <  4:111397/133‑1 (MQ=255)
                                                                    tGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCa                                                                                     >  3:111397/1‑133 (MQ=255)
                                                                            cGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATgcag                                                               <  8:387934/147‑1 (MQ=255)
                                                                              ctgctgACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATgcagcct                                                            >  1:310482/1‑146 (MQ=255)
                                                                                gctgACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTg                                                          >  3:411032/1‑148 (MQ=255)
                                                                                    aCCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCAGTTCTGGGCGCTGGGGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGcc                                                        <  8:159141/146‑1 (MQ=255)
                                                                                    aCCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGcc                                                        <  6:144152/146‑1 (MQ=255)
                                                                                                     acacGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGc                                                                                         <  8:168337/96‑1 (MQ=255)
                                                                                                     acacGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGc                                                                                         >  7:168337/1‑96 (MQ=255)
                                                                                                             attaCATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAAc                                                                              >  7:218840/1‑99 (MQ=255)
                                                                                                             attaCATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAAc                                                                              <  8:218840/99‑1 (MQ=255)
                                                                                                               taCATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGgcg                                                                                                                            <  5:68860/51‑1 (MQ=255)
                                                                                                               taCATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGgcg                                                                                                                            >  6:68860/1‑51 (MQ=255)
                                                                                                                  aTCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGat                        >  6:161590/1‑148 (MQ=255)
                                                                                                                   tCGACTATCGCAATGCACGTCCTGGCTACCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGat                        >  2:132604/1‑147 (MQ=255)
                                                                                                                   aagACTATCCCAATGCACGTCCTGGCTACCTGGAGCACTTCTGGGCGCTGGTGCACTGGGAATTCGTAGCGAAAAATCTCGCTACATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACTCTTATCATGACTACGCGGat                        >  8:181313/3‑147 (MQ=255)
                                                                                                                      aCTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCTTAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGatcatc                    >  8:123689/1‑148 (MQ=255)
                                                                                                                            gCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGt               <  3:155034/147‑1 (MQ=255)
                                                                                                                            gCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGt               <  5:161590/147‑1 (MQ=255)
                                                                                                                            gCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTa              <  4:411032/148‑1 (MQ=255)
                                                                                                                                 gCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCGGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCg         >  8:225738/1‑148 (MQ=255)
                                                                                                                                 gCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGat                        >  6:76958/1‑133 (MQ=255)
                                                                                                                                 gCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGat                        <  5:76958/133‑1 (MQ=255)
                                                                                                                                    cGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGAt      <  8:153262/148‑1 (MQ=255)
                                                                                                                                    cGTCCTGGCTACCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGAt      <  7:181313/148‑1 (MQ=255)
                                                                                                                                    cGTCCTGGCTACCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGAt      <  1:132604/148‑1 (MQ=255)
                                                                                                                                         tGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCaa                                                                     >  1:403269/1‑80 (MQ=255)
                                                                                                                                         tGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCaa                                                                     <  2:403269/80‑1 (MQ=255)
                                                                                                                                         tGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTCCGCGGATCATCGATGTAGGTCGGATAAgg  >  1:97703/1‑147 (MQ=255)
                                                                                                                                          ggCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATc                                      >  2:220462/1‑110 (MQ=255)
                                                                                                                                          ggCTATCTGGAGCACTTCTGGGCGCTGGGGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATc                                      <  1:220462/110‑1 (MQ=255)
                                                                                                                                               tCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCAtaataa                                                                               >  4:211768/1‑64 (MQ=255)
                                                                                                                                               tCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCAtaataa                                                                               <  3:211768/64‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                            
GTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGG  >  NC_000913/1735756‑1736039

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