Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 261,102 C→T 33.2% P73P (CCG→CCA yajI ← predicted lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE261,1020CT33.2% 50.2 / 31.0 63P73P (CCG→CCAyajIpredicted lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (31/11);  new base T (6/15);  total (37/26)
Fisher's exact test for biased strand distribution p-value = 9.65e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.62e-01

CCATCTGCATCAGGCGTAATATTCACCAGCGACATACGTAAAGTACCGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCAGGCAGAAGGTAAACCCCACTGCTGGATTTTGCATTCAGC  >  minE/261034‑261156
                                                                    |                                                      
ccATCTGCATCAGGCGTAATATTCACCAGCGACATACGTAAAGTACCGATCTGGCTTTCCAGTCTTGcc                                                        >  1:852128/1‑69 (MQ=255)
     tGCATCAGGCGTAATATTCACCAGCGACATACGTAAATTACCGATCTGGCTTTCCAGTCTTGCCGgtgt                                                   >  1:2742813/1‑69 (MQ=255)
     tGCATCAGGCGTAATATTCACCAGCGACATACGTAAAGTACCGATCTGGCTTTCCAGTCTTGCCGgtgt                                                   >  1:2673310/1‑69 (MQ=255)
           aGGCGTAATATTCACCAGCGACATACGTAAAGTACCGATCTGGCTTTCCAGTCTTGCTGGTGTTTTTGc                                             >  1:15710/1‑69 (MQ=255)
           aGGCGTAATATTCACCAGCGACATACGTAAAGTACCGATCTGGCTTTCCAGTCTTGCTGGTGTTTTTGc                                             >  1:1940665/1‑69 (MQ=255)
           aGGCGTAATATTCACCAGCGACATACGTAAAGTACCGATCTGGCTTTCCAGTCTTGCTGGTGTTTTTGc                                             >  1:152231/1‑69 (MQ=255)
           aGGCGTAATATTCACCAGCGACATACGTAAAGTACCGATCTGGCTTTCCAGTCTTGCTGGTGTTTTTGc                                             >  1:399644/1‑69 (MQ=255)
           aGGCGTAATATTCACCAGCGACATACGTAAAGTACCGATCTGGCTTTCCAGTCTTGCTGGTGTTTTTGc                                             >  1:523321/1‑69 (MQ=255)
           aGGCGTAATATTCACCAGCGACATACGTAAAGTACCGATCTGGCTATCCAGTCTTGCTGGTGTTTTTGc                                             >  1:3379486/1‑69 (MQ=255)
                  atatTCACCAGCGACATACGTAAAGTACCGATCTGGCTTTCCAGTCTTGCCGGTGTTTTc                                               >  1:2354718/1‑60 (MQ=255)
                  atatTCACCAGCGACATACGTAAAGTACCGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCa                                         >  1:347939/1‑66 (MQ=255)
                  atatTCACCAGCGACATACGTAAAGTACCGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCa                                         >  1:74392/1‑66 (MQ=255)
                  atatTCACCAGCGACATACGTAAAGTACCGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCa                                         >  1:997930/1‑66 (MQ=255)
                     ttCACCAGCGACATACGTAAAGTACCGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCAGGCAGa                                   >  1:1541961/1‑69 (MQ=255)
                     ttCACCAGCGACATACGTAAAGTACCGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCAGGCAGa                                   >  1:2154569/1‑69 (MQ=255)
                          cAGCGACATACGTAAAGTACCGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCAGGCAGAAGGTa                              <  1:2832718/69‑1 (MQ=255)
                          cAGCGACATACGTAAAGTACCGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCAGGCAGAAGGTa                              <  1:1746527/69‑1 (MQ=255)
                             cGACATACGTAAAGTACCGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCAGGCAGAAGGTAAAc                           <  1:3229922/69‑1 (MQ=255)
                                cATACGTAAAGTACCGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCAGGCAGAAGGTAAAcccc                        <  1:1232966/69‑1 (MQ=255)
                                cATACGTAAAGTACCGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCAGGCAGAAGGTAAAcccc                        <  1:187004/69‑1 (MQ=255)
                                    cGTAAAGTACCGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCAGGCAGAAGGTAAAcccc                        >  1:539520/1‑65 (MQ=255)
                                    cGTAAAGTACCGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCAGGCAGAAGGTAAAcccc                        >  1:1382605/1‑65 (MQ=255)
                                    cGTAAAGTACCGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCAGGCAGAAGGTAAAcccc                        >  1:2393135/1‑65 (MQ=255)
                                    cGTAAAGTACCGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCAGGCAGAAGGTAAAcccc                        >  1:2794669/1‑65 (MQ=255)
                                     gTAAAGTACCGATCTGGCTTTCCAGTCTTGCTGGTGtttttg                                              <  1:2017719/42‑3 (MQ=38)
                                     gTAAAGTACCGATCTGGCTTTCCAGTCTTGCTGGTGtttttg                                              <  1:763606/42‑3 (MQ=38)
                                     gTAAAGTACCGATCTGGCTTTCCAGTCTTGCTGGTGtttttg                                              <  1:40480/42‑3 (MQ=38)
                                     gTAAAGTACCGATCTGGCTTTCCAGTCTTGCTGGTGtttttg                                              <  1:1433129/42‑3 (MQ=38)
                                     gTAAAGTACCGATCTGGCTTTCCAGTCTTGCTGGTGtttttg                                              <  1:3085049/42‑3 (MQ=38)
                                     gTAAAGTACCGATCTGGCTTTCCAGTCTTGCTGGTGtttttg                                              <  1:2980926/42‑3 (MQ=38)
                                     gTAAAGTACCGATCTGGCTTTCCAGTCTTGCTGGTGtttttg                                              <  1:2941386/42‑3 (MQ=38)
                                     gTAAAGTACCGATCTGGCTTTCCAGTCTTGCTGGTGtttttg                                              <  1:2883900/42‑3 (MQ=38)
                                     gTAAAGTACCGATCTGGCTTTCCAGTCTTGCTGGTGtttttg                                              <  1:1595177/42‑3 (MQ=38)
                                     gTAAAGTACCGATCTGGCTTTCCAGTCTTGCTGGTGtttttg                                              <  1:1602836/42‑3 (MQ=38)
                                     gTAAAGTACCGATCTGGCTTTCCAGTCTTGCTGGTGtttttg                                              <  1:161744/42‑3 (MQ=38)
                                     gTAAAGTACCGATCTGGCTTTCCAGTCTTGCTGGTGtttttg                                              <  1:2490795/42‑3 (MQ=38)
                                     gTAAAGTACCGATCTGGCTTTCCAGTCTTGCTGGTGtttttg                                              <  1:2406054/42‑3 (MQ=38)
                                     gTAAAGTACCGATCTGGCTTTCCAGTCTTGCTGGTGtttttg                                              <  1:2139525/42‑3 (MQ=38)
                                     gTAAAGTACCGATCTGGCTTTCCAGTCTTGCTGGTGtttttg                                              <  1:226423/42‑3 (MQ=38)
                                        aaGTACCGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCAGGCAGAAGGTa                              >  1:576485/1‑55 (MQ=255)
                                        aaGTACCGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCAGGCAGAAGGTa                              >  1:2035505/1‑55 (MQ=255)
                                        aaGTACCGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCAGGCAGAAGGTa                              >  1:491026/1‑55 (MQ=255)
                                        aaGTACCGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCAGGCAGAAGGTa                              >  1:1275639/1‑55 (MQ=255)
                                        aaGTACCGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCAGGCAGAAGGTa                              >  1:86279/1‑55 (MQ=255)
                                        aaGTACCGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCAGGCAGAAGGTa                              >  1:2304652/1‑55 (MQ=255)
                                        aaGTACCGATCTGGCTTTCCACTCTTGCCGGTGTTTTCGCTCCAGGCAGAAg                                 >  1:2470358/1‑52 (MQ=255)
                                             ccGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCAGGCAGAAGGTAAACCCCACTGCTGGAttt            >  1:1592062/1‑68 (MQ=255)
                                             ccGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCAGGCAGAAGGTAAACCCCACTGCTGGAttt            >  1:979630/1‑68 (MQ=255)
                                             ccGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCAGGCAGAAGGTAAACCCCACTGCTGGAttt            >  1:1275057/1‑68 (MQ=255)
                                             ccGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCAGGCAGAAGGTAAACCCCACTGCTGGAttt            >  1:1972407/1‑68 (MQ=255)
                                             ccGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCAGGCAGAAGGTAAACCCCACTGCTGGAttt            >  1:150985/1‑68 (MQ=255)
                                             ccGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCAGGCAGAAGGTAAACCCCACTGCTGGAttt            >  1:447110/1‑68 (MQ=255)
                                             ccGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCAGGCAGAAGGTAAACCCCACTGCTGGAttt            >  1:2089607/1‑68 (MQ=255)
                                             ccGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCAGGCAGAAGGTAAACCCCACTGCTGGAttt            >  1:1087553/1‑68 (MQ=255)
                                             ccGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCAGGCAGAAGGTAAACCCCACTGCTGGAtt             >  1:996878/1‑67 (MQ=255)
                                                  cTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCAGGCAGAAGGTAAACCCCACTGCTGGATTTTGCa        >  1:2750622/1‑67 (MQ=255)
                                                               cTTGCCGGTGTTTTCGCTCCAGGCAGAAGGTAAACCCCACTGCTGGATTTTGCATTc     >  1:3230968/1‑57 (MQ=255)
                                                               cTTGCCGGTGTTTTCGCTCCAGGCAGAAGGTAAACCCCACTGCTGGATTTTGCATTc     >  1:2890210/1‑57 (MQ=255)
                                                                  gCCGGTGTTTTCGCTCCAGGCAGAAGGTAAACCCCACTGCTGGATTTTGCATTCAGc  <  1:1888243/57‑1 (MQ=255)
                                                                  gCCGGTGTTTTCGCTCCAGGCAGAAGGTAAACCCCACTGCTGGATTTTGCATTCAGc  <  1:1575087/57‑1 (MQ=255)
                                                                  gCCGGTGTTTTCGCTCCAGGCAGAAGGTAAACCCCACTGCTGGATTTTGCATTCAGc  <  1:3259284/57‑1 (MQ=255)
                                                                  gCCGGTGTTTTCGCTCCAGGCAGAAGGTAAACCCCACTGCTGGATTTTGCATTCAGc  <  1:862038/57‑1 (MQ=255)
                                                                  gCCGGTGTTTTCGCTCCAGGCAGAAGGTAAACCCCACTGCTGGATTTTGCATTCAGc  <  1:2260548/57‑1 (MQ=255)
                                                                  gCCGGTGTTTTCGCTCCAGGCAGAAGGTAAACCCCACTGCTGGATTTTGCATTCAGc  <  1:2785024/57‑1 (MQ=255)
                                                                    |                                                      
CCATCTGCATCAGGCGTAATATTCACCAGCGACATACGTAAAGTACCGATCTGGCTTTCCAGTCTTGCCGGTGTTTTCGCTCCAGGCAGAAGGTAAACCCCACTGCTGGATTTTGCATTCAGC  >  minE/261034‑261156

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