Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,187,308 T→C 34.6% P183P (CCA→CCG nanT ← sialic acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,187,3080TC34.6% 22.3 / 23.1 26P183P (CCA→CCGnanTsialic acid transporter
Reads supporting (aligned to +/- strand):  ref base T (14/3);  new base C (4/5);  total (18/8)
Fisher's exact test for biased strand distribution p-value = 7.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGAC  >  minE/2187244‑2187370
                                                                |                                                              
cGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCaaa                                                          <  1:466545/71‑1 (MQ=255)
              cAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCAAAATGcc                                                     >  1:93162/1‑62 (MQ=255)
                  cttccgcttccgGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCAAAATGc                                                      >  1:781662/1‑57 (MQ=255)
                     ccgcttccgGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCaaaa                                                         >  1:1096012/1‑51 (MQ=255)
                             gggATGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCAAAATGCCGATAAAGAACAGCGCa                                     >  1:1884972/1‑63 (MQ=255)
                                aTGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCaa                                                           <  1:2112608/38‑1 (MQ=255)
                                aTGTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGccc                            >  1:232810/1‑69 (MQ=255)
                                  gTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGa                        >  1:413820/1‑71 (MQ=255)
                                    ttttACGCAGCCAGAGAGCAAAGATGATTGGCaaaa                                                         >  1:55073/1‑36 (MQ=255)
                                     tttACGCAGCCAGAGAGCAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCg                     >  1:2089105/1‑71 (MQ=255)
                                             gCCAGAGAGCAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGc                                  <  1:142262/50‑1 (MQ=255)
                                                  agagCAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAaccac              <  1:2143354/65‑1 (MQ=255)
                                                  agagCAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAaccac              <  1:431798/65‑1 (MQ=255)
                                                   gagCAAAGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCta       <  1:666745/71‑1 (MQ=255)
                                                   gagCAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCta       <  1:221450/71‑1 (MQ=255)
                                                      cAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAg           <  1:590402/64‑1 (MQ=255)
                                                        aaGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAgcc         >  1:999982/1‑62 (MQ=255)
                                                        aaGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGc         >  1:592097/1‑64 (MQ=255)
                                                        aaGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGc         >  1:652922/1‑64 (MQ=255)
                                                        aaGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGc         >  1:136974/1‑64 (MQ=255)
                                                        aaGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGAc  >  1:2095945/1‑71 (MQ=255)
                                                          gatgatTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAac                 >  1:129231/1‑54 (MQ=255)
                                                            tgatTGGCAAAATGCCGATAAAGAACAGCGCACGc                                  >  1:1408791/1‑35 (MQ=255)
                                                             gatTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCa                          >  1:1616168/1‑42 (MQ=255)
                                                             gatTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCa                          >  1:1291602/1‑42 (MQ=255)
                                                              atTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCa                          >  1:543916/1‑41 (MQ=255)
                                                                |                                                              
CGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGAC  >  minE/2187244‑2187370

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