Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,462,758 G→A 25.0% A70T (GCA→ACA)  mag → DNA‑3‑methyladenine glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,462,7580GA25.0% 84.3 / 24.1 44A70T (GCA→ACA) magDNA‑3‑methyladenine glycosylase
Reads supporting (aligned to +/- strand):  ref base G (11/22);  new base A (5/6);  total (16/28)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

CTCAAACGTATACAGGCTACATCGTGGAAACGGAAGCTTACTTAGGTTTGAATGATCGTGCGGCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCCGCGCAATTGAACCAGAAGAAGGTTTATCCGCTATGTTCCGTAACAGAGGTAAGAAAGGCTACGA  >  CP000730/2462627‑2462897
                                                                                                                                   |                                                                                                                                           
cTCAAACGTATACAGGCTAAATCATGGAAACGGAAGCTTACTTAGGTTTGAATGATCGTGCGGCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCa                                                                                                                                    <  2:175147/141‑1 (MQ=255)
       gTATACAGGCTACATCGTGGAAACGGAAGCTTACTTAGGTTTGAATGATCGTGCGGCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATAc                                                                                                                             >  2:169248/1‑141 (MQ=255)
       gTATACAGGCTACATCGTGGAAACGGAAGCTTACTTAGGTTTGAATGATCGTGCGGCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATAc                                                                                                                             >  2:219250/1‑141 (MQ=255)
            cAGGCTACATCGTGGAAACGGAAGCTTACTTAGGTTTTAATGATCGTGCGGCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCAtt                                                                                                                        <  2:357182/141‑1 (MQ=255)
            cAGGCTACATCGTGGAAACGGAAGCTTACTTAGGTTTGAATGATCGTGCGGCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCAtt                                                                                                                        <  2:448806/141‑1 (MQ=255)
            cAGGCTACATCGTGGAAACGGAAGCTTACTTAGGTTTGAATGATCGTGCGGCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCAtt                                                                                                                        <  1:300236/141‑1 (MQ=255)
            cAGGCTACATCGTGGAAACGGAAGCTTACTTAGGTTTGAATGATCGTGCGGCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCAtt                                                                                                                        <  2:275732/141‑1 (MQ=255)
            cAGGCTACATCGTGGAAACGGAAGCTTACTTAGGTTTGAATGATCGTGCGGCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCAtt                                                                                                                        <  2:260501/141‑1 (MQ=255)
             aGGCTACATCGTGGAAACGGAAGCTTACTTAGGTTTGAATGATCGTGCGGCTCATGGCTATGGCGGTACAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCAttt                                                                                                                       <  1:365183/141‑1 (MQ=255)
             aGGCTACATCGTGGAAACGGAAGCTTACTTAGGTTTGAATGATCGTGCGGCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCAttt                                                                                                                       <  1:258360/141‑1 (MQ=255)
             aGGCTACATCGTGGAAACGGAAGCTTACTTAGGTTTGAATGATCGTGCGGCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCAttt                                                                                                                       <  1:269526/141‑1 (MQ=255)
             aGGCTACATCGTGGAAACGGAAGCTTACTTAGGTTTGAATGATCGTGCGGCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCAttt                                                                                                                       <  1:145397/141‑1 (MQ=255)
                    aTCGTGGAAACGGAAGCTTACTTAGGTTTGAATGATCGTGCGGCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCATTTACTCAtt                                                                                                                >  2:255383/1‑141 (MQ=255)
                      cGTGGAAACGGAAGCTTACTTAGGTTTGAATGATCGTGCGGCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCATTTACTCATTaa                                                                                                              >  1:330778/1‑141 (MQ=255)
                      cGTGGAAACGGAAGCTTACTTAGGTTTGAATGATCGTGCGGCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCATTTACTCATTaa                                                                                                              >  2:295759/1‑141 (MQ=255)
                      cGTGGAAACGGAAGCTTACTTAGGTTTGAATGATCGTGCGGCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCATTTACTCATTaa                                                                                                              >  2:379558/1‑141 (MQ=255)
                         ggAAACGGAAGCTTACTTAGGTTTGAATGATCGTGCGGCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCATTTACTCATTAAttt                                                                                                           <  1:323902/141‑1 (MQ=255)
                         ggAAACGGAAGCTTACTTAGGTTTGAATGATCGTGCGGCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCATTTACTCATTAAttt                                                                                                           <  1:400612/141‑1 (MQ=255)
                            aaCGGAAGCTTACTTAGGTTTGAATGATCGTGCGGCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCATTTACTCATTAATTTTGt                                                                                                        >  1:182835/1‑141 (MQ=255)
                                        cttaGGTTTGAATGATCGTGCGGCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGa                                                                                            <  1:219250/141‑1 (MQ=255)
                                                  aaTGATCGTGCGGCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATACACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCt                                                                                  >  2:243018/1‑141 (MQ=255)
                                                  aaTGATCGTGCGGCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATACACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCt                                                                                  >  1:348368/1‑141 (MQ=255)
                                                       tCGTGCGGCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAgg                                                                             >  2:97525/1‑141 (MQ=255)
                                                       tCGTGCGGCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAgg                                                                             >  2:44228/1‑141 (MQ=255)
                                                                    ggCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCcgc                                                                <  1:255383/141‑1 (MQ=255)
                                                                      cTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCCgcgc                                                              >  1:348934/1‑141 (MQ=255)
                                                                                     aaCACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCCGCGCAATTGAACCagaaga                                               <  1:284384/141‑1 (MQ=255)
                                                                                     aaCACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCCGCGCAATTGAACCagaaga                                               <  2:128678/141‑1 (MQ=255)
                                                                                     aaCACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCCGCGCAATTGAACCagaaga                                               <  1:155357/141‑1 (MQ=255)
                                                                                      acacCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATACACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCCGCGCAATTGAACCAgaagaa                                              <  1:243018/141‑1 (MQ=255)
                                                                                      acacCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATACACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCCGCGCAATTGAACCAgaagaa                                              <  2:348368/141‑1 (MQ=255)
                                                                                                tCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCCGCGCAATTGAACCAGAAGAAGGTTTATCCg                                    <  2:392438/141‑1 (MQ=255)
                                                                                                tCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCCGCGCAATTGAACCAGAAGAAGGTTTATCCg                                    <  1:327112/141‑1 (MQ=255)
                                                                                                   cGTCATTATATAAACGTGGTGGTACAATTTATACACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTCATCCGCGCAATTGAACCAGAAGAAGGTTTATCAGCTa                                 >  1:246696/1‑141 (MQ=255)
                                                                                                           tataAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCCGCGCAATTGAACCAGAAGAAGGTTTATCCGCTATGTTCCGt                         <  2:182835/141‑1 (MQ=255)
                                                                                                             taAACGTGGTGGTACAATTTATACACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCCGCGCAATTGCACCAGAAGAAGGTTTATCCGCTATGTTCCGTaa                       <  2:366687/141‑1 (MQ=255)
                                                                                                             taAACGTGGTGGTACAATTTATACACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAAACTGAAGGTATACCTGAAGGCGTACTTACCCGCGCAATTGAACCAGAAGAAGGTTTATCCGCTATGTCCCGTaa                       <  2:402189/141‑1 (MQ=255)
                                                                                                              aaaCGTGGTGGTACAATTTATACACATGTCATGCATGCGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCGGAAGGCGTACTTATCCGCGCAATTGAACCAGAAGAAGGTTTATCCGCTATGTTCCGTAAc                      <  2:249414/141‑1 (MQ=255)
                                                                                                                   tggtggTACAATTTATACACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCCGCGCAATTGAACCAGAAGAAGGTTTATCCGCTATGTTCCGTAACAGAgg                 <  2:104262/141‑1 (MQ=255)
                                                                                                                      tggtACAATTTATGCACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCCGCGCAATTGAACCAGAAGAAGGTTTATCCGCTATGTTCCGTAACAGAGGTaa              >  1:189194/1‑141 (MQ=255)
                                                                                                                      tggtACAATTTATACACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCCGCGCAATTGAACCAGAAGAAGGTTTATCCGCTATGTTCCGTAACAGAGGTaa              >  1:187165/1‑141 (MQ=255)
                                                                                                                      tggtACAATTTATACACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCCGCGCAATTGAACCAGAAGAAGGTTTATCCGCTATGTTCCGTAACAGAGGTaa              >  1:42237/1‑141 (MQ=255)
                                                                                                                                  tGCACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCCGCGCAATTGAACCAGAAGAAGGTTTATCCGCTATGTTCCGTAACAGAGGTAAGAAAGGCTACGa  <  1:44228/141‑1 (MQ=255)
                                                                                                                                  tGCACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCCGCGCAATTGAACCAGAAGAAGGTTTATCCGCTATGTTCCGTAACAGAGGTAAGAAAGGCTACGa  <  1:97525/141‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                           
CTCAAACGTATACAGGCTACATCGTGGAAACGGAAGCTTACTTAGGTTTGAATGATCGTGCGGCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCCGCGCAATTGAACCAGAAGAAGGTTTATCCGCTATGTTCCGTAACAGAGGTAAGAAAGGCTACGA  >  CP000730/2462627‑2462897

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