Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,441,038 G→A intergenic (‑254/‑20) 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,441,0380GA81.0% 43.7 / ‑1.7 22intergenic (‑254/‑20)ydbK/ydbJputative pyruvate‑flavodoxin oxidoreductase/DUF333 domain‑containing protein YdbJ
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (6/11);  minor base T (4/0);  total (10/12)
Fisher's exact test for biased strand distribution p-value = 3.51e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTGAAGATCAGAGGAAG‑AAAA‑GATGCGAGCAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGT‑AGTAGTGAACCTGTTCAGCAGGCGACTGCCGCGCACGTAGCG‑CCAGGTTTAAAAGCGTCGATGTCCAGTAGTGGAGAAGCAAATTGTGCAATGATCGGCGGTTCGCTTTCTGTTG‑CCCGTCAACTGGATGGTACGGCGATTGGGATGTGTGCATTACCCAACGGCAAACGCTGT  >  NC_000913/1440836‑1441282
                                                                                                                                                                                                          |                                                                                                                                                                                                                                                          
cAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGATCAGAGGAAG‑AAAA‑GATGCGAGCAGCG‑TTTTGGGTAGGGTGTGc                                                                                                                                                                                                         <  1:108668/251‑1 (MQ=255)
cAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGATCAGAGGAAG‑AAAA‑GATGCGAGCAGCG‑TTTTGGGTAGGGTGTGc                                                                                                                                                                                                         <  2:201824/251‑1 (MQ=255)
              cAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTTGCTTTAAGATCAGAGGAAG‑AAAA‑AAAGCGGGCAGCG‑TTTTGGGGAggggggggcgctttttttttgt                                                                                                                                                                                           >  1:273404/1‑232 (MQ=255)
                aTTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGAAATGAATGTTTTGATCAAACAGAGGGCAAACAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTTAAGATCAGAGGAAGGAAAA‑G‑TGCGGGCAACG‑CTTTTGGTAGGGTGTGCCGCTTTAttgtttt                                                                                                                                                                                           >  1:192977/1‑247 (MQ=255)
                     aGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCAATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGGTTGGCTTAAAGATCAGAGGAAG‑AAAAAGA‑GCGAGCCACG‑GTTTGGGGGGGGGGTGCCGCTTTATTTTTTTCGGCgtg                                                                                                                                                                                    >  1:196900/1‑251 (MQ=255)
                                ttcttcTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGATCAGAGGAAG‑AAAA‑GATGCGAGCAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGT‑AGTAGTGaa                                                                                                                                                                         <  2:273404/250‑1 (MQ=255)
                                      tGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTGAAGATCAGAGGAAG‑AAAA‑GATGCGAGCAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGT‑AGTAGTGAACCTGtt                                                                                                                                                                   <  2:91785/250‑1 (MQ=255)
                                         aTTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTTAAGATCAGAGGGAG‑AAAA‑GATGCGAGGAGCG‑GTTTGGGTAGGGGGTGGCGCTTTTTTGTTGTCGGGGGGTGAG‑AGTGAACCTGTTtagc                                                                                                                                                               >  1:252536/1‑251 (MQ=255)
                                                    tgtgTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGATCAGAGGAAG‑AAAA‑GATGCGAGCAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGT‑AGTAGTGAACCTGTTCAGCAGGCGACTGc                                                                                                                                                     <  1:291803/250‑1 (MQ=255)
                                                                        aatgaatgTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGATCAGAGGAAG‑AAAA‑GATGCGAGCAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGT‑AGTAGTGAACCTGTTCAGCAGGCGACTGcc                                                                                                                                                    <  1:17986/231‑1 (MQ=255)
                                                                          tgaatgTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGATCAGAGGAAG‑AAAA‑GATGCGAGCAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTTTTTTCGGCGGGT‑AGTAGTGAGCCTGGTTAGTAGGCGACTGCCGCGCACGCAGCG‑CCAGGtttt                                                                                                                              >  1:37285/1‑249 (MQ=255)
                                                                            aatgTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGATCAGAGGAAG‑AAAA‑GATGCGAGCAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGT‑AGTAGTGAACCTGTTCAGCAGGCGACTGCCGCGCACGTAGCG‑CCAGGTTTaaaa                                                                                                                           <  1:897/251‑1 (MQ=255)
                                                                                                       aaaTAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGATCAGAGGAAG‑AAAA‑AATGCGAGCAGCG‑TTTTGGGGAGGGGGTGCCGCTTTATTGTTGTCGGCGTGT‑AGTAGGGGACCTGTTCAGCAGGCGGCTGCCGCGCACGTAGCG‑CCAGGTTTAAAAGCGTCGATGTCCCGTAGTGGGGAAgg                                                                                                 >  1:224363/1‑249 (MQ=255)
                                                                                                        aaTAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCGCTCTGTGCTGCCTATTATGTATGGCGGTGTCGTTTGGCTTAAAGATCAGAGGAAT‑AAAA‑GATGCGAGCAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTGGTCGGCGTGT‑AGTAGTGAACCTGTTCAGCAGGCGACTGCGGCGCACGTAGCG‑CCAGGTTTAAAAGCGTCGATGTCCAGTAGTGGAGAAGCaa                                                                                               <  1:139046/251‑1 (MQ=255)
                                                                                                            gCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGATCAGAGGAAG‑AAAA‑GATGCGAGCAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGT‑AGTAGTGAACCTGTTCAGCAGGCGACTGCCGCGCACGTAGCG‑CCAGGTTTAAAAGCGTCGATGTCCAGTAGTGGAGAAGCAAATtg                                                                                           <  1:279534/251‑1 (MQ=255)
                                                                                                            accATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGATCAGAGGAAG‑AAAA‑GATGCGAGCAGCG‑TTTTGGGTAGGGTGTTCCGCTTTATTGTTGTCGTCGTGT‑AGTAGTGAACCTGTTCAGCAGGCGACTGCCGCGCACGTAGCG‑CCAGGTTTAAAAGCGTCGATGTCCAGTAGTGGAGAAGCAAATtg                                                                                           <  1:266121/250‑1 (MQ=255)
                                                                                                                                    gCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTTAAGATCAGAGGAAG‑AAAA‑GATGCGGGGAGCG‑TTTTGGGGAGGGGGTGCCGCTTTATTGTTTTCGGCGTGG‑AGTAGTGAAACTGGTTAGCAGGCGACTGCCGCGCACGGAGCG‑CCCGGGTTAAAAAAGTCGGAGTTCCGTAGAGGGGAAGGAAAATGTGCAATGagtggcgggtcgggttt                                                                   >  1:287757/1‑235 (MQ=255)
                                                                                                                                    gCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGATCAGAGGAGG‑AAAA‑GATGCGGGCAGCGGTTTTGGG‑AGGGGGTGCCGCTTTATTGTTGTCGGGGGGG‑AGTAGGGAACCTGTTCAGCAGGCGGCTGCCGCGCACGGAGCGGCCGGGTTAAAAAGCG‑CGGTGGTCAGGAGGGGGGGAGCAAATTTTGCAATTATCGGCGGGTTGCtttt                                                                   >  1:112962/1‑250 (MQ=255)
                                                                                                                                                                        gctgcCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGATCAGAGGAAG‑AAAA‑GATGCGAGCAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGT‑AGTAGTGAACCTGTTCAGCAGGCGACTGCCGCGCACGTAGCG‑CCAGGTTTAAAAGCGTCGATGTCCAGTAGTGGAGAAGCAAATTGTGCAATGATCGGCGGTTCGCTTTCTGTTG‑CCCGTCAACTGGATGGTACGGCGATTGGGAt                              <  1:74533/251‑1 (MQ=255)
                                                                                                                                                                                                  tttGGCTTAAAGATCAGCGGAAG‑CAAA‑GATGCGAGCAGCG‑TTTTGGGTAGGGCGTGCCGCTTTATTGTTGTCGGCGTGT‑AGTAGTGCACCTGTTCAGCAGGCGACTGCCGCGCACGTAGCG‑CAAGGTTTAAAAGAGTCGTTGTCCAGTAGTGGAGAAGCAAATTGTGCAATGCTCGGCGGTTCGCTTTCTTTTGACCTGTCAGTTGG‑TGGTACGGGCATTCGGTTGTTTGCATTACCCAACgtcacaaccc    >  1:203953/1‑242 (MQ=255)
                                                                                                                                                                                                  tttGGCTTAAAGATCAGAGGAAG‑AAAA‑GATGCGAGCAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGT‑AGTAGTGAACCTGTTCAGCAGGCGACTGCCGCGCACGTAGCG‑CCAGGTTTAAAAGCGTCGATGTCCAGTAGTGGAGAAGCAAATTGTGCAATGATCGGCGGTTCGCTTTTTTTTG‑CCCGTCAACTGGGTGGGACGGCGGTTTGGGTGTGTGACTTACCCAAcagccaaagcg    >  1:42696/1‑241 (MQ=255)
                                                                                                                                                                                                     ggCTTAAAGATCAGAGGAAG‑AAAA‑GATGCGAGCAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGT‑AGTAGTGAACCTGTTCAGCAGGCGACTGCCGCGCACGTAGCG‑CCAGGTTTAAAAGCGTCGATGTCCAGTAGTGGAGAAGCAAATTGTGCAATGATCGGCGGTTCGCTTTCTGTTG‑CCCGTCAACTGGATGGTACGGCGATTGGGATGTGTGCATTACCCAACGGCAAACGCTGt  <  1:153576/250‑1 (MQ=255)
                                                                                                                                                                                                          |                                                                                                                                                                                                                                                          
CAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTGAAGATCAGAGGAAG‑AAAA‑GATGCGAGCAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGT‑AGTAGTGAACCTGTTCAGCAGGCGACTGCCGCGCACGTAGCG‑CCAGGTTTAAAAGCGTCGATGTCCAGTAGTGGAGAAGCAAATTGTGCAATGATCGGCGGTTCGCTTTCTGTTG‑CCCGTCAACTGGATGGTACGGCGATTGGGATGTGTGCATTACCCAACGGCAAACGCTGT  >  NC_000913/1440836‑1441282

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