Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,399,985 (C)7→8 coding (252/516 nt) ogt ← O‑6‑alkylguanine‑DNA:cysteine‑protein methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,399,9781.C97.2% 121.5 / ‑3.2 36coding (259/516 nt)ogtO‑6‑alkylguanine‑DNA:cysteine‑protein methyltransferase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (17/18);  total (17/19)
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.79e-01

GGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAA‑TGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑T‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCC  >  NC_000913/1399842‑1400111
                                                                                                                                           |                                                                                                                                     
ggATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAA‑TGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑T‑cc                                                                                                                                     <  2:206579/139‑1 (MQ=255)
       aTCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAA‑TGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑T‑ccccccccg                                                                                                                              <  1:1029468/139‑3 (MQ=255)
       aTCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAA‑TGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑T‑ccccccccg                                                                                                                              >  1:116132/1‑137 (MQ=255)
       aTCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAA‑TGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑T‑ccccccccg                                                                                                                              <  2:821313/139‑3 (MQ=255)
          cGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAA‑TGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑T‑ccccccccgtag                                                                                                                           <  2:1281857/139‑6 (MQ=255)
           gTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAA‑TGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TCCCCCCCCGTAGc                                                                                                                           <  1:736316/139‑1 (MQ=255)
              tGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAA‑TGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TCCCCCCCCGTAGCAGt                                                                                                                        >  1:409478/1‑139 (MQ=255)
                 cGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAA‑TGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑T‑ccccccccg                                                                                                                              <  2:346319/129‑3 (MQ=255)
                 cGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAA‑TGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑T‑ccccccccg                                                                                                                              >  1:346319/1‑127 (MQ=255)
                  gCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAA‑TGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TCCCCCCCCGTAGCAGTGGGa                                                                                                                    <  1:579001/139‑1 (MQ=255)
                  gCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAA‑TGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TCCCCCCCCGTAGCAGTGGGa                                                                                                                    <  2:1385383/139‑1 (MQ=255)
                       aaCGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAA‑TGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TCCCCCCCCGTAGCAGTGGGAAGCGt                                                                                                               <  2:1255710/139‑1 (MQ=255)
                              cgcgccgcgcCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAA‑TGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TCCCCCCCCGTAGCAGTGGGAAGCGTATCaata                                                                                                        >  2:818890/1‑139 (MQ=255)
                                   cgggccggaaCGGCCCAATTGCTAAGCCATTTGG‑CGTAATTGCATTCCCTGCCGCAAGGGGTTAGGGCGTATTTTTTTCCAGACTTCGCGCTGAATTGCCGCC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCt                                                                                                   <  1:959767/130‑1 (MQ=255)
                                                cccAATTGCTCAGCCAGTTGGCCGTAA‑TGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCa                                                                                      >  2:180056/1‑139 (MQ=255)
                                                cccAATTGCTCAGCCAGTTGGCCGTAA‑TGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑T‑cc                                                                                                                                     <  1:1009892/91‑1 (MQ=255)
                                                cccAATTGCTCAGCCAGTTGGCCGTAA‑TGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑T‑cc                                                                                                                                     >  2:1009892/1‑91 (MQ=255)
                                                          cAGCCAGTTGGCCGTAA‑TGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCAc                                                                            >  2:641848/1‑139 (MQ=255)
                                                                                 tACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGcc                                                      <  1:180056/139‑1 (MQ=255)
                                                                                      gCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTGAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGAt                                                 <  2:834706/139‑1 (MQ=255)
                                                                                      gCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGAt                                                 >  1:534804/1‑139 (MQ=255)
                                                                                      gCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑T‑ccccccccg                                                                                                                              >  1:471464/1‑59 (MQ=255)
                                                                                      gCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑T‑ccccccccg                                                                                                                              <  2:471464/61‑3 (MQ=255)
                                                                                       cccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGAtt                                                >  1:1359763/1‑139 (MQ=255)
                                                                                        ccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATtg                                               >  2:1366219/1‑139 (MQ=255)
                                                                                           cAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATtggtg                                            >  2:318962/1‑139 (MQ=255)
                                                                                           cAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATtggtg                                            >  2:1429106/1‑139 (MQ=255)
                                                                                            aGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATtggtgg                                           <  1:318962/139‑1 (MQ=255)
                                                                                             ggggATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TCCCCCCCCGTAGCAGTGGGAAGc                                                                                                                 >  1:538349/1‑68 (MQ=255)
                                                                                             ggggATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TCCCCCCCCGTAGCAGTGGGAAGc                                                                                                                 <  2:538349/68‑1 (MQ=255)
                                                                                                     tGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑T‑cccc                                                                                                                                   >  2:870532/1‑41 (MQ=255)
                                                                                                     tGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑T‑cccc                                                                                                                                   <  1:870532/41‑1 (MQ=255)
                                                                                                       cGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATgcgc                                <  1:1429106/139‑1 (MQ=255)
                                                                                                          aGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGc                                                                                                    <  1:1162898/68‑1 (MQ=255)
                                                                                                          aGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGc                                                                                                    >  2:1162898/1‑68 (MQ=255)
                                                                                                          aGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCa                             >  2:1019142/1‑139 (MQ=255)
                                                                                                            tgtTTTCCAGACTTCGCGCTGAAATGGCG‑TCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATa                           >  2:198796/1‑139 (MQ=255)
                                                                                                               tttCCAGACTTCGCGCTGAAATGGCG‑TCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATg                                   >  2:1440156/1‑128 (MQ=255)
                                                                                                               tttCCAGACTTCGCGCTGAAATGGCG‑TCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATg                                   <  1:1440156/128‑1 (MQ=255)
                                                                                                                 tCCAGACTTCGCGCTGAAATGGCG‑TCCCCCCCCGTAGCAGTggg                                                                                                                     >  2:211584/1‑44 (MQ=255)
                                                                                                                 tCCAGACTTCGCGCTGAAATGGCG‑TCCCCCCCCGTAGCAGTggg                                                                                                                     <  1:211584/44‑1 (MQ=255)
                                                                                                                    aGACTTCGCGCTGAAATGGCG‑TCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCcttctt                   <  1:641848/139‑1 (MQ=255)
                                                                                                                        ttCGCGCTGAAATGGCG‑TCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATGGCCTTCTTTGCg               >  2:1404218/1‑139 (MQ=255)
                                                                                                                          cgcgCTGAAATGGCG‑TCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTaa                                                            <  1:837571/92‑1 (MQ=255)
                                                                                                                          cgcgCTGAAATGGCG‑TCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTaa                                                            >  2:837571/1‑92 (MQ=255)
                                                                                                                              cTGAAATGGCG‑TCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATGGCCTTCTTTGCGATAATg         <  1:1404218/139‑1 (MQ=255)
                                                                                                                                     ggCG‑TCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTcc  <  2:409478/139‑1 (MQ=255)
                                                                                                                                     ggCG‑TCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTcc  <  1:1366219/139‑1 (MQ=255)
                                                                                                                                     ggCG‑TCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTcc  <  1:1176539/139‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                     
GGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAA‑TGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑T‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCC  >  NC_000913/1399842‑1400111

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