Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,767,214 C→A intergenic (+35/‑7) yibA → / → yibJ putative lyase containing HEAT‑repeat/putative RHS domain‑containing protein YibJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,767,2140CA96.8% 102.4 / ‑2.6 31intergenic (+35/‑7)yibA/yibJputative lyase containing HEAT‑repeat/putative RHS domain‑containing protein YibJ
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (18/12);  total (19/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

AGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGATGGACCCGGTGTACACG  >  NC_000913/3767080‑3767349
                                                                                                                                      |                                                                                                                                       
aGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCAGGCAc                                                                                                                                    <  1:1154/140‑1 (MQ=255)
                                 tATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCAGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAAt                                                                                                   <  2:26342/140‑1 (MQ=255)
                                 tATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCAGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAAt                                                                                                   >  2:322182/1‑140 (MQ=255)
                                   tgttgtACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCAGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATcc                                                                                                 >  1:128700/1‑140 (MQ=255)
                                   tgttgtACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCAGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATcc                                                                                                 >  2:166614/1‑140 (MQ=255)
                                    gttgtACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCAGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCt                                                                                                <  1:245584/140‑1 (MQ=255)
                                      tgtACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCAGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTgg                                                                                              <  1:227677/140‑1 (MQ=255)
                                        tACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCAGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCt                                                                                            <  1:329836/140‑1 (MQ=255)
                                        tACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCAGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCt                                                                                            >  2:11804/1‑140 (MQ=255)
                                          cAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCAGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGa                                                                                          >  1:72431/1‑140 (MQ=255)
                                              tttGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCAGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCgg                                                                                      <  1:322182/140‑1 (MQ=255)
                                              tttGATGACAATGAAATTATAACTTCCGCTATTGATAAGCCGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCAGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCgg                                                                                      <  1:25287/140‑1 (MQ=255)
                                               ttGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCAGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCggg                                                                                     >  2:283570/1‑140 (MQ=255)
                                                  atgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCAGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGgtga                                                                                  >  2:58661/1‑140 (MQ=255)
                                                  atgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCAGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGgtga                                                                                  >  1:173090/1‑140 (MQ=255)
                                                                  aaCTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCAGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGccgccg                                                                  <  1:247415/140‑1 (MQ=255)
                                                                    cTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCAGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAgc                                                                       >  2:227244/1‑133 (MQ=255)
                                                                    cTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCAGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAgc                                                                       <  1:227244/133‑1 (MQ=255)
                                                                    cTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCAGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCgct                                                                >  2:311142/1‑140 (MQ=255)
                                                                                       gAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCAGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCt                                             <  1:58661/140‑1 (MQ=255)
                                                                                       gAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCAGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCt                                             >  2:186650/1‑140 (MQ=255)
                                                                                        aaGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCAGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTg                                            >  2:274/1‑140 (MQ=255)
                                                                                          gCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCCTGCGGCCTGAc                                          >  2:337256/1‑140 (MQ=255)
                                                                                             ttCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCAGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGt                                       >  2:84305/1‑140 (MQ=255)
                                                                                               cATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCAGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTgg                                     >  1:280014/1‑140 (MQ=255)
                                                                                                             ccTATGCACAGTGGGGGATTCCCGCAGGCACGGTGCTGATGCTCGACCGGCTGCAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAAc                       <  2:115231/140‑1 (MQ=255)
                                                                                                                 tGCACAGTGGGGGATTCCCGCAGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATTAAAAACCAGa                   >  2:224936/1‑140 (MQ=255)
                                                                                                                  gCACAGTGGGGGATTCCCGCAGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAt                  >  1:53126/1‑140 (MQ=255)
                                                                                                                     caGTGGGGGATTCCCGCAGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGATGGa               <  1:274/140‑1 (MQ=255)
                                                                                                                       gTGGGGGATTCCCGCAGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGATGGAcc             <  1:84305/140‑1 (MQ=255)
                                                                                                                        tGGGGGATTCCCGCAGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGATGGAccc            >  1:174487/1‑140 (MQ=255)
                                                                                                                          ggggATTCCCGCAGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGATGGACCCgg          >  2:74394/1‑140 (MQ=255)
                                                                                                                                  ccgcagGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGATGGACCCGGTGTACACg  >  2:82868/1‑140 (MQ=37)
                                                                                                                                  ccgcagGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGAGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGGTGGACCCGGTTTaccca  >  2:336467/1‑137 (MQ=11)
                                                                                                                                      |                                                                                                                                       
AGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGATGGACCCGGTGTACACG  >  NC_000913/3767080‑3767349

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