Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,007,693 G→A *610* (TAG→TAA)  recQ → ATP‑dependent DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,007,6930GA100.0% 129.7 / NA 41*610* (TAG→TAA) recQATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/21);  total (20/21)

GATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACT‑CCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTCCCGATTTCTT  >  NC_000913/4007558‑4007827
                                                                                                                                       |                                                                                                                                       
gaTGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGACGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATca                                                                                                                                      >  2:1262916/1‑139 (MQ=255)
             aTGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATaaaaaa                                                                                                                         >  1:232706/1‑139 (MQ=255)
               gCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGt                                                                                                                       >  1:21698/1‑139 (MQ=255)
                ccGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTg                                                                                                                      <  1:1114750/139‑1 (MQ=255)
                ccGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTg                                                                                                                      <  1:2480/139‑1 (MQ=255)
                  gATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATACGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                                    >  2:176115/1‑139 (MQ=255)
                        cGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                                    <  2:993643/133‑1 (MQ=255)
                        cGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                                    >  1:993643/1‑133 (MQ=255)
                          ccAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcag                                                                                                                                     >  1:1090036/1‑114 (MQ=255)
                          ccAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcag                                                                                                                                     <  2:1090036/114‑1 (MQ=255)
                            aGCGAAATGCTCAGCGTTAACGGCGTTGGGATACGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGa                                                                                                          <  1:176115/139‑1 (MQ=255)
                                         gCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACT‑CCGTAAACGtt                                                                                            <  2:21698/139‑1 (MQ=255)
                                                    gTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCAGGTAT‑AAGAATCCCGTAAAGGtctgtccttattc                                                                                 >  1:114538/1‑127 (MQ=255)
                                                     ttGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACT‑CCGTAAACGTTTCCCCCGCGAGt                                                                                >  1:343271/1‑139 (MQ=255)
                                                     ttGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACT‑CCGTAAACGTTTCCCCCGCGAGt                                                                                <  2:232706/139‑1 (MQ=255)
                                                     ttGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACT‑CCGTAAACGTTTCCCCCGCGAGt                                                                                <  1:876338/139‑1 (MQ=255)
                                                        ggATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACT‑CCGTAAACGTTTCCCCCGCGAGTCaa                                                                             <  1:202899/139‑1 (MQ=255)
                                                             cgcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGa                                                                                                          >  2:1315638/1‑106 (MQ=255)
                                                             cgcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGa                                                                                                          <  1:1315638/106‑1 (MQ=255)
                                                              gcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACT‑CCGTAAACGTTTCCCCCGCGAGTCAAatgtat                                                                       <  2:1225301/139‑1 (MQ=255)
                                                                  gCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAg                                                                                                                        >  2:388761/1‑87 (MQ=255)
                                                                  gCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAg                                                                                                                        <  1:388761/87‑1 (MQ=255)
                                                                               ggCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACT‑CCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCa                                                      <  2:343271/139‑1 (MQ=255)
                                                                                      cGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGAc                                                                                                         <  1:436240/82‑1 (MQ=255)
                                                                                      cGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGAc                                                                                                         >  2:436240/1‑82 (MQ=255)
                                                                                            tGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACT‑CCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGt                                         <  1:618388/139‑1 (MQ=255)
                                                                                            tGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACT‑CCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGt                                         <  1:1262916/139‑1 (MQ=255)
                                                                                                  tGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGa                                                                                                          <  2:194885/69‑1 (MQ=255)
                                                                                                  tGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGa                                                                                                          >  1:194885/1‑69 (MQ=255)
                                                                                                  tGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACT‑CCGTAAACGTTTCCCCCGCGAGTCAAat                                                                           <  2:1006239/99‑1 (MQ=255)
                                                                                                  tGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACT‑CCGTAAACGTTTCCCCCGCGAGTCAAat                                                                           >  1:1006239/1‑99 (MQ=255)
                                                                                                              aTGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACT‑CCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTAt                       >  2:582309/1‑139 (MQ=255)
                                                                                                                tttGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACT‑CCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGa                     >  1:828136/2‑139 (MQ=255)
                                                                                                                          aTGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACT‑CCGTAAACGTTTCCCCcgc                                                                                    <  2:1310088/66‑1 (MQ=255)
                                                                                                                          aTGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACT‑CCGTAAACGTTTCCCCcgc                                                                                    >  1:1310088/1‑66 (MQ=255)
                                                                                                                           tGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACT‑CCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCAc                                                     <  1:1188135/96‑1 (MQ=255)
                                                                                                                           tGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACT‑CCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCAc                                                     >  2:1188135/1‑96 (MQ=255)
                                                                                                                            gACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACT‑CCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTCCCg         >  2:380719/1‑139 (MQ=255)
                                                                                                                             aCGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACT‑CCGTAAACGTTTCCCCCGCGAGt                                                                                <  2:233826/67‑1 (MQ=255)
                                                                                                                             aCGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACT‑CCGTAAACGTTTCCCCCGCGAGt                                                                                >  1:233826/1‑67 (MQ=255)
                                                                                                                                 agagTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACT‑CCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTCCCGATTTc    >  2:715656/1‑139 (MQ=255)
                                                                                                                                   agtaatCAGCAGCATAAAAAAGTGCCAGTATGAAGACT‑CCGTAAACGTTTCCCCCTCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTCCCGATTTCtt  <  1:715656/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
GATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACT‑CCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTCCCGATTTCTT  >  NC_000913/4007558‑4007827

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