Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,299,999 G→A A12T (GCC→ACC)  yhbQ → GIY‑YIG nuclease superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,299,9990GA96.0% 147.1 / ‑1.6 50A12T (GCC→ACC) yhbQGIY‑YIG nuclease superfamily protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (28/20);  total (29/21)
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.46e-01

GTTGTTTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCG  >  NC_000913/3299873‑3300121
                                                                                                                              |                                                                                                                          
gttgttTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAAt                                                                                                                    >  1:732146/1‑135 (MQ=255)
gttgttTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAAt                                                                                                                    <  2:732146/135‑1 (MQ=255)
                                         aTTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCt                                                                            <  1:811646/134‑1 (MQ=255)
                                         aTTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATcag                                                                       <  1:464778/139‑1 (MQ=255)
                                         aTTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTAGAACGCCGCt                                                                            >  2:811646/1‑134 (MQ=255)
                                                ttCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGaa                                                                                   >  1:683260/1‑120 (MQ=255)
                                                ttCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGaa                                                                                   <  2:683260/120‑1 (MQ=255)
                                                ttCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACa                                                                                                                      >  2:25444/1‑85 (MQ=255)
                                                ttCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACa                                                                                                                      <  1:25444/85‑1 (MQ=255)
                                                ttCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAAt                                                                                                                    <  1:602519/87‑1 (MQ=255)
                                                ttCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTa                                                                           <  2:642615/128‑1 (MQ=255)
                                                ttCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTa                                                                           >  1:642615/1‑128 (MQ=255)
                                                ttCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCAcc                                                                >  1:839103/1‑139 (MQ=255)
                                                ttCTCTTTTATACTGTGGGCCTGACCTTAACATGAGACGCAGACGATGAAACCCTGGTTTCTTTACTTGATCCGTACCACCGACAAt                                                                                                                    >  2:602519/1‑87 (MQ=255)
                                                  ctctTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCaa                                                              >  2:682351/1‑139 (MQ=255)
                                                          taCTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATAc                                                                                                         >  2:213152/1‑88 (MQ=255)
                                                          taCTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATAc                                                                                                         <  1:213152/88‑1 (MQ=255)
                                                             tgtgGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAgg                                                   >  1:408710/1‑139 (MQ=255)
                                                               tgGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTaccacc                                                                                                                          <  1:591069/66‑1 (MQ=255)
                                                               tgGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACa                                                                                                                      <  1:784230/70‑1 (MQ=255)
                                                               tgGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACa                                                                                                                      >  2:784230/1‑70 (MQ=255)
                                                               tgGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATa                                                                                                                   >  1:191708/1‑73 (MQ=255)
                                                               tgGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATa                                                                                                                   <  2:191708/73‑1 (MQ=255)
                                                               tgGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCt                                                                                                               <  1:437206/77‑1 (MQ=255)
                                                               tgGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCt                                                                                                               >  2:437206/1‑77 (MQ=255)
                                                               tgGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGa                                                                                                    <  1:47316/88‑1 (MQ=255)
                                                               tgGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGa                                                                                                    >  2:47316/1‑88 (MQ=255)
                                                               tgGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAAcg                                                                                 <  1:316291/107‑1 (MQ=255)
                                                               tgGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAAcg                                                                                 >  2:316291/1‑107 (MQ=255)
                                                               tgGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAgggg                                                 >  1:98373/1‑139 (MQ=255)
                                                               tgGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAgggg                                                 >  2:180690/1‑139 (MQ=255)
                                                               tgGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGCTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAgggg                                                 >  2:62521/1‑139 (MQ=255)
                                                                 ggCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCg                                               >  2:599506/1‑139 (MQ=255)
                                                                 ggCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCg                                               >  1:517549/1‑139 (MQ=255)
                                                                 ggCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCg                                               >  1:675410/1‑139 (MQ=255)
                                                                 ggCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCg                                               >  2:198913/1‑139 (MQ=255)
                                                                            aaCATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAAt                                                                                                                    <  2:487525/59‑1 (MQ=255)
                                                                            aaCATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAAt                                                                                                                    >  1:487525/1‑59 (MQ=255)
                                                                            aaCATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTgcg                                    <  2:675410/139‑1 (MQ=255)
                                                                               aTGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCgg                                 <  1:682351/139‑1 (MQ=255)
                                                                                  cGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGaa                              >  2:284159/1‑139 (MQ=255)
                                                                                               gACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACg                 <  2:408710/139‑1 (MQ=255)
                                                                                                 cacCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCt               >  1:495447/1‑139 (MQ=255)
                                                                                                    ccTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAg                            <  2:554499/123‑1 (MQ=255)
                                                                                                    ccTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAg                            >  1:554499/1‑123 (MQ=255)
                                                                                                        gTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGttt        >  2:867776/1‑139 (MQ=255)
                                                                                                           tCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTgcg                                    >  2:321199/1‑108 (MQ=255)
                                                                                                           tCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTgcg                                    <  1:321199/108‑1 (MQ=255)
                                                                                                            cTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCcg    <  1:867776/139‑1 (MQ=255)
                                                                                                             tttACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCcgc   >  2:889258/1‑139 (MQ=255)
                                                                                                              ttACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCcgcg  <  1:889258/139‑1 (MQ=255)
                                                                                                                           accaccGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCa                                                       >  1:396814/1‑73 (MQ=255)
                                                                                                                           accaccGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCa                                                       <  2:396814/73‑1 (MQ=255)
                                                                                                                              |                                                                                                                          
GTTGTTTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCG  >  NC_000913/3299873‑3300121

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