Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,261,358 T→A 16.0% intergenic (+52/‑311) yjbM → / → dusA uncharacterized protein/tRNA‑dihydrouridine synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,261,3580TA16.0% 61.9 / 6.5 25intergenic (+52/‑311)yjbM/dusAuncharacterized protein/tRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base T (10/11);  new base A (2/2);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01

ATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATTTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACAT  >  NC_000913/4261175‑4261627
                                                                                                                                                                                       |                                                                                                                                                                                                                                                                             
aTGCGAAATGCTTAACCAAAATTAATTAGNTTCGCGATGTCGTTTGACCCAACCTNAAATGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGNGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTg                                                                                                                                                                                    <  3:90291/275‑1 (MQ=255)
       aTGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGaa                                                                                                                                                                                        >  1:90291/1‑264 (MQ=255)
                                                                               ctgACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTCGGTTCTTTTTTGTCTATCAATACTATTAGCAGGgaga                                                                                                                                                                                                                                    <  7:958854/146‑1 (MQ=255)
                                                                                   ctcTTCTAAATTAATTTATGATTCGTTCTGGGTTAAATCGGAGAAATAATATGTTAAGTTCTTCGATCGTATGGAAGGGATTATATAAAAGGACCCAATATTTAGTGGGTTATTTTTTCTCTATCAATACTATTAGCAGGGAGataag                                                                                                                                                                                                                                >  8:958854/3‑146 (MQ=255)
                                                                                     cTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtattta                                                                                                                                                                                                                                                                            >  8:823002/1‑102 (MQ=255)
                                                                                     cTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtattta                                                                                                                                                                                                                                                                            <  7:823002/102‑1 (MQ=255)
                                                                                                        tttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatttat                                                                                                                                                                                                                                                                           >  8:736396/1‑84 (MQ=255)
                                                                                                        tttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatttat                                                                                                                                                                                                                                                                           <  7:736396/84‑1 (MQ=255)
                                                                                                                    tAAATCGGAGAAATAATATGTGCAGTTCTTCGATAGTATGGAAGGCAGTATAAATAAGTGCCCATGATTTATTAGGTTCTTATATCTCTCTCAATAGTATTAGCGGGGTGGTATATCAGGAGAGATTAAGCTGTGATTTtttattttt                                                                                                                                                                                               <  8:106051/148‑3 (MQ=255)
                                                                                                                           gagaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGc                                                                                                                                                                                                                                           <  7:587407/97‑1 (MQ=255)
                                                                                                                           gagaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGc                                                                                                                                                                                                                                           >  8:587407/1‑97 (MQ=255)
                                                                                                                                ataaaaTGTGAAGTTTTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATTTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTa                                                     <  2:491665/274‑1 (MQ=255)
                                                                                                                                                                    tatataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCCCCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTCATAAGATTATTTCCATCACTTCGTTAGAGCTTCATGCCTTGAGAGGATCTCATTTTTATTTTGCAATGAGTCAGGCGCTTCCTGTTGTTATGGTATAGTTCTCCGATATTGACCCTCCTGCATtgtg                >  4:486487/1‑275 (MQ=255)
                                                                                                                                                                    tatataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTaa                                                                                                                                                                                                                  <  7:1114376/81‑1 (MQ=255)
                                                                                                                                                                    tatataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTaa                                                                                                                                                                                                                  >  8:1114376/1‑81 (MQ=255)
                                                                                                                                                                       ataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttttt                                                                                                                                                                                                     <  7:304949/91‑1 (MQ=255)
                                                                                                                                                                       ataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttttt                                                                                                                                                                                                     >  8:304949/1‑91 (MQ=255)
                                                                                                                                                                          aaaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGt                                                                                                                                                                              >  7:451406/1‑111 (MQ=255)
                                                                                                                                                                          aaaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGt                                                                                                                                                                              <  8:451406/111‑1 (MQ=255)
                                                                                                                                                                              aaCCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcg                                                                                                                                                                                             <  8:830833/91‑1 (MQ=255)
                                                                                                                                                                              aaCCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcg                                                                                                                                                                                             >  7:830833/2‑92 (MQ=255)
                                                                                                                                                                               aCCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATc                                                                                                                                                                                                                              <  7:147092/58‑1 (MQ=255)
                                                                                                                                                                               aCCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATc                                                                                                                                                                                                                              >  8:147092/1‑58 (MQ=255)
                                                                                                                                                                                cccAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAttt                                                                                                                                                                                                        <  6:1063305/79‑1 (MQ=255)
                                                                                                                                                                                cccAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAttt                                                                                                                                                                                                        >  5:1063305/1‑79 (MQ=255)
                                                                                                                                                                                  ccatatTTATTGGGTTCTTTTTTCTCAATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATTTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTAAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAataacat  <  2:486488/273‑1 (MQ=255)
                                                                                                                                                                                  cAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtc                                                                                                                                                                                              <  5:713886/87‑1 (MQ=255)
                                                                                                                                                                                  cAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtc                                                                                                                                                                                              >  6:713886/1‑87 (MQ=255)
                                                                                                                                                                                  cAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTa                                                                                                                                                               <  5:381507/118‑1 (MQ=255)
                                                                                                                                                                                  cAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTa                                                                                                                                                               >  6:381507/1‑118 (MQ=255)
                                                                                                                                                                                       |                                                                                                                                                                                                                                                                             
ATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATTTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACAT  >  NC_000913/4261175‑4261627

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