Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 52,517 T→C 23.5% E757E (GAA→GAG imp ← exported protein required for envelope biosynthesis and integrity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE52,5170TC23.5% 92.2 / 29.7 51E757E (GAA→GAGimpexported protein required for envelope biosynthesis and integrity
Reads supporting (aligned to +/- strand):  ref base T (19/20);  new base C (7/5);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 7.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

GATACGGCAGAATGTTCGAACGCAGCATCTCTTGCGTACCCAGACCGTAGTTGGAGCTCAGGCCGCGAAGTTCGATGTTAAAGCCGATTGCGTTGTCATATACCGCATGTTGTTTATCGTTATCCCAACCGTT  >  minE/52447‑52579
                                                                      |                                                              
gATACGGCAGAATGTTCGAACGCAGCATCTCTTGCGTACCCAGACCGTAGTTGGAGCTCAGGCCGCGAAGt                                                                >  1:990859/1‑71 (MQ=255)
  tACGGCAGAATGTTCGAACGCAGCATCTCTTGCGTACCCAGACCGTAGTTGGAGCTCAGGCCGCGAAGTTc                                                              <  1:1654464/71‑1 (MQ=255)
    cGGCAGAATGTTCGAACGCAGCATCTCTTGCGTACCCAGACCGTAGTTGGAGCTCAGGCCGCGAAGTTCGa                                                            <  1:249028/71‑1 (MQ=255)
         gAATGTTCGAACGCAGCATCTCTTGCGTACCCAGACCGTAGTTGGAGCTCAGGCCGCGAAGCTCGATGTTa                                                       <  1:1725648/71‑1 (MQ=255)
         gAATGTTCGAACGCAGCATCTCTTGCGTACCCAGACCGTAGTTGGAGCTCAGGACGCGAAGTTCGATGTTa                                                       <  1:2003137/71‑1 (MQ=255)
           aTGTTCGAACGCAGCATCTCTTGCGTACCCAGACCGTAGTTGGAGCTCAGGCCGCGAAGTTCGAt                                                           >  1:1011762/1‑65 (MQ=255)
           aTGTTCGAACGCAGCATCTCTTCGTACCCAGACCGTAGTTGGAGCTCAGGCCGCGAAGTTCGATGTTAAAg                                                    >  1:2039411/1‑71 (MQ=255)
               tCGAACGCAGCATCTCTTGCGTACCCAGACCGTAGTTGGAGCTCAGGCCGCGAAGTTCGATGTTAAAGcc                                                  >  1:1681455/1‑70 (MQ=255)
                cGAACGCAGCATCTCTTGCGTACCCAGACCGTAGTTGGAGCTCAGGCCGCGAAGTTCGATGTTAAAGCCGa                                                >  1:2026833/1‑71 (MQ=255)
                      cagcaTCTCTTGCGTACCCAGACCGTAGTTGGAGCTCAGGCCGCGAAGTTCGATGTTAAAGCCGATTGCGt                                          >  1:77385/1‑71 (MQ=255)
                      cagcaTCTCTTGCGTACCCAGACCGTAGTTGGAGCTCAGGCCGCGAAGCTCGATGTTAAAg                                                    >  1:2073931/1‑61 (MQ=255)
                            ctctTGCGTACCCAGACCGTAGTTGGAGCTCAGGCCGCGAAGCTCGATGTTAAAGCCGATTGCGTTGTCa                                     >  1:78091/1‑70 (MQ=255)
                            ctctTGCGTACCCAGACCGTAGTTGGAGCTCAGGCCGCGAAGCTCGATGTTAAAGCCGATTGCGTTGTCa                                     >  1:2148871/1‑70 (MQ=255)
                             tctTGCGTACCCAGACCGTAGTTGGAGCTCAGGCCGCGAAGTTCGATGTTAAAGCCGATTGCGTTGTCata                                   >  1:1432880/1‑71 (MQ=255)
                             tctTGCGTACCCAGACCGTAGTTGGAGCTCAGGCCGCGAAGTTCGATGTTAAAGCCGATTGCGTTGTCata                                   >  1:1593832/1‑71 (MQ=255)
                             tctTGCGTACCCAGACCGTAGTTGGAGCTCAGGCCGCGAAGTTCGATGTTAAAGCCGATTGCGTTGTCata                                   >  1:2150314/1‑71 (MQ=255)
                             tctTGCGTACCCAGACCGTAGTTGGAGCTCAGGCCGCGAAGTTCGATGTTAAAGCCGATTGCGTTGTCat                                    >  1:708360/1‑70 (MQ=255)
                               ttGCGTACCCCGACCGTAGTTGGAGCTCAGGCCGCGAAGCTCGATGTTAAAGCCGATTGCGTTGTCa                                     <  1:2182266/67‑1 (MQ=255)
                               ttGCGTACCCAGACCGTAGTTGGAGCTCAGGCCGCGAAGTTCGATGTTAAAGCCGATTGCGTTGTCa                                     <  1:1220242/67‑1 (MQ=255)
                               ttGCGTACCCAGACCGTAGTTGGAGCTCAGGCCGCGAAGTTCGATGTTAAAGCCGATTGCGTTGTCa                                     <  1:1689929/67‑1 (MQ=255)
                               ttGCGTACCCAGACCGTAGTTGGAGCTCAGGCCGCGAAGTTCGATGTTAAAGCCGATTGCGTTGTCa                                     <  1:1603001/67‑1 (MQ=255)
                               ttGCGTACCCAGACCGTAGTTGGAGCTCAGGCCGCGAAGTTCGATGTTAAAGCCGATTGCGTTGTCAtata                                 >  1:1111864/1‑71 (MQ=255)
                               ttGCGTACCCAGACCGTAGTTGGAGCTCAGGCCGCGAAGTGCGATGTTAAAGCCGATTGCGTTGTCa                                     <  1:1398260/67‑1 (MQ=255)
                               ttGCGTACCCAGACCGTAGTTGGAGCTCAGGCCGCGAAGCTCGATGTTAAAGCCGATTGCGTTGTCa                                     <  1:1962529/67‑1 (MQ=255)
                               ttGCGTACCCAGACCGTAGTTGGAGCTCAGGCCGCGAAGCTCGATGTTAAAGCCGATTGCGTTGTCa                                     <  1:1375538/67‑1 (MQ=255)
                                 gCGTACCCAGACCTTAGTTGGAGCTCAGGCCGCGAAGTTCGATGTTAAAGCCGATTGCGTTGTCatat                                  <  1:1174340/68‑1 (MQ=255)
                                 gCGTACCCAGACCGTAGTTGGAGCTCAGGCCGCGAAGTTCGATGTTAAAGCCGATTGCGTTGTCATATAcc                               >  1:1358985/1‑71 (MQ=255)
                                 gCGTACCCAGACCGTAGTTGGAGCTCAGGCCGCGAAGCTCGATGTTAAAGCCGATTGCGTTGTCATATAcc                               >  1:1083912/1‑71 (MQ=255)
                                 gCGTACCCAGACCGTAGTTGGAGCTCAGGCCGCGAAGCTCGATGTTAAAGCCGATTGCGTTGTCATATAcc                               >  1:1922944/1‑71 (MQ=255)
                                      cccAGACCGTAGTTGGAGCTCAGGCCGCGAAGTTCGATGTTAAAGCCGATTGCGTTGTCATATACCGCAt                           <  1:35221/70‑1 (MQ=255)
                                        cAGACCGTAGTTGGAGCTCAGGCCGCGAAGTTCGATGTTAAAGcc                                                  <  1:656561/45‑1 (MQ=255)
                                        cAGACCGTAGTTGGAGCTCAGGCCGCGAAGTTCGATGTTAAAGcc                                                  <  1:2121760/45‑1 (MQ=255)
                                        cAGACCGTAGTTGGAGCTCAGGCCGCGAAGTTCGATGTTAAAGCCGATTGCGTTGTCAtata                                 >  1:1009020/1‑62 (MQ=255)
                                        cAGACCGTAGTTGGAGCTCAGGCCGCGAAGCTCGATGTTAAAGCCGATTGCGTTGTCATATACCGCAtgtt                        >  1:473129/1‑71 (MQ=255)
                                         aGACCGTAGTTGGAGCTCAGGCCGCGAAGTTCGATGTTAAAGCCGATTGCGTTGTCATATACCGCAtgttg                       >  1:2013794/1‑71 (MQ=255)
                                            ccGTAGTTGGAGCTCAGGCCGCGAAGCTCGATGTTAAAGCCGATTGCGTTGTCATATACCGCATgttgtt                     >  1:1369464/1‑70 (MQ=255)
                                             cGTAGTTGGAGCTCAGGCCGCGAAGTTCGATGTTAAAGCCGATTGCGTTGTCat                                    <  1:1609521/54‑1 (MQ=255)
                                             cGTAGTTGGAGCTCAGGCCGCGAAGTTCGATGTTAAAGCCGATTGCGTTGTCATATACCGCATGTTGTTTa                   <  1:2212623/71‑1 (MQ=255)
                                             cGTAGTTGGAGCTCAGGCCGCGAAGCTCGATGTTAAAGCCGATTGCGTTGTCATATACCGCATGTTGTTTa                   <  1:243400/71‑1 (MQ=255)
                                                 gTTGGAGCTCAGGCCGCGAAGTTCGATGTTAAAGCCGATTGc                                            <  1:951658/42‑1 (MQ=255)
                                                 gTTGGAGCTCAGGCCGCGAAGTTCGATGTTAAAGCCGATTGCGTTGTCATATACCGCATGTTGTTTATCGt               <  1:1118631/71‑1 (MQ=255)
                                                 gTTGGAGCTCAGGCCGCGAAGTTCGATGTTAAAGCCGATTGCGTTGTCATATACCGCATGTTGTTTATCGt               <  1:2099909/71‑1 (MQ=255)
                                                    ggAGCTCAGGCCGCGAAGTTCGATGTTAAAGCCGATTGCGTTGTCATATACCGCATGTTGTTTATCGTTAt            >  1:581260/1‑71 (MQ=255)
                                                    ggAGCTCAGGCCGCGAAGTTCGATGTTAAAGCCGATTGCGTTGTCATATACCGCATGTTGTTTATCGTTAt            >  1:1910183/1‑71 (MQ=255)
                                                    ggAGCTCAGGCCGCGAAGTTCGATGTTAAAGCCGATTGCGTTGTCATATAC‑‑CATGTTGTTTATCGTTATc           <  1:1315466/70‑1 (MQ=255)
                                                     gAGCTCAGGCCGCGAAGTTCGATGTTAAAGCCGATTGCGTTGTCATATACCGCATGTTGTTTATCGTTAt            <  1:819143/70‑1 (MQ=255)
                                                             gCCGCGAAGTTCGATGTTAAAGCCGATTGCGTTGTCATATACCGCAtgttgt                      <  1:1396922/52‑1 (MQ=255)
                                                             gCCGCGAAGTTCGATGTTAAAGCCGATTGCGTTGTCATATACCGCAtgt                         <  1:2213139/49‑1 (MQ=255)
                                                              ccGCGAAGTTCGATGTTAAAGCCGATTGCGTTGTCATATACCGCATGTTGTTTATCGTTATCCCAACCGtt  >  1:1870559/1‑71 (MQ=255)
                                                              ccGCGAAGTTCGATGTTAAAGCCGATTGCGTTGTCATATACCGCATGTTGTTTATCGTTATCCCAACCGtt  >  1:1504248/1‑71 (MQ=255)
                                                              ccGCGAAGTTCGATGTTAAAGCCGATTGCGTTGTCATATACCGCATGTTGTTTATCGTTATCCCAACCGtt  >  1:1896067/1‑71 (MQ=255)
                                                               cgcgAAGTTCGATGTTAAAGCCGATTGCGTTGTCATATACCGCATGTTGTTTATCGTTATCCCAACCGtt  >  1:226993/1‑70 (MQ=255)
                                                                      |                                                              
GATACGGCAGAATGTTCGAACGCAGCATCTCTTGCGTACCCAGACCGTAGTTGGAGCTCAGGCCGCGAAGTTCGATGTTAAAGCCGATTGCGTTGTCATATACCGCATGTTGTTTATCGTTATCCCAACCGTT  >  minE/52447‑52579

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