Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,147,052 T→C A163A (GCA→GCG)  alkA ← DNA‑3‑methyladenine glycosylase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,147,0520TC100.0% 143.7 / NA 41A163A (GCA→GCG) alkADNA‑3‑methyladenine glycosylase 2
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/17);  total (24/17)

GCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCACGCTCACCAGTTGGCCTAAAATCGCCCGC  >  NC_000913/2146916‑2147187
                                                                                                                                        |                                                                                                                                       
gCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTctgccgct                                                                                                                                      <  1:10614/140‑1 (MQ=255)
   ttttCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGcc                                                                                                                                   <  1:321061/140‑1 (MQ=255)
                cTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCg                                                                                                                      >  2:684464/1‑140 (MQ=255)
                cTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCg                                                                                                                      >  2:587172/1‑140 (MQ=255)
                  ccACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGcc                                                                                                                                   <  2:82687/125‑1 (MQ=255)
                  ccACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGcc                                                                                                                                   >  1:82687/1‑125 (MQ=255)
                      aTCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGa                                                                                                                <  1:684464/140‑1 (MQ=255)
                      aTCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGa                                                                                                                <  2:918286/140‑1 (MQ=255)
                        cGCCCGGTATTGTCATTGGTAAGGTGCCCACCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGTGAAg                                                                                                              <  1:358022/140‑1 (MQ=255)
                                     cATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCg                                                                                                 <  2:951742/140‑1 (MQ=255)
                                         ggTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGaaa                                                                                             >  2:592743/1‑140 (MQ=255)
                                         ggTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGaaa                                                                                             >  2:452094/1‑140 (MQ=255)
                                         ggTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGaaa                                                                                             >  2:382135/1‑140 (MQ=255)
                                                tGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATcc                                                                                      >  2:358107/1‑140 (MQ=255)
                                                  cccTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAatcatc                                                                                       <  1:1081072/137‑1 (MQ=255)
                                                  cccTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAatcatc                                                                                       >  2:1081072/1‑137 (MQ=255)
                                                  cccTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAg                                                                                    >  1:925100/1‑140 (MQ=255)
                                                  cccTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAg                                                                                    <  1:358107/140‑1 (MQ=255)
                                                       cAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTTCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGtt                                                                               <  2:737993/140‑1 (MQ=255)
                                                       cAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGtt                                                                               <  1:845740/140‑1 (MQ=255)
                                                       cAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGtt                                                                               <  1:592743/140‑1 (MQ=255)
                                                         gcgcCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCg                                                                             >  2:175363/1‑140 (MQ=255)
                                                             cgcATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCAt                                                                         <  2:925100/140‑1 (MQ=255)
                                                              gcATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATa                                                                        >  2:198456/1‑140 (MQ=255)
                                                              gcATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATa                                                                        >  2:392264/1‑140 (MQ=255)
                                                               cATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATaa                                                                       >  2:726420/1‑140 (MQ=255)
                                                               cATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATaa                                                                       <  1:762881/140‑1 (MQ=255)
                                                               cATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATaa                                                                       >  2:843489/1‑140 (MQ=255)
                                                                                       ctGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAAt                                               >  2:1098256/1‑140 (MQ=255)
                                                                                           ccGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTg                                           <  1:1149705/140‑1 (MQ=255)
                                                                                              tttCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTgccg                                        >  2:1045093/1‑140 (MQ=255)
                                                                                               ttCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTgccgc                                       >  2:120251/1‑140 (MQ=255)
                                                                                                 cAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCa                                     <  2:1156411/140‑1 (MQ=255)
                                                                                                    cGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCg                                  <  1:452094/140‑1 (MQ=255)
                                                                                                      gCATACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGcc                                <  1:1045093/140‑1 (MQ=255)
                                                                                                        aTACCTAACGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCAc                              >  2:1108990/1‑140 (MQ=255)
                                                                                                                cGCTTTTAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCTGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCACGCTCACCa                      >  2:62402/1‑140 (MQ=255)
                                                                                                                     ttAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCACGCTCACCAGTTgg                 >  2:1168374/1‑140 (MQ=255)
                                                                                                                     ttAATGCCTGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCACGCTCACCAGTTgg                 >  2:797812/1‑140 (MQ=255)
                                                                                                                             tGCGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCACGCTCACCAGTTGGCCTAAAAt         >  1:833777/1‑140 (MQ=255)
                                                                                                                               cGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCACGCTCACCAGTTGGCCTAAAATCg       >  2:375153/1‑140 (MQ=255)
                                                                                                                               cGGGTCTGCCGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCACGCTCACCAGTTGGCCTAAAATCg       >  2:963856/1‑140 (MQ=255)
                                                                                                                                    ctgccgctgCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCACGCTCACCAGTTGGCCTAAAATCGCCCGc  <  1:587172/140‑1 (MQ=255)
                                                                                                                                    ctgccgctgCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCACGCTCACCAGTTGGCCTAAAATCGCCCGc  <  1:390043/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                       
GCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCACGCTCACCAGTTGGCCTAAAATCGCCCGC  >  NC_000913/2146916‑2147187

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