Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,067,323 A→T 25.5% intergenic (‑150/+17) ygfI ← / ← yggE putative DNA‑binding transcriptional regulator/oxidative stress defense protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,067,3230AT25.5% 92.6 / 33.8 47intergenic (‑150/+17)ygfI/yggEputative DNA‑binding transcriptional regulator/oxidative stress defense protein
Reads supporting (aligned to +/- strand):  ref base A (20/15);  new base T (6/6);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 7.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.57e-01

CCTCCCCCGCATCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCG  >  NC_000913/3067177‑3067469
                                                                                                                                                  |                                                                                                                                                  
cctcccccGCATCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACg                                                                                                                                                  >  3:282090/1‑149 (MQ=255)
 ctcccccGCATCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGt                                                                                                                                                 <  4:282090/149‑1 (MQ=255)
      ccGCATCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTc                                                                                                                                                >  3:1102212/1‑145 (MQ=255)
      ccGCATCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTc                                                                                                                                                <  4:1102212/145‑1 (MQ=255)
              ttGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCgg                                                                                                                                              >  2:879568/1‑139 (MQ=255)
              ttGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCgg                                                                                                                                              <  2:1454609/139‑1 (MQ=255)
              ttGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCgg                                                                                                                                              >  1:1454609/1‑139 (MQ=255)
              ttGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCgg                                                                                                                                              <  1:879568/139‑1 (MQ=255)
                   cGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCggg                                                                                                                                             <  3:1282776/135‑1 (MQ=255)
                   cGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCggg                                                                                                                                             >  4:1282776/1‑135 (MQ=255)
                    gCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCg                                                                                                                                               <  1:1090986/132‑1 (MQ=255)
                    gCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCg                                                                                                                                               >  2:1090986/1‑132 (MQ=255)
                          ccGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGtgctgc                                                                                                                        >  4:1029802/1‑149 (MQ=255)
                           cGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCa                                                                                                                       >  4:378205/1‑149 (MQ=255)
                                     agCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAg                                                                                                             <  3:1029802/149‑1 (MQ=255)
                                        cTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGgtgt                                                                                                                  <  2:1123961/141‑1 (MQ=255)
                                        cTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGgtgt                                                                                                                  >  1:1123961/1‑141 (MQ=255)
                                              ttGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGGATTAGCCCGACGTATTCGTCgg                                                                                                                                              >  3:1382944/1‑107 (MQ=255)
                                              ttGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCgg                                                                                                                                              <  4:1382944/107‑1 (MQ=255)
                                              ttGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCgg                                                                                                                                              <  3:984393/107‑1 (MQ=255)
                                              ttGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCgg                                                                                                                                              >  4:984393/1‑107 (MQ=255)
                                                                    tAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGaa                                                                              <  3:378205/149‑1 (MQ=255)
                                                                               gTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAata                                                                                                                                    <  2:790532/84‑1 (MQ=255)
                                                                               gTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAata                                                                                                                                    >  1:790532/1‑84 (MQ=255)
                                                                                      ttttACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCggg                                                                                                                                             >  4:924077/1‑68 (MQ=255)
                                                                                      ttttACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCggg                                                                                                                                             <  3:924077/68‑1 (MQ=255)
                                                                                             gaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACcac                                                                        >  2:138035/1‑130 (MQ=255)
                                                                                             gaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACcac                                                                        <  1:138035/130‑1 (MQ=255)
                                                                                                       gaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCgg                                           >  1:1031232/1‑149 (MQ=255)
                                                                                                       gaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGAGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCgg                                           >  1:280780/1‑149 (MQ=255)
                                                                                                                                     tAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGGTGCAGGTGTTTTAGCg                                                                                                           >  4:1210399/1‑55 (MQ=255)
                                                                                                                                     tAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCg                                                                                                           <  3:1210399/55‑1 (MQ=255)
                                                                                                                                      aGCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTaa                                >  4:552056/1‑129 (MQ=255)
                                                                                                                                      aGCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTaa                                <  3:552056/129‑1 (MQ=255)
                                                                                                                                         ccGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgttgtt                                                                                                    >  4:378372/1‑58 (MQ=255)
                                                                                                                                         ccGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgttgtt                                                                                                    <  3:378372/58‑1 (MQ=255)
                                                                                                                                         ccGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATc                                                                                                <  2:1221831/62‑1 (MQ=255)
                                                                                                                                         ccGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATc                                                                                                >  1:1221831/1‑62 (MQ=255)
                                                                                                                                         ccGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGaa                                                                              <  3:1042571/80‑1 (MQ=255)
                                                                                                                                         ccGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGaa                                                                              >  4:1042571/1‑80 (MQ=255)
                                                                                                                                          cgacgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTcc                          <  1:854830/131‑1 (MQ=255)
                                                                                                                                          cgacgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTcc                          >  2:854830/1‑131 (MQ=255)
                                                                                                                                           gacgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCgg       >  2:1020879/1‑149 (MQ=255)
                                                                                                                                              gaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCAt    >  1:1111872/1‑149 (MQ=255)
                                                                                                                                                aTACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCt                         >  3:731437/1‑126 (MQ=255)
                                                                                                                                                aTACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCt                         <  4:731437/126‑1 (MQ=255)
                                                                                                                                                aTACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCg  >  3:293843/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CCTCCCCCGCATCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCG  >  NC_000913/3067177‑3067469

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