Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,139,692 T→C 12.1% intergenic (+99/+24) yghU → / ← hybG putative S‑transferase/hydrogenase 2 accessory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,139,6920TC12.1% 134.3 / 7.9 49intergenic (+99/+24)yghU/hybGputative S‑transferase/hydrogenase 2 accessory protein
Reads supporting (aligned to +/- strand):  ref base T (22/21);  new base C (4/2);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCC  >  NC_000913/3139562‑3139830
                                                                                                                                  |                                                                                                                                          
agaCGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTAt                                                                                                                                          <  3:506311/133‑1 (MQ=255)
agaCGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTAt                                                                                                                                          >  4:506311/1‑133 (MQ=255)
        aCGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTa                                                                                                                         >  4:246258/1‑142 (MQ=255)
        aCGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTa                                                                                                                         <  3:246258/142‑1 (MQ=255)
           gAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAat                                                                                                                       <  3:831591/141‑1 (MQ=255)
           gAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAat                                                                                                                       >  4:831591/1‑141 (MQ=255)
                 aaGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGGCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTa                                                                                                                         <  3:1063160/133‑1 (MQ=255)
                 aaGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCGGCCGCGTCTTATCAGGCCTACGAAAGCTa                                                                                                                         >  4:1063160/1‑133 (MQ=255)
                  aGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTaa                                                                                                                        >  4:285274/1‑133 (MQ=255)
                  aGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTaa                                                                                                                        <  3:285274/133‑1 (MQ=255)
                  aGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGt                                                                                                         >  1:125831/1‑148 (MQ=255)
                                 ggTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTaa                                                                                                                        >  1:341706/1‑118 (MQ=255)
                                 ggTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTaa                                                                                                                        <  2:341706/118‑1 (MQ=255)
                                         ttCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTAt                                                                                                                                          >  3:350453/1‑92 (MQ=255)
                                         ttCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTAt                                                                                                                                          <  4:350453/92‑1 (MQ=255)
                                          tcgtcgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATc                                                                                                                                         <  3:392142/92‑1 (MQ=255)
                                          tcgtcgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATc                                                                                                                                         >  4:392142/1‑92 (MQ=255)
                                             tcgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGc                                                                              >  3:1102766/1‑148 (MQ=255)
                                                 agcaATCCTTCAAGGCCGGATAAGGCGTTCACGCCGCATGCGGCAATCGTGCGCAATGCCTGATGCGGCGCTTCCTCTTCTTATCATGCGTACGCAAGCTAATCCTCACGCGCGGGTAACGCCCTAACCCAACTGGCGCAGtttgcct                                                                          <  2:125831/148‑7 (MQ=12)
                                                      gCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatcatc                                                                                                                   <  1:825321/102‑1 (MQ=255)
                                                      gCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatcatc                                                                                                                   >  2:825321/1‑102 (MQ=255)
                                                           ccAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTaa                                                                        >  3:498477/1‑140 (MQ=255)
                                                           ccAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTaa                                                                        <  4:498477/140‑1 (MQ=255)
                                                                ccGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGTACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTa                                                                         >  1:919617/1‑134 (MQ=255)
                                                                ccGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTa                                                                         <  2:919617/134‑1 (MQ=255)
                                                                      gcggCGTGAACGCCTTATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCAtt                                                                                         <  1:352184/110‑1 (MQ=37)
                                                                      gcggCGTGAACGCCTTATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCAtt                                                                                         >  2:352184/3‑112 (MQ=37)
                                                                          cGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGt                                                                                                 >  4:365441/1‑100 (MQ=255)
                                                                          cGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGt                                                                                                 <  3:365441/100‑1 (MQ=255)
                                                                                  ccgcATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtggctttg                                                              <  1:516947/127‑1 (MQ=255)
                                                                                  ccgcATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtggctttg                                                              >  2:516947/1‑127 (MQ=255)
                                                                                         cGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTa                                                                                                                                  >  3:113597/1‑52 (MQ=255)
                                                                                         cGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTa                                                                                                                                  <  4:113597/52‑1 (MQ=255)
                                                                                           gCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGt                                                                            >  1:814204/1‑104 (MQ=255)
                                                                                           gCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGt                                                                            <  2:814204/104‑1 (MQ=255)
                                                                                               tCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTtcgtcg                                               >  1:548650/1‑129 (MQ=255)
                                                                                               tCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTtcgtcg                                               <  2:548650/129‑1 (MQ=255)
                                                                                                      caATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACgtg                     >  1:358232/1‑148 (MQ=255)
                                                                                                       aaTGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACgtgt                    <  2:358232/148‑1 (MQ=255)
                                                                                                           ccTGATAAGACGCGGCAGCGTCGCATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCa                >  2:964888/1‑148 (MQ=255)
                                                                                                               ataagACGCGGCAGCGTCGCATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTTCGTCgatga                                           <  1:907713/113‑1 (MQ=255)
                                                                                                               ataagACGCGGCAGCGTCGCATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTTCGTCgatga                                           <  3:621921/113‑1 (MQ=255)
                                                                                                               ataagACGCGGCAGCGTCGCATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTTCGTCgatga                                           >  2:907713/5‑117 (MQ=255)
                                                                                                               ataagACGCGGCAGCGTCGCATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTTCGTCgatga                                           >  4:621921/5‑117 (MQ=255)
                                                                                                                   gACGCGGCAGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtggctt                                                                >  4:883582/1‑92 (MQ=255)
                                                                                                                   gACGCGGCAGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtggctt                                                                <  3:883582/92‑1 (MQ=255)
                                                                                                                       cTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCa                                                        >  3:820232/1‑96 (MQ=255)
                                                                                                                       cTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCa                                                        <  4:820232/96‑1 (MQ=255)
                                                                                                                         gcagcgTCGCATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGccc  >  1:908206/4‑148 (MQ=255)
                                                                                                                                  tATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtgg                                                                   >  3:513610/1‑74 (MQ=255)
                                                                                                                                  tATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtgg                                                                   <  4:513610/74‑1 (MQ=255)
                                                                                                                                  tATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtgg                                                                   <  2:553643/74‑1 (MQ=255)
                                                                                                                                  tATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTc                                      <  4:152213/103‑1 (MQ=255)
                                                                                                                                  tATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTc                                      >  3:152213/1‑103 (MQ=255)
                                                                                                                                  tATCAGGCCTACGAAAGCTAATCATCACGCGCCGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtgg                                                                   >  1:553643/1‑74 (MQ=255)
                                                                                                                                  |                                                                                                                                          
AGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCC  >  NC_000913/3139562‑3139830

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