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,9990GA98.2% 171.0 / NA 56A12T (GCC→ACC) yhbQGIY‑YIG nuclease superfamily protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (30/25);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGTTTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGCCAGTGGG  >  NC_000913/3299873‑3300129
                                                                                                                              |                                                                                                                                  
gttgttTTGCCAGCCAACGGCTGATGGCGATGTGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGc                                                                                                                        >  2:1101617/1‑139 (MQ=255)
 ttgttTTGCCAGCCAACGGCTGATGGCGATGTGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCt                                                                                                                       <  1:1101617/139‑1 (MQ=255)
 ttgttTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACa                                                                                                                              >  2:957120/1‑132 (MQ=255)
 ttgttTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACa                                                                                                                              <  1:957120/132‑1 (MQ=255)
                       aTGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAAt                                                                                                                            <  1:305550/112‑1 (MQ=255)
                       aTGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAAt                                                                                                                            >  2:305550/1‑112 (MQ=255)
                        tGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGCTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATg                                                                                                >  1:191437/1‑139 (MQ=255)
                            gATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCggg                                                                                                             >  2:1171622/1‑122 (MQ=255)
                            gATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCggg                                                                                                             <  1:1171622/122‑1 (MQ=255)
                                  tgtTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAAcgccg                                                                                      <  2:191437/139‑1 (MQ=255)
                                    tttCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATAcc                                                                                                                >  1:205894/1‑111 (MQ=255)
                                    tttCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATAcc                                                                                                                <  2:205894/111‑1 (MQ=255)
                                          ttGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATcagc                                                                              >  2:1101401/1‑139 (MQ=255)
                                                               tgGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCt                                                                                                                       >  1:578721/1‑77 (MQ=255)
                                                               tgGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCt                                                                                                                       <  2:578721/77‑1 (MQ=255)
                                                               tgGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAgggg                                                         >  1:245821/1‑139 (MQ=255)
                                                               tgGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAgggg                                                         >  1:86779/1‑139 (MQ=255)
                                                                 ggCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACgg                                                                                                   >  2:99924/1‑95 (MQ=255)
                                                                 ggCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACgg                                                                                                   <  1:99924/95‑1 (MQ=255)
                                                                 ggCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAgg                                                           <  2:140055/135‑1 (MQ=255)
                                                                 ggCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAgg                                                           >  1:140055/1‑135 (MQ=255)
                                                                 ggCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCg                                                       >  2:875323/1‑139 (MQ=255)
                                                                       aCCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATc                                                                                                          <  1:82003/82‑1 (MQ=255)
                                                                       aCCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAAt                                                                                                                            <  1:1309349/64‑1 (MQ=255)
                                                                       aCCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAAt                                                                                                                            >  2:1309349/1‑64 (MQ=255)
                                                                       aCCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATc                                                                                                          >  2:82003/1‑82 (MQ=255)
                                                                         cTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCaa                                                                      <  1:869724/116‑1 (MQ=255)
                                                                         cTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCaa                                                                      >  2:869724/1‑116 (MQ=255)
                                                                          ttAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTg                                              <  1:875323/139‑1 (MQ=255)
                                                                             aCATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCAc                                                                         <  1:1138196/109‑1 (MQ=255)
                                                                             aCATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCAc                                                                         >  2:1138196/1‑109 (MQ=255)
                                                                             aCATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTgcgc                                           <  1:1101401/139‑1 (MQ=255)
                                                                                  cGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCggg                                                                                                             <  2:580236/68‑1 (MQ=255)
                                                                                  cGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCggg                                                                                                             >  1:580236/1‑68 (MQ=255)
                                                                                   gACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGaaa                                     >  1:222546/1‑139 (MQ=255)
                                                                                    aCGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAg                                    >  1:834278/1‑139 (MQ=255)
                                                                                       cAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATAc                                                                                                                 <  1:1245465/59‑1 (MQ=255)
                                                                                       cAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATAc                                                                                                                 >  2:1245465/1‑59 (MQ=255)
                                                                                       cAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGgag                                 >  1:247429/1‑139 (MQ=255)
                                                                                       cAGACGATGACAACCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGgag                                 >  1:578584/1‑139 (MQ=255)
                                                                                             aTGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTaa                           <  2:86779/139‑1 (MQ=255)
                                                                                               gACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATaa                                                                                                                          <  1:758438/42‑1 (MQ=255)
                                                                                               gACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATaa                                                                                                                          >  2:758438/1‑42 (MQ=255)
                                                                                               gACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACg                         <  2:834278/139‑1 (MQ=255)
                                                                                                    ccTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACgg                                                                                                   >  1:115827/1‑60 (MQ=255)
                                                                                                    ccTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACgg                                                                                                   <  2:115827/60‑1 (MQ=255)
                                                                                                    ccTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCggg                                        >  2:870960/1‑119 (MQ=255)
                                                                                                    ccTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCggg                                        <  1:870960/119‑1 (MQ=255)
                                                                                                    ccTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATATCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGgt                    >  1:799596/1‑139 (MQ=255)
                                                                                                      tGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCgg                                         >  1:969938/1‑116 (MQ=255)
                                                                                                      tGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCgg                                         <  2:969938/116‑1 (MQ=255)
                                                                                                               tACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCgcgc         <  2:247429/139‑1 (MQ=255)
                                                                                                               tACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCgcgc         <  2:222546/139‑1 (MQ=255)
                                                                                                                      tCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGCCAGTggg  >  1:938068/1‑139 (MQ=255)
                                                                                                                          tACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTg                                              >  1:1270297/1‑91 (MQ=255)
                                                                                                                          tACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTg                                              <  2:1270297/91‑1 (MQ=255)
                                                                                                                           accaccGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGCCAGt     <  2:984560/131‑1 (MQ=255)
                                                                                                                              |                                                                                                                                  
GTTGTTTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGCCAGTGGG  >  NC_000913/3299873‑3300129

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: