Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,439,845 T→A 10.1% intergenic (+52/+27) ytfK → / ← ytfL DUF1107 family protein/UPF0053 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,439,8450TA10.1% 189.4 / 12.2 60intergenic (+52/+27)ytfK/ytfLDUF1107 family protein/UPF0053 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (26/27);  new base A (3/3);  total (30/30)
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.06e-01

ATTCTTATCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCGATATCCACCA  >  NC_000913/4439701‑4439990
                                                                                                                                                |                                                                                                                                                 
aTTCTTATCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAg                                                                                                                                               <  1:125919/149‑1 (MQ=255)
aTTCTTATCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAg                                                                                                                                               <  2:775315/149‑1 (MQ=255)
         ccAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCa                                                                                                                                                >  4:56134/1‑139 (MQ=255)
         ccAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCa                                                                                                                                                <  3:56134/139‑1 (MQ=255)
          cacaaGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATCAGACGGCTCCTGATTCAGGAGCCGTTGa                                                                                                                                     >  2:276363/4‑149 (MQ=255)
          cAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGa                                                                                                                                     <  1:276363/149‑1 (MQ=255)
          cAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAATCAGGAGcc                                                                                                                                          <  2:933900/144‑1 (MQ=255)
          cAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAATCAGGAGcc                                                                                                                                          >  1:933900/1‑144 (MQ=255)
                     gATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGa                                                                                                                                             <  3:411658/130‑1 (MQ=255)
                     gATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGa                                                                                                                                             >  4:411658/1‑130 (MQ=255)
                     gATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTgg                                                                                                                            >  4:64468/1‑147 (MQ=255)
                     gATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTgg                                                                                                                            <  3:64468/147‑1 (MQ=255)
                       tAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGcc                                                                                                                                          <  4:1190522/131‑1 (MQ=255)
                       tAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGcc                                                                                                                                          >  3:1190522/1‑131 (MQ=255)
                             gCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGa                                                                                                                                             <  1:419687/122‑1 (MQ=255)
                             gCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGa                                                                                                                                             >  2:419687/1‑122 (MQ=255)
                                    tCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGAt                                                                                                                                    >  1:1275777/1‑124 (MQ=255)
                                    tCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGAt                                                                                                                                    <  2:1275777/124‑1 (MQ=255)
                                      aGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGa                                                                                                                                     <  3:570805/121‑1 (MQ=255)
                                      aGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGa                                                                                                                                     >  4:570805/1‑121 (MQ=255)
                                                   gTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatcttt                                                                                              <  3:584940/147‑1 (MQ=255)
                                                   gTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatcttt                                                                                              >  4:584940/1‑147 (MQ=255)
                                                    ttAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatcttt                                                                                              <  4:944920/146‑1 (MQ=255)
                                                    ttAACCGTCAGATTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatcttt                                                                                              >  3:944920/1‑146 (MQ=255)
                                                                ttATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAAACGGCTCCTGAATCAGGAGCCGtctt                                                                                                                                     <  2:14467/95‑4 (MQ=255)
                                                                ttATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAAACGGCTCCTGAATCAGGAGCCGtctt                                                                                                                                     >  1:14467/1‑92 (MQ=39)
                                                                  aTGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGtt                                                                                                                                       <  3:931818/91‑1 (MQ=255)
                                                                  aTGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGtt                                                                                                                                       >  4:931818/1‑91 (MQ=255)
                                                                                    aTGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAtt                                                                                                                                                  <  2:1412591/62‑1 (MQ=255)
                                                                                    aTGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAtt                                                                                                                                                  >  1:1412591/1‑62 (MQ=255)
                                                                                      ggCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGa                                                                                                                                             >  2:1099483/1‑65 (MQ=255)
                                                                                      ggCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGa                                                                                                                                             <  1:1099483/65‑1 (MQ=255)
                                                                                              cAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGc                                                             <  1:1035353/137‑1 (MQ=255)
                                                                                              cAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGc                                                             >  1:1391574/1‑137 (MQ=255)
                                                                                              cAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGc                                                             >  2:1035353/1‑137 (MQ=255)
                                                                                              cAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGc                                                             <  2:1391574/137‑1 (MQ=255)
                                                                                              cAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGAt                                                 >  3:1071297/1‑149 (MQ=255)
                                                                                                aagtgctatgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTgggg                                                                                                                          >  3:1083020/1‑74 (MQ=255)
                                                                                                aagtgctatgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCATTGATGTTTCTgggg                                                                                                                          <  4:1083020/74‑1 (MQ=255)
                                                                                                       atgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatc                                                                                                 >  4:442774/1‑92 (MQ=255)
                                                                                                       atgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatc                                                                                                 <  3:442774/92‑1 (MQ=255)
                                                                                                               aaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGt                                                     <  2:48806/128‑1 (MQ=255)
                                                                                                               aaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGt                                                     >  1:48806/1‑128 (MQ=255)
                                                                                                                            agtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGAAGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTg                                                         >  4:723804/1‑111 (MQ=255)
                                                                                                                            agtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGAAGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTg                                                         <  3:723804/111‑1 (MQ=255)
                                                                                                                             gtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCgg                                                                 <  4:672854/102‑1 (MQ=255)
                                                                                                                             gtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGACGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCgg                                                                 >  3:672854/1‑102 (MQ=255)
                                                                                                                                      cGGCTCCTGAATCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatc                                                                                                 <  1:612599/61‑1 (MQ=255)
                                                                                                                                      cGGCTCCTGAATCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatc                                                                                                 >  2:612599/1‑61 (MQ=255)
                                                                                                                                          tCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTa                   >  1:97006/1‑135 (MQ=255)
                                                                                                                                          tCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTa                   <  2:97006/135‑1 (MQ=255)
                                                                                                                                          tCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCGATATcca     >  1:354672/1‑149 (MQ=255)
                                                                                                                                           ccTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCGATATccac    <  4:1071297/149‑1 (MQ=255)
                                                                                                                                            cTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGtt                <  1:250082/136‑1 (MQ=255)
                                                                                                                                            cTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGtt                >  2:250082/1‑136 (MQ=255)
                                                                                                                                             tGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGaaa                                                                       <  1:347395/80‑1 (MQ=255)
                                                                                                                                             tGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCGATATccacca  >  3:103777/1‑149 (MQ=255)
                                                                                                                                             tGAGTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGaaa                                                                       >  2:347395/1‑80 (MQ=255)
                                                                                                                                              gATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACcagca                                   >  4:333805/1‑115 (MQ=255)
                                                                                                                                              gATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACcagca                                   <  3:333805/115‑1 (MQ=255)
                                                                                                                                               aTTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCt                                                           <  1:983989/90‑1 (MQ=255)
                                                                                                                                               aTTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCt                                                           >  2:983989/1‑90 (MQ=255)
                                                                                                                                                ttCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCagcagc                                 >  1:799742/1‑115 (MQ=255)
                                                                                                                                                ttCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCAACagcagc                                 <  2:799742/115‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
ATTCTTATCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCGATATCCACCA  >  NC_000913/4439701‑4439990

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: