Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,706,647 G→T 28.7% V66L (GTG→TTG)  tag → 3‑methyl‑adenine DNA glycosylase I, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,706,6470GT28.7% 35.3 / 22.2 35V66L (GTG→TTG) tag3‑methyl‑adenine DNA glycosylase I, constitutive
Reads supporting (aligned to +/- strand):  ref base G (14/11);  new base T (8/2);  total (22/13)
Fisher's exact test for biased strand distribution p-value = 2.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.27e-01

TGATAATGAGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTT  >  CP009273/3706502‑3706795
                                                                                                                                                 |                                                                                                                                                    
tGATAATGAGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGa                                                                                                                                                   >  2:658552/1‑149 (MQ=255)
tGATAATGAGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGATTATCGTGGATCACCGTCCTCAAAAAACGAGTGGCCTTTCGCGCCTGCTTTCATCAGTTCGATCCGTTGa                                                                                                                                                   <  2:425195/149‑1 (MQ=255)
           ggggCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGATTATCAAGGCTCCCCATCCACAAAAAACGCGAAAACTCTCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCgcagca                                                                                                                                        >  1:95258/1‑149 (MQ=255)
            gggCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCaa                                                                                                                                       >  2:104722/1‑149 (MQ=255)
             ggCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATAGCGCCTGCTTTCATCAGTTCGATCCGTGGAAGGTCGAAGCAAt                                                                                                                                      >  2:303712/1‑149 (MQ=255)
                       aaaCTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAagag                                                                                                                            <  2:609265/149‑1 (MQ=255)
                            gACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGt                                                                                                                       >  1:687464/1‑149 (MQ=255)
                                    aaaaCTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGCAAGAc                                                                                                               >  1:237624/1‑149 (MQ=255)
                                         tGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGGGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTa                                                                                                          >  1:834162/1‑149 (MQ=255)
                                         tGTTCGAAATGATCTGCCGTGAAGGGCAGCAGGCTGGATTATCGGGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTa                                                                                                          <  2:99176/149‑1 (MQ=255)
                                           ttCGAAATGATCTGCCTTGAAGGACAGCCGGCTGAATTATCGTGGATCACCGTACTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTCCa                                                                                                        >  2:479778/1‑149 (MQ=255)
                                               aaaTGATCTGCCTTGAAGGGCAGCAGGCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCGGTTCGATCCGGTGCAGGTCGCACCACTGCAGGGAAAGGATGTCGAAAGACTGGTACCGGAc                                                                                                    >  2:199777/1‑149 (MQ=255)
                                                      cTGCCTTGAAGGGCAGCAGGCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGTTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGa                                                                                             >  1:715061/1‑149 (MQ=255)
                                                      cTGCCTTGAAGGGCAGCAGGCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGa                                                                                             <  1:104722/149‑1 (MQ=255)
                                                          cTTGAAGGGCAGCAGGCTGGATTATCGTGGATCACCGTCCGCGAAAAACGAGAAAAGTATCGCGCCTGCTTTCATCAGTTCGATCAGGTGGAGGACGCAGCAATGCAGGAAGAGGGTGCCGAAAGACTGGTACAGGCGGCCGGGATTaa                                                                                         >  1:423117/1‑148 (MQ=255)
                                                               aGGGCAGCAGGCTGGATTATCGTGGATCAGCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGTTGAAGGACGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGCCGCCGGGCTTATCCGcc                                                                                    <  1:640650/149‑1 (MQ=255)
                                                                     gcagGCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGCAGAGGATGTCGAAAGACTGGTCCAGGACGCCCGGATTCTCCGCCATCTAg                                                                              >  2:799884/1‑149 (MQ=255)
                                                                               ttATCGTGGTTCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGGTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCa                                                                    <  1:374937/149‑1 (MQ=255)
                                                                                tagcGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTGTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCTATCGAGGGAAAATTCAg                                                                   <  1:510229/146‑1 (MQ=255)
                                                                                tATCGTGGATGACGCTCCGTAAAAAACGCGAAAACTATCGCGCCGGCGTTCCTCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGAGGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAg                                                                   <  1:509670/149‑1 (MQ=255)
                                                                                  tCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGTTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGc                                                                 >  1:660564/1‑149 (MQ=255)
                                                                                        aTCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGCGGGAAAATTCAGGCAattat                                                           >  1:374965/1‑149 (MQ=255)
                                                                                        aTCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAattat                                                           >  2:787893/1‑149 (MQ=255)
                                                                                                  tCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGTTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATgcg                                                 >  2:280260/1‑149 (MQ=255)
                                                                                                           gcgAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTAc                                        <  2:788081/149‑1 (MQ=255)
                                                                                                            cgAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGTTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGCTTATCCGCCATCGAGGGAAAATTCAGGCAATTCTTGGTAATGCGCGGGctcacc                                       >  1:72319/1‑144 (MQ=255)
                                                                                                             gAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGTTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCt                                      >  1:70660/1‑149 (MQ=255)
                                                                                                                      cgcgCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATGCAGGCAATTATTGGTAATGCGGGGGCGTACCTGCAAATGGa                             >  1:656139/1‑149 (MQ=255)
                                                                                                                       gcgcCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGaa                            <  2:687464/149‑1 (MQ=255)
                                                                                                                       gcgcCTGCTTTCATCAGGTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGaa                            <  1:436624/149‑1 (MQ=255)
                                                                                                                          ccTGCTTTCATCAGTTCGATCCGTTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAg                         >  2:731915/1‑149 (MQ=255)
                                                                                                                          ccTGCTTTCATCAGTTCGATCCGTTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATAATTGGTAATGCGCGGGCGTACCTCCAAATGGAACAg                         >  2:172654/1‑149 (MQ=255)
                                                                                                                                       gTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGt            <  2:75889/149‑1 (MQ=255)
                                                                                                                                         tCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGttt          <  2:320705/149‑1 (MQ=255)
                                                                                                                                         tCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGttt          >  2:662495/1‑149 (MQ=255)
                                                                                                                                                ggTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGGCt   >  2:375775/1‑149 (MQ=255)
                                                                                                                                                 gTGAAGGTCGCAGCAATGAAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGAGTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACtt  <  1:418072/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TGATAATGAGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTT  >  CP009273/3706502‑3706795

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