Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,229,323 C→G 14.3% 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,3230CG14.3% 145.7 / 17.3 49intergenic (+258/+115)yohP/dusCuncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base C (23/19);  new base G (4/3);  total (27/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.16e-01

GGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCTGAATA  >  NC_000913/2229181‑2229468
                                                                                                                                              |                                                                                                                                                 
ggcggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTgg                                                                                                                                             >  4:677395/1‑149 (MQ=32)
    gcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAt                                                                                                                                                 >  2:893818/1‑141 (MQ=32)
    gcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAt                                                                                                                                                 <  1:893818/141‑1 (MQ=32)
                      cTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAGCGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAAc                                                                                                                       >  3:1209568/1‑149 (MQ=17)
                                      tGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGa                                                                                                       <  4:1209568/149‑1 (MQ=35)
                                          gAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTggg                                                                                                                                            <  4:1236389/108‑1 (MQ=32)
                                          gAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTggg                                                                                                                                            >  3:1236389/1‑108 (MQ=32)
                                                                          ttATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAagag                                                                   <  3:677395/149‑1 (MQ=255)
                                                                                                         ggTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGt                                    >  1:245698/1‑149 (MQ=255)
                                                                                                          gTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAaga                                                                    >  1:1088971/1‑116 (MQ=255)
                                                                                                          gTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAaga                                                                    <  2:1088971/116‑1 (MQ=255)
                                                                                                            aaCATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCa                                           <  4:1464570/139‑1 (MQ=255)
                                                                                                            aaCATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCa                                           >  3:1464570/1‑139 (MQ=255)
                                                                                                             aCATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAaga                                                                    <  2:463606/113‑1 (MQ=255)
                                                                                                             aCATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAaga                                                                    >  1:463606/1‑113 (MQ=255)
                                                                                                             aCATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAag                                                                     >  3:919441/1‑112 (MQ=255)
                                                                                                             aCATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAag                                                                     <  4:919441/112‑1 (MQ=255)
                                                                                                               aTCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAagag                                                                   >  3:153700/1‑112 (MQ=255)
                                                                                                               aTCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAagag                                                                   <  4:153700/112‑1 (MQ=255)
                                                                                                                  ccTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTagag                           >  3:1361736/1‑149 (MQ=255)
                                                                                                                        aaaCGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAAAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTaga                            <  1:906059/142‑1 (MQ=255)
                                                                                                                        aaaCACGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTaga                            >  2:906059/1‑142 (MQ=255)
                                                                                                                          aCGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGa                   >  3:428358/1‑149 (MQ=255)
                                                                                                                            gcgcCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACAt                                             <  1:740013/121‑1 (MQ=255)
                                                                                                                            gcgcCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACAt                                             >  2:740013/1‑121 (MQ=255)
                                                                                                                            gcgcCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGata                 >  1:23031/1‑149 (MQ=255)
                                                                                                                                 cTCAGGCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGc                                                                  <  1:547039/95‑1 (MQ=255)
                                                                                                                                 cTCAGCCCACATCCATGTGGGCTTCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAAt            >  2:582033/1‑149 (MQ=255)
                                                                                                                                 cTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGc                                                                  >  2:547039/1‑95 (MQ=255)
                                                                                                                                 cTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAAt            <  2:23031/149‑1 (MQ=255)
                                                                                                                                  tCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACAt                                             >  2:154608/1‑115 (MQ=255)
                                                                                                                                  tCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACAt                                             <  1:154608/115‑1 (MQ=255)
                                                                                                                                    aGCCCACATGGATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGc         >  2:205223/1‑149 (MQ=255)
                                                                                                                                     gCCCACATGGATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGAc                                                                                                      >  4:716197/1‑55 (MQ=11)
                                                                                                                                     gCCCACATGGATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGAc                                                                                                      <  3:716197/55‑1 (MQ=14)
                                                                                                                                     gCCCACATGGATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATaaa                                                                      <  4:1363108/87‑1 (MQ=255)
                                                                                                                                     gCCCACATGGATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATaaa                                                                      >  3:1363108/1‑87 (MQ=255)
                                                                                                                                     gCCCACATGGATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAag                                                                     >  2:826243/1‑88 (MQ=255)
                                                                                                                                     gCCCACATGGATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAag                                                                     <  1:826243/88‑1 (MQ=255)
                                                                                                                                     gCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGcc        <  4:1361736/149‑1 (MQ=255)
                                                                                                                                      cccACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCt       <  2:245698/149‑1 (MQ=255)
                                                                                                                                       ccACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGttttt                                <  1:727344/123‑1 (MQ=255)
                                                                                                                                       ccACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGttttt                                >  2:727344/1‑123 (MQ=255)
                                                                                                                                       ccACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAAt            >  2:1054731/1‑143 (MQ=255)
                                                                                                                                       ccACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAAt            <  1:1054731/143‑1 (MQ=255)
                                                                                                                                         acaTCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGtttttt                               <  4:1253826/122‑1 (MQ=255)
                                                                                                                                         acaTCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGtttttt                               >  3:1253826/1‑122 (MQ=255)
                                                                                                                                           aTCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCTGAATa  >  1:1320536/1‑149 (MQ=255)
                                                                                                                                           aTCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCTGAATa  >  2:721590/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
GGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCTGAATA  >  NC_000913/2229181‑2229468

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