Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,151,211 A→G 22.2% V2V (GTA→GTG ydjQ → endonuclease of nucleotide excision repair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,151,2110AG22.2% 68.0 / 18.0 36V2V (GTA→GTGydjQendonuclease of nucleotide excision repair
Reads supporting (aligned to +/- strand):  ref base A (19/9);  new base G (5/3);  total (24/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.65e-01

TCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTAC  >  minE/1151146‑1151278
                                                                 |                                                                   
tctTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCACTCGGAGTCAACAGTGGTACGGCg                                                                >  1:188421/1‑71 (MQ=255)
  ttttttAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGtt                                                              >  1:1819401/1‑71 (MQ=255)
  ttttttAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGtt                                                              >  1:1507678/1‑71 (MQ=255)
            gCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTGCGGCGTTTAACTTCTcc                                                    >  1:947152/1‑71 (MQ=255)
                cccGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGttt                                                             <  1:1273222/58‑1 (MQ=255)
                 ccGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGc                                               <  1:1601948/71‑1 (MQ=255)
                   gCCTGTTACACTGGATATATAACCAGCATTCGGAGTCAACAGTGGTACGGCGttt                                                             <  1:1253363/55‑1 (MQ=255)
                   gCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTc                                             <  1:1704000/71‑1 (MQ=255)
                                gatagataACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTaa                                                           >  1:2127191/1‑44 (MQ=255)
                                gatagataACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTaa                                                           >  1:140605/1‑44 (MQ=255)
                                    gataACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAAttt                                       <  1:1524862/60‑1 (MQ=255)
                                       aaCCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAAt                         >  1:1038050/1‑71 (MQ=255)
                                          cAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTa                      <  1:2005187/71‑1 (MQ=255)
                                           aGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTAt                     <  1:1666316/71‑1 (MQ=255)
                                            gCATTCGGAGTCAACAGTGGTGCGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATg                    <  1:334848/71‑1 (MQ=255)
                                            gCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCcg                                                   >  1:414621/1‑40 (MQ=255)
                                            gCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATg                    <  1:1087255/71‑1 (MQ=255)
                                               ttCGGAGTCAACAGTGGTGCGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGaa                  <  1:1200654/70‑1 (MQ=255)
                                               ttCGGAGTCAACAGTGGTGCGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAat                 <  1:2071010/71‑1 (MQ=255)
                                                 cGGAGTCAACAGTGGTACGGCGTTTAAGTTCTCCGCGGCTc                                             <  1:1176109/41‑1 (MQ=255)
                                                   gAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATcc             >  1:1939880/1‑71 (MQ=255)
                                                        aaCAGTGGTGCGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGaa         >  1:2022286/1‑70 (MQ=255)
                                                        aaCAGTGGTGCGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGaa         >  1:1288689/1‑70 (MQ=255)
                                                        aaCAGTGGTGCGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGaa         >  1:542818/1‑70 (MQ=255)
                                                        aaCAGTGGTGCGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGaa         >  1:213374/1‑70 (MQ=255)
                                                        aaCAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCg           >  1:920696/1‑68 (MQ=255)
                                                        aaCAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGa          >  1:1915679/1‑69 (MQ=255)
                                                        aaCAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGa          >  1:1781616/1‑69 (MQ=255)
                                                         aCAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACa       >  1:2160080/1‑71 (MQ=255)
                                                         aCAGTGGTAAGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACa       >  1:13468/1‑71 (MQ=255)
                                                          cAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACAt      >  1:1690272/1‑71 (MQ=255)
                                                           aGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACAtt     >  1:2184491/1‑71 (MQ=255)
                                                            gTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCGGAACAtt     >  1:2132920/1‑70 (MQ=255)
                                                            gTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACAtt     >  1:1818369/1‑70 (MQ=255)
                                                              ggTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTAc  >  1:1987292/1‑71 (MQ=255)
                                                              ggTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTAc  >  1:192268/1‑71 (MQ=255)
                                                                 |                                                                   
TCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTAC  >  minE/1151146‑1151278

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