Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,229,323 C→G 11.6% intergenic (+258/+115) yohP → / ← dusC uncharacterized protein/tRNA‑dihydrouridine synthase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,229,3230CG11.6% 140.2 / 9.6 43intergenic (+258/+115)yohP/dusCuncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base C (23/15);  new base G (2/3);  total (25/18)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCTGAA  >  NC_000913/2229185‑2229466
                                                                                                                                          |                                                                                                                                               
gcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTg                                                                                                                                       >  1:155748/1‑149 (MQ=32)
  cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGcc                                                                                                                                     >  1:163045/1‑149 (MQ=32)
     tCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATGGATGTGGGCTga                                                                                                                                      >  2:287347/1‑144 (MQ=12)
     tCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATGGATGTGGGCTga                                                                                                                                      <  1:287347/145‑2 (MQ=12)
          gccgccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATGGATGTgg                                                                                                                                           <  3:297725/135‑1 (MQ=17)
                     cACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGcc                                                                                                                  >  2:501373/1‑149 (MQ=32)
                             ggCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCa                                                                                                                           >  2:524707/1‑132 (MQ=32)
                             ggCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCa                                                                                                                           <  1:524707/132‑1 (MQ=32)
                                            aaTCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCaaa                                                                                           <  2:155748/149‑1 (MQ=35)
                                                                cTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGAt                                                                       <  2:163045/149‑1 (MQ=255)
                                                                 tCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATa                                                                      >  4:599989/1‑149 (MQ=255)
                                                                    gCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAag                                                                   >  2:227176/1‑149 (MQ=255)
                                                                     cTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAaga                                                                  >  1:342773/1‑149 (MQ=255)
                                                                      ttATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAagag                                                                 <  2:391883/149‑1 (MQ=255)
                                                                      ttATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAagag                                                                 >  1:133871/1‑149 (MQ=255)
                                                                      ttATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAagag                                                                 <  2:133871/149‑1 (MQ=255)
                                                                      ttATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAagag                                                                 >  1:391883/1‑149 (MQ=255)
                                                                          cccTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCaa                                                                                                           >  1:16976/1‑103 (MQ=31)
                                                                          cccTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCaa                                                                                                           <  2:16976/103‑1 (MQ=31)
                                                                          cccTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTaa                                                             >  2:389142/1‑149 (MQ=255)
                                                                           ccTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAac                                                            >  3:558255/1‑149 (MQ=255)
                                                                            cTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAaca                                                           >  2:611925/1‑149 (MQ=255)
                                                                            cTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAaca                                                           >  2:605788/1‑149 (MQ=255)
                                                                             ttACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTCCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAacac                                                          >  4:221243/1‑149 (MQ=255)
                                                                                    tACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAaacaa                                                                                        <  2:567375/112‑1 (MQ=38)
                                                                                    tACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAaacaa                                                                                        >  1:567375/1‑112 (MQ=38)
                                                                                            cGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACAt                                           <  3:599989/149‑1 (MQ=255)
                                                                                             gggCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATc                                          >  3:314616/1‑149 (MQ=255)
                                                                                                   gAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGc                                                                <  4:270329/121‑1 (MQ=255)
                                                                                                   gAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGc                                                                >  3:270329/1‑121 (MQ=255)
                                                                                                        aaCATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGtttt                               >  1:392426/1‑149 (MQ=255)
                                                                                                        aaCATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGtttt                               <  1:501373/149‑1 (MQ=255)
                                                                                                        aaCATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGtttt                               >  2:613664/1‑149 (MQ=255)
                                                                                                              ccTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTagag                         >  2:157071/1‑149 (MQ=255)
                                                                                                                  tAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTtc                     >  4:221923/1‑149 (MQ=255)
                                                                                                                    aaaCGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTtctc                   <  1:389142/149‑1 (MQ=255)
                                                                                                                        gcgcCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGata               <  4:314616/149‑1 (MQ=255)
                                                                                                                              tCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAAtt         <  1:227176/149‑1 (MQ=255)
                                                                                                                                 gCCCACATGGATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTtct                    <  2:469560/135‑1 (MQ=255)
                                                                                                                                 gCCCACATGGATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTtct                    >  1:469560/1‑135 (MQ=255)
                                                                                                                                     acaTCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCTGaa  <  1:605788/149‑1 (MQ=255)
                                                                                                                                      caTCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATg                                    <  4:213799/114‑1 (MQ=255)
                                                                                                                                      caTCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATg                                    >  3:213799/1‑114 (MQ=255)
                                                                                                                                          |                                                                                                                                               
GCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCTGAA  >  NC_000913/2229185‑2229466

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: