Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,440,791 C→T 40.9% intergenic (‑7/‑267) ydbK ← / → ydbJ putative pyruvate‑flavodoxin oxidoreductase/DUF333 domain‑containing protein YdbJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,440,7910CT40.9% 18.1 / 35.0 43intergenic (‑7/‑267)ydbK/ydbJputative pyruvate‑flavodoxin oxidoreductase/DUF333 domain‑containing protein YdbJ
Reads supporting (aligned to +/- strand):  ref base C (14/10);  new base T (9/8);  total (25/18)
Fisher's exact test for biased strand distribution p-value = 7.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGTTTCCGGCCCAGGCATCAGCCTGTTCTGCCATCGTGGAACTGGGGGTAATAGGGTAGATGGCGATAACTTCACTGGTGCGAAATGCGACCGAAGCAACCGCGCCATTACCGTCAATAGTAATCATATGACACCCTTACATTGCGCAAATGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGA  >  NC_000913/1440656‑1440930
                                                                                                                                       |                                                                                                                                           
gCCGTTTCCGGCCCAGGCATCAGCCTGTTCTGCCATCGTGGAACTGGGGGTAATAGGGTAGATGACGTTAACTTCACTGGTGCGAAATGCGACCGAAGCAACCGCGCCATTACCGTCAATAGTAATCATATGACATACTTACa                                                                                                                                      >  1:344522/1‑143 (MQ=255)
 ccGTTTCCGGCCCAGGCATCAGCCTGTTCTGCCATCGTGGAACTGGGGGTAATAGGGTAGATGGAGATAACTTCACTGGTGCGAAATGCGACCGGAGCAAAGAGGCCATTCCCATCAATGATAAGCCTGTGACAACCCTACag                                                                                                                                     >  2:545469/1‑142 (MQ=255)
     ttCCGGCTCAGGCATCAGCCTGTTCTGCCATCGTGGAACTGGGGGGAATAGGGTGGATGGCGATAACTTCACTGGTGCGAAATGCGACCGAAGCAACCGCGCCATTACCTTCAATAGTAATCATATGACACCCTTACATTgcg                                                                                                                                 <  2:215347/143‑1 (MQ=255)
          gCCCAGGCATCAGCCTGTTCTGCCATCGTGGAACTGGGGGTAATAGGGTAGATGGCGATAACTTCACTGGTGCGAAATGCGACCGAAGCAACCGCGCCATTACCGTCAATAGTAATCATATGACACCCTTACATTGCGCAAAt                                                                                                                            >  1:190307/1‑143 (MQ=255)
                 taTCAGCCTGTTCTGCCATCGTGGAACTGAGGGTATTTTGTTTTGTGGTGATAACTTCACTGGTGCGAAATGCGACCGAAGCAACCGCGCCATTACCGTCAATAGTAATCATATGACATCCTTACATTGCGCAAATGAGGGgc                                                                                                                     <  2:112153/142‑1 (MQ=255)
                     aGCCTGTTCTGCCATCGTGGAACTGGGGGTAATAGGGTAGATGGCGATAACTTCACTGGTGCGAAATGCGACCGAAGCAACCGCGCCATTACCGTCAATAGTAATCATATGACACCCTTACATTGCGCAAATGAGGGGCGCAc                                                                                                                 >  2:149687/1‑143 (MQ=255)
                            tCTGCCATCGTGGGTCTGGGGGTAATAGGGTAGATGGCGATAACTTCACTGGTGCGAAATGCGACCGAAGCAATCGCGCCATTACCGTCAATAGTAATCATATGACATCCTTACATTGCGCAAATGAGGGGCGCACGAAATtg                                                                                                          <  2:442362/143‑1 (MQ=255)
                              tGCCATCGTGGAACTGGGGGTAATAGGGTAGATGGCGATAACCTCATTGGTGCGAAATGCGACCGAAGCAACCGCGCCATTACCGTCAATAGTAATCATATTATATCCTTACATTGCGCAAATGAGGGGCGCACGAAATtgct                                                                                                        <  1:134505/143‑1 (MQ=255)
                                  aTCGTGGAAGTGGGTGTAATAGGGTAGATGTCGATAACTTCACTGGTGCGAAATGCCACCGAAGCAACCGCGCCATTACCGTCAATAGGAATCATATGACACCCTTACATTGCGCAAATGAGGGGCGCACGAAATTGCTgcgc                                                                                                    <  2:350469/143‑1 (MQ=255)
                                   tCGTGGAACTGGGGGTAATAGGGTAGATGGCGATAACTTCACTGGTGCGAAATGCGACCGAAGCAACCGCGCCATTACCGTCAATAGTAATCATATGACATCCTTACATTGCGCAAATGAGGGGCGCACGAAATTGCTgcgcg                                                                                                   >  1:84515/1‑143 (MQ=255)
                                         aaCTGGGGGTAATAGGGTAGTTGGCGATAACTTCACTGGTGCGAAATGCGACCGAAGCAACCGCGCCATTACCGTCAATAGTAATCATATGACACCCTTACATTGCGCAAATGAGGGGCGCACGAAATTGCTGCGCGCCCagt                                                                                             >  1:94825/1‑143 (MQ=255)
                                         aaCTGGGGGTAATAGGGTAGATGGCGATAACTTCACTGGTGCGAAATGCGACCGAAGCAACCGCGCCATTACCGTCAATAGTAATCATATGACATCCTTACATTGCGCACATGAGGGGCGCACGAAATTGCTGCGCGCCCagt                                                                                             >  2:228673/1‑143 (MQ=255)
                                              ggggTAAGAGGGTAGATGGCGATAACTGCACTGGTGCGAAATGCGACCGAAGCAACCGCGCCATTACCGTCAATAGTAATCATATGACACCCTTACATTGCGCAAATGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTaa                                                                                        <  2:78771/143‑1 (MQ=255)
                                                  tAATAGGGTAGATGGCGATAACTTCACTGGTGCGAAATGCGACCGAAGCAACCGCGCCGTTACCGTCAATAATAATCATATGACATCCTTACATTGCGCAAATGAGGGGGGCACGAAACTGCTGCCGGCCCCGGAGCAATCtt                                                                                    >  2:542743/1‑143 (MQ=255)
                                                  tAATAGGGTAGATGGCGATAACTTCACTGGTGCGAAATGCGACCGAAGCAACCGCGCCATTACCGTCAATAGTAATCATATGACACCCTTACATTGCGCAAATGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCtt                                                                                    <  1:59050/143‑1 (MQ=255)
                                                  tAAGAGGGTAGCGGGCGATAACGTCACGGGTGCGAAATGCGACCGAAGCAACCGCGCCATTACCGTCAATAGTAATCATATGACATCCTTACATGGCGCAAATGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCtt                                                                                    <  1:357760/143‑1 (MQ=255)
                                                          tAGATGGCGATAACTTCACTGGTGCGAAATGCGACCGAAGCAACCGCGCCATTACCGTCAATAGTAATCATATGACACCCTTACATTGCGCAAATGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAAtttt                                                                            >  2:547108/1‑143 (MQ=255)
                                                              tGGCGATAACTTCACTGGTGCGAAATGCGACCGAAGCAACCGCGCCATTACCGTCAATAGTAATCATATGACATCCTTACATTGCGCAAATGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTGAGCa                                                                        >  1:410962/1‑143 (MQ=255)
                                                                 cGATAACTTCACTGGTGCGAACTGCGCCCGAAGCAACCGCGCCCTAACCGTGAATAGTCAGTATTTGATACCCTTACATTGCGCAAACGAGGGGCGCACGAACTTGCTGCTCGCCCAGTAGTAGTCGTTCAATTTTAGCaagt                                                                     >  1:177494/1‑141 (MQ=255)
                                                                          cACTGGTGCGAAATGCGACCGAAGCCACCGCGCCATTACCGTCAATAGTAATCATATGACATCCTTACATTGCGCAAATGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTttctt                                                            <  2:304887/143‑1 (MQ=255)
                                                                          cACTGGTGCGAAATGCGACCGAAGCAACCGCGCCATTACCGTCAAGAGTAATCATAGGACACCCTTACATTGCGCAAAGGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTttctt                                                            <  1:285806/143‑1 (MQ=255)
                                                                            cTGGTGCGAAATGCGACCGCAGCAACCGCGCCATTACCGTCAATAGTAATCATATGACACCCTTACAATGCGCATATGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTtcttct                                                          <  2:65167/143‑1 (MQ=255)
                                                                               gTGCGAAATGCGACCGAAGCAACCGCGCCATTACCGTCAATAGTCATCATATGACAGCCTTACATTGCGCAAATGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAAGTTTAGCAAATGGCTTGCTTCTGCa                                                       >  2:203832/1‑143 (MQ=255)
                                                                                tGCGAAATGCGACCGAAGCAACCGCGCCATTACCGTCAATAGTAATCATATGACATCCTTACATTGCGCAAATGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTCATCTTTCCATTTTAGCACACGGCTCGCTACTGCaa                                                      >  2:57942/1‑142 (MQ=255)
                                                                                   gAAATGCGACCGGAGCAACCGCGCCATTACCGTCAATAGTAATCATATGACACCCTTACATTGCGCAAATGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCCATTTTTGCAATTGGCTTTCTTCTGCAtttt                                                   <  2:333163/143‑1 (MQ=255)
                                                                                        gCGACCGAAGCAACCGCGCCATTACCGTCAATAGTAATCATATGACATCCTTACATTGCGCAAATGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCtt                                              <  2:80815/143‑1 (MQ=255)
                                                                                         cGACCGAAGCAACCGCGCCACTACCGTCAATGGCCATCATATGACATCCAAACATTGCGCAAATGAGGGGTGGACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCACATGGCTTTCTTCTGCATCTTCGCttt                                             >  1:479790/1‑143 (MQ=255)
                                                                                          gACCGAAGCAACCGCGCCATTACCGTCAATAGTAATCATATGACACCCTTACATTGCGCAAATGAGGGGCGCACGAAATTGGTGCGCGCCCAGTAGTAATCTTTCAATCTTAGCAAAGGGCTGCGTTCTGCATTTTCGGtttt                                            >  1:19332/1‑143 (MQ=255)
                                                                                            ccGAAGCAACCGCGCCATTACCGTCAATAGTAATCATATGACACCCTTACATTGCGCAAATGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATCTTCGCGTTtgt                                          >  2:69730/1‑143 (MQ=255)
                                                                                              gAAGCAACCGCGCCATTACCGTCAATAGTAATCATATGACGTCCGTACATTGCGCAAATGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCCTTTTCGCTTTtgtgg                                        >  1:205836/1‑142 (MQ=255)
                                                                                                      cgcgCCATTACCGTCAATAGTAATCATATGACATCCTTACATTGCGCAAATGAGACGCGCAAGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCACATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCGCcaca                                >  2:482944/1‑143 (MQ=255)
                                                                                                      cgccCCATTACCGTCAATGGTAATCATATGACATCCTTACATTGCGCAAATGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCcaca                                <  2:411943/143‑1 (MQ=255)
                                                                                                       gcgcCATTACCGTCAATAGTAATCATATGACACCCTTACATTGCGCAAATGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACAt                               >  2:161701/1‑143 (MQ=255)
                                                                                                            aTTACCGTCAATAGTAATCATATAACACCCTTACATTGCACAAATGAGGGGCGCACAAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTGTGCATTTTCGCTAGTGTGTCCCCCACAACAGCg                          >  2:84292/1‑143 (MQ=255)
                                                                                                            aTTACCCTCAATTGTTATGATATGACCTCCTTTCATTGCGCAAATGACGGGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTCTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCg                          <  2:471017/143‑1 (MQ=255)
                                                                                                                   tCAATAGTAATCATATGACACCCTTACATTGCGCAAATGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTaatgaa                   >  2:49743/1‑143 (MQ=255)
                                                                                                                       tAGTAATCATATGACACCCTTACATTGCGCAAATGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCa                             <  1:132252/129‑1 (MQ=255)
                                                                                                                       tAGTAATCATATGACACCCTTACATTGCGCAAATGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCa                             >  2:132252/1‑129 (MQ=255)
                                                                                                                          tAATCATATGACACCCTTACATTGCGCAAATGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTg            >  2:434341/1‑143 (MQ=255)
                                                                                                                          tAATCATATGACACCCTTACATTGCGCAAATGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTACCTGATGTGTCCCCCACATCAGCGTAATGAATGTTTTg            >  1:370125/1‑143 (MQ=255)
                                                                                                                            aTCAGATGACACCCTTACATTGCGCAAATGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGAt          <  1:482305/143‑1 (MQ=255)
                                                                                                                             tCATATGACACCCTTACATTGCGCAAATGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAAAGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACCTCAGCGTAATGAATTTTTTGATc         >  1:160217/1‑143 (MQ=255)
                                                                                                                                tatGACACCCTTACATTGCGCAAATGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCaaa      <  1:530835/143‑1 (MQ=255)
                                                                                                                                    acacCCTTACATTGCGCAAATGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACaga  <  1:417985/143‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                           
GCCGTTTCCGGCCCAGGCATCAGCCTGTTCTGCCATCGTGGAACTGGGGGTAATAGGGTAGATGGCGATAACTTCACTGGTGCGAAATGCGACCGAAGCAACCGCGCCATTACCGTCAATAGTAATCATATGACACCCTTACATTGCGCAAATGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGA  >  NC_000913/1440656‑1440930

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