Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,261,358 T→A 19.6% 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,3580TA19.6% 151.0 / 30.3 61intergenic (+52/‑311)yjbM/dusAuncharacterized protein/tRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base T (29/20);  new base A (6/6);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

CCAACCTCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATTTCC  >  NC_000913/4261223‑4261500
                                                                                                                                       |                                                                                                                                              
ccAACCTCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCtt                                                                                                                                   >  1:609843/1‑149 (MQ=255)
  aaCCTCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCtttt                                                                                                                                 <  3:496542/149‑1 (MQ=255)
      tCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCt                                                                                                                                    <  3:305739/142‑1 (MQ=255)
      tCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCt                                                                                                                                    >  4:305739/1‑142 (MQ=255)
              agagATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAAt                                                                                                                     >  1:210818/1‑149 (MQ=255)
               gagaTTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATa                                                                                                                    <  2:609843/149‑1 (MQ=255)
               gagaTTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATAAATATTGGGTCCtttt                                                                                                                                 <  3:288544/136‑1 (MQ=255)
               gagaTTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATAAATATTGGGTCCtttt                                                                                                                                 >  4:288544/1‑136 (MQ=255)
                   ttAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTAt                                                                                                                <  2:45488/149‑1 (MQ=255)
                      aaaTACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAg                                                                                                             >  4:552943/1‑149 (MQ=255)
                        atgcTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAACa                                                                                                           >  3:458996/4‑149 (MQ=255)
                            ttACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAAt                                                                                                                     >  3:666700/1‑135 (MQ=255)
                            ttACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAAt                                                                                                                     <  4:666700/135‑1 (MQ=255)
                            ttACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGga                                                                                                       >  2:651483/1‑149 (MQ=255)
                            ttACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGga                                                                                                       >  1:238832/1‑149 (MQ=255)
                                aGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                                                                                   >  2:209810/1‑149 (MQ=255)
                                 gACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCtt                                                                                                                                   <  4:391057/116‑1 (MQ=255)
                                 gACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCtt                                                                                                                                   >  3:391057/1‑116 (MQ=255)
                                      ttCTAAATTAATTTGTGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCAc                                                                                             >  1:518951/1‑149 (MQ=255)
                                            attaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGat                                                                                                    >  3:96515/1‑136 (MQ=255)
                                            attaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGat                                                                                                    <  4:96515/136‑1 (MQ=255)
                                            attaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCttttat                                                                                                                               >  4:664894/1‑107 (MQ=38)
                                            attaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCttttat                                                                                                                               <  3:664894/109‑3 (MQ=38)
                                             ttaattTATGAGTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGtt                                                                                      >  1:514450/1‑149 (MQ=255)
                                                   tATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgt                                                                                <  3:552943/149‑1 (MQ=255)
                                                       aTTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGa                                                                            >  4:14508/1‑149 (MQ=255)
                                                        tttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAt                                                                           >  1:123570/1‑149 (MQ=255)
                                                            ttCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAttttt                                                                       <  2:321287/149‑1 (MQ=255)
                                                              cGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTttta                                                                     <  1:609249/149‑1 (MQ=255)
                                                               ggggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTCTCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttat                                                                    >  3:445356/1‑149 (MQ=255)
                                                                   ttAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTaa                                                                                   >  1:265638/1‑130 (MQ=255)
                                                                   ttAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTaa                                                                                   <  2:265638/130‑1 (MQ=255)
                                                                     aaaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCa                                                                                                           >  2:34205/1‑104 (MQ=255)
                                                                     aaaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCa                                                                                                           <  1:34205/104‑1 (MQ=255)
                                                                          ggAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGaa                                                         >  3:168816/1‑149 (MQ=255)
                                                                             gaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCt                                                      <  1:651483/149‑1 (MQ=255)
                                                                                ataataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGat                                                                                                    >  1:615197/1‑100 (MQ=255)
                                                                                ataataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGat                                                                                                    <  2:615197/100‑1 (MQ=255)
                                                                                   atatGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGCGTGATTTTTTATTTATCGTCGAACCTGGATTg                                                >  3:103862/1‑149 (MQ=255)
                                                                                       gtgAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTa                                            <  2:238832/149‑1 (MQ=255)
                                                                                        tgAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTAt                                           >  1:328798/1‑149 (MQ=255)
                                                                                          aaGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCa                                         <  2:328798/149‑1 (MQ=255)
                                                                                                ttcGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGcc                                   >  3:134219/1‑149 (MQ=255)
                                                                                                     tAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAAc                              >  2:419792/1‑149 (MQ=255)
                                                                                                        tATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACaaa                           <  2:514450/149‑1 (MQ=255)
                                                                                                          tGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGt                         >  1:419560/1‑149 (MQ=255)
                                                                                                              aGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAAc                     >  4:353177/1‑149 (MQ=255)
                                                                                                                gCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACgg                   <  4:458996/149‑1 (MQ=255)
                                                                                                                       ataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCtaataa            >  4:492158/1‑149 (MQ=255)
                                                                                                                        taAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAg           >  3:510376/1‑149 (MQ=255)
                                                                                                                         aaaaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGa          <  2:518951/149‑1 (MQ=255)
                                                                                                                         aaaaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGa          >  4:18272/1‑149 (MQ=255)
                                                                                                                         aaaaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGa          >  4:292691/1‑149 (MQ=255)
                                                                                                                          aaaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGat         >  4:43117/1‑149 (MQ=255)
                                                                                                                          aaaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGat         <  2:645083/149‑1 (MQ=255)
                                                                                                                          aaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATCTCACCAGAGTTCAATGTGTGATTTTTTATTTATCGTCGAACCTGGAATGTTTATCATTGGCCTTaa                               <  4:478056/127‑1 (MQ=255)
                                                                                                                          aaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTaa                               >  3:478056/1‑127 (MQ=255)
                                                                                                                             gaaCCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTgg                                     <  1:246114/116‑1 (MQ=255)
                                                                                                                             gaaCCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTgg                                     >  2:246114/3‑118 (MQ=255)
                                                                                                                              gACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGattatt     <  3:14508/149‑1 (MQ=255)
                                                                                                                                 ccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATTTcc  <  3:292691/149‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                              
CCAACCTCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATTTCC  >  NC_000913/4261223‑4261500

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