Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,476,586 T→C 13.3% intergenic (‑146/+21) yheO ← / ← fkpA putative PAS domain‑containing DNA‑binding transcriptional regulator/FKBP‑type peptidyl‑prolyl cis‑trans isomerase (rotamase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,476,5860TC13.3% 120.5 / 10.0 45intergenic (‑146/+21)yheO/fkpAputative PAS domain‑containing DNA‑binding transcriptional regulator/FKBP‑type peptidyl‑prolyl cis‑trans isomerase (rotamase)
Reads supporting (aligned to +/- strand):  ref base T (19/20);  new base C (4/2);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

CATGAAAAAAAACACTCCTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGA  >  NC_000913/3476438‑3476730
                                                                                                                                                    |                                                                                                                                                
cATGAAAAAAAACACTCCTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGt                                                                                                                                                  <  2:304019/149‑1 (MQ=255)
 aTGAAAAAAAACACTCCTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTc                                                                                                                                                 <  1:487358/149‑1 (MQ=255)
 aTGAAAAAAAACACTCCTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTc                                                                                                                                                 <  1:883535/149‑1 (MQ=255)
  tGAAAAAAAACACTCCTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCa                                                                                                                                                <  1:803474/149‑1 (MQ=255)
      aaaaaaCACTCCTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAggcg                                                                                                                                            >  1:448872/1‑149 (MQ=255)
        aaaaCACTCCTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAggcggc                                                                                                                                          <  2:536342/149‑1 (MQ=255)
             acTCCTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGACCAGGCGGCGGttt                                                                                                                                     <  3:497614/149‑2 (MQ=255)
               tCCTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCtt                                                                                                                                   >  1:682190/1‑149 (MQ=255)
                 cTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTaa                                                                                                                                 <  2:682190/149‑1 (MQ=255)
                    ttCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGc                                                                                                                              <  4:781198/149‑1 (MQ=255)
                                 cGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCa                                                                                                                 <  3:381419/149‑1 (MQ=255)
                                   ttAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGa                                                                                                               >  1:46534/1‑149 (MQ=255)
                                        ttttCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTg                                                                                                          >  1:722817/1‑149 (MQ=255)
                                           tCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGACCAGgcggcg                                                                                                                                         <  1:319850/115‑1 (MQ=255)
                                           tCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGACCAGgcggcg                                                                                                                                         >  2:319850/1‑115 (MQ=255)
                                                    gCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCAt                                                                                              <  2:722817/149‑1 (MQ=255)
                                                     cTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATc                                                                                             >  2:179475/1‑149 (MQ=255)
                                                      tCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCa                                                                                            >  4:859526/1‑149 (MQ=255)
                                                            atgatgCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGc                                                                                                         >  2:387915/1‑130 (MQ=255)
                                                            atgatgCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGc                                                                                                         <  1:387915/130‑1 (MQ=255)
                                                                 gCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCtt                                                                                 >  2:724378/1‑149 (MQ=255)
                                                                   gggTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCtt                                                                                                    <  3:769015/128‑1 (MQ=255)
                                                                   gggTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCtt                                                                                                    >  4:769015/1‑128 (MQ=255)
                                                                        aCTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCa                                                                          <  1:724378/149‑1 (MQ=255)
                                                                          tAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGc                                                                        >  2:234412/1‑149 (MQ=255)
                                                                          tAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGc                                                                        <  4:258812/149‑1 (MQ=255)
                                                                          tAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGc                                                                        >  2:576213/1‑149 (MQ=255)
                                                                                           ggTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTca                                                       <  1:179475/149‑1 (MQ=255)
                                                                                              tCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACAt                                                    <  2:448872/149‑1 (MQ=255)
                                                                                                        aGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGctc                                          >  3:872227/1‑149 (MQ=255)
                                                                                                           aCTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCcagca                                              <  4:574157/142‑1 (MQ=255)
                                                                                                           aCTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCcagca                                              >  3:574157/1‑142 (MQ=255)
                                                                                                           aCTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTAc                                       <  4:872227/149‑1 (MQ=255)
                                                                                                                ttATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCaa                                  >  3:853154/1‑149 (MQ=255)
                                                                                                                 tataTCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCaaa                                 <  2:642065/149‑1 (MQ=255)
                                                                                                                     tCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACAcc                             >  2:755226/1‑149 (MQ=255)
                                                                                                                      ccGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCa                            <  1:576213/149‑1 (MQ=255)
                                                                                                                       cGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAg                           <  1:226630/149‑1 (MQ=255)
                                                                                                                                  ttaagaaCCGCCGCCTGACCAGGCGGCGGTTTTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCg                >  4:778189/6‑149 (MQ=255)
                                                                                                                                  aaaaaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGTCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCg                >  2:331733/1‑149 (MQ=255)
                                                                                                                                  aaaaaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCg                <  1:755226/149‑1 (MQ=255)
                                                                                                                                    aaaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGt              >  2:529196/1‑149 (MQ=255)
                                                                                                                                        aCCGCCGCCTGACCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGa          >  3:479077/1‑149 (MQ=255)
                                                                                                                                          cgccgcCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAACCACCAGGGTCGAATTCGGTGGGATc        >  3:339355/1‑149 (MQ=255)
                                                                                                                                            ccgccTGACCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGGTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATccc      >  3:278044/1‑149 (MQ=255)
                                                                                                                                                cTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGa  >  1:122078/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CATGAAAAAAAACACTCCTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGA  >  NC_000913/3476438‑3476730

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: ‑