Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 339,856 T→C intergenic (+113/‑309) yahG → / → yahI DUF1116 family protein/carbamate kinase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913339,8560TC92.7% 115.0 / 1.1 42intergenic (+113/‑309)yahG/yahIDUF1116 family protein/carbamate kinase‑like protein
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (16/22);  total (17/25)
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.80e-01

ATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGC  >  NC_000913/339762‑339951
                                                                                              |                                                                                               
aTGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACgcggc                                                                                            <  2:324502/100‑1 (MQ=255)
         gAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGcc                                                                                   <  1:953578/100‑1 (MQ=255)
          ccgatcTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCt                                                                                  <  2:568062/95‑1 (MQ=255)
          aaCGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGTCGGACGCGGCGTGAACGCCt                                                                                  <  2:794355/100‑1 (MQ=255)
          aaCGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCt                                                                                  <  1:932945/100‑1 (MQ=255)
          aaCGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCt                                                                                  <  2:280999/100‑1 (MQ=255)
               cTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATcc                                                                             <  1:636934/100‑1 (MQ=255)
                 tATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCgg                                                                           >  1:97816/1‑100 (MQ=255)
                  aTCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTCCCGGACGCGGCGTGAACGCCTTATcc                                                                             >  2:87650/1‑97 (MQ=255)
                        cttACGAATAGCGCAAGAATCTGTAGGCCTGATAAGTGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGa                                                                    <  2:337917/98‑1 (MQ=11)
                          tACGAATGGCGCAAGAATCTTTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCAGCGTGAACGCCTTATCCGGCCTACGa                                                                    <  2:337865/98‑1 (MQ=255)
                          tACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAAt                                                                  <  1:1297363/100‑1 (MQ=255)
                           aCGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATg                                                                 <  1:162027/100‑1 (MQ=255)
                              aaTGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGgcg                                                              >  2:1264721/1‑100 (MQ=255)
                                 ggCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGgcg                                                              >  2:1076892/1‑97 (MQ=38)
                                       aGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCGAGAATCt                                                     >  2:340387/1‑100 (MQ=9)
                                          aTCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTa                                                  >  1:698617/1‑100 (MQ=255)
                                             tGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGc                                               <  1:566659/100‑1 (MQ=255)
                                             tGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGAAGGc                                               >  2:278964/1‑100 (MQ=255)
                                               tAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCt                                             <  1:637670/100‑1 (MQ=255)
                                                 ggCCGGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGa                                           <  2:903292/100‑1 (MQ=255)
                                                   ccTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATa                                         <  1:558139/100‑1 (MQ=255)
                                                      gATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCGAGAATCTGAAGGCCGGCTAAGc                                      >  2:444842/1‑100 (MQ=11)
                                                      gATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAg                                       >  1:1050699/1‑99 (MQ=255)
                                                        naaGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGt                                    >  1:385087/2‑100 (MQ=255)
                                                           gCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAgc                                 >  2:686355/1‑100 (MQ=255)
                                                           gCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAgc                                 <  1:1264732/100‑1 (MQ=255)
                                                             gTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAgcgc                               >  1:1088762/1‑100 (MQ=255)
                                                             gTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAgcgc                               <  1:779121/100‑1 (MQ=255)
                                                                 cgcATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCa                           >  1:964991/1‑100 (MQ=255)
                                                                  gcATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAg                          <  1:87652/100‑1 (MQ=255)
                                                                        ggCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAttt                    <  1:1228367/100‑1 (MQ=255)
                                                                           aTTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGcc                 <  1:990928/100‑1 (MQ=255)
                                                                               gTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTGGGCCTGTTAACCGTAGGGCATCAGGCATTTGccacc              >  2:1133486/1‑99 (MQ=14)
                                                                                tCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCcaccac            <  1:686359/100‑1 (MQ=255)
                                                                                 cACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACt           <  1:104687/100‑1 (MQ=255)
                                                                                    cATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGcc        >  1:1186182/1‑100 (MQ=255)
                                                                                     aTTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGccacc              >  2:754308/1‑93 (MQ=25)
                                                                                       tGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCg       <  2:875027/98‑1 (MQ=255)
                                                                                       ggCCGGATAAGGCGTTCACGCCGCGTCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCAGGa     >  2:1119951/2‑100 (MQ=11)
                                                                                        gCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGACATTTGCCACCACTGCCGGAt    <  1:1195970/100‑1 (MQ=255)
                                                                                          cGGAGGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATgc  <  2:1203810/100‑1 (MQ=255)
                                                                                              |                                                                                               
ATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGC  >  NC_000913/339762‑339951

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: