Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,440,251 G→A 7.8% G17G (GGC→GGT ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,440,2510GA 7.8% 208.6 / 6.5 77G17G (GGC→GGTompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base G (44/27);  new base A (2/4);  total (46/31)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGTTATTAACCCTCTGTTAT  >  minE/1440181‑1440320
                                                                      |                                                                     
ggCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTgcg                                                                       <  1:767531/71‑1 (MQ=255)
     tCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTg                                                                   <  1:3581500/70‑1 (MQ=255)
       tACTTTTCCGTACAAATCCAGTTTATTACCATCTTTATTATAGACTTCAGCCGCATTTGCTGCGCCAGc                                                                  <  1:3016962/69‑1 (MQ=2)
       tACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTa                                                                >  1:1528568/1‑71 (MQ=255)
        aCTTTTCCGTACAAATCCAGTTTATTACCATCTTTATTATAGACTTCAGCCGCATTTGCTGCGcc                                                                     >  1:2294653/1‑65 (MQ=2)
        aCTTTTCCGTACAAATCCAGTTTATTACCATCTTTATTATAGACTTCAGCCGCATTTGCTGCGcc                                                                     >  1:1285980/1‑65 (MQ=2)
        aCTTTTCCGTACAAATCCAGTTTATTACCATCTTTATTATAGACTTCAGCCGCATTTGCTGCGcc                                                                     >  1:470653/1‑65 (MQ=2)
        aCTTTTCCGTACAAATCCAGTTTATTACCATCTTTATTATAGACTTCAGCCGCATTTGCTGCGcc                                                                     >  1:165270/1‑65 (MQ=2)
         cTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTa                                                                <  1:4598928/69‑1 (MQ=255)
          tttACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTgcgc                                                                      <  1:934179/62‑1 (MQ=255)
          tttACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTgcgc                                                                      <  1:4836368/62‑1 (MQ=255)
          tttACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTAc                                                               >  1:4481721/1‑69 (MQ=255)
            tACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcag                                                            >  1:496206/1‑70 (MQ=255)
            tACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcag                                                            >  1:1621446/1‑70 (MQ=255)
            tACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcag                                                            >  1:4389517/1‑70 (MQ=255)
            tACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcag                                                            >  1:520484/1‑70 (MQ=255)
            tACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcag                                                            >  1:1937160/1‑70 (MQ=255)
              ccatACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGTGCAGCACCTGCGGCgagcag                                                         <  1:3618084/68‑1 (MQ=255)
                  aCAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcag                                                            >  1:2338551/1‑64 (MQ=255)
                  aCAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGTGCAGCACCTGcggc                                                               <  1:1619161/61‑4 (MQ=255)
                  aCAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGTGCAGCACCTGCGGCgagcag                                                         <  1:235588/67‑1 (MQ=255)
                   cAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTg                                                    <  1:1884941/71‑1 (MQ=255)
                    aGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTgg                                                   >  1:4972598/1‑71 (MQ=255)
                    aGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTgg                                                   >  1:1228954/1‑71 (MQ=255)
                        cTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGTGCAGCACCTGCGGCGAgcagggcaggaat                                                >  1:2648995/1‑67 (MQ=255)
                        cTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGTGCAGCACCTGCGGCGAgcagggcaggaat                                                >  1:4935129/1‑67 (MQ=255)
                         tAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGTGCAGCACCt                                                                    <  1:1763009/49‑1 (MQ=12)
                          aaTTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCag                                             >  1:2153372/1‑71 (MQ=255)
                           aTTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGAc                                                <  1:788333/67‑1 (MQ=255)
                            tttGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGAcc                                               <  1:5335710/67‑1 (MQ=255)
                            tttGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCagga                                           <  1:3937048/71‑1 (MQ=255)
                            tttGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCagga                                           <  1:204270/71‑1 (MQ=255)
                            tttGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCag                                             >  1:2103598/1‑69 (MQ=255)
                              tgttgCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCaggagg                                         >  1:1721441/1‑71 (MQ=255)
                               gttgCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCaggagg                                         <  1:3760441/70‑1 (MQ=255)
                               gttgCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCaggagg                                         <  1:988937/70‑1 (MQ=255)
                                   ccGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAg                                    <  1:75376/71‑1 (MQ=255)
                                     gTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTa                                  >  1:3093427/1‑71 (MQ=255)
                                      tCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGCGCTGGGACCAGGAGGGACAGTAc                                 >  1:2184546/1‑71 (MQ=255)
                                      tCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGa                                                 <  1:2971903/55‑1 (MQ=255)
                                      tCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGAcc                                               <  1:3814144/57‑1 (MQ=255)
                                      tCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTAc                                 >  1:1652277/1‑71 (MQ=255)
                                       cTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTa                                                                <  1:4306071/39‑1 (MQ=255)
                                        tttGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCag                                             <  1:4770638/57‑1 (MQ=255)
                                        tttGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACtt                               >  1:4708676/1‑71 (MQ=255)
                                        tttGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACtt                               >  1:3535897/1‑71 (MQ=255)
                                        tttGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACtt                               >  1:5492746/1‑71 (MQ=255)
                                        tttGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACtt                               >  1:2461263/1‑71 (MQ=255)
                                          tgttgtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGATCTGGGACCAGGAGGGACAGTACttt                              <  1:1410667/70‑1 (MQ=255)
                                          tgttgtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACttt                              >  1:4949594/1‑70 (MQ=255)
                                          tgttgtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACttt                              >  1:3792894/1‑70 (MQ=255)
                                          tgttgtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACttt                              <  1:169379/70‑1 (MQ=255)
                                          tgttgtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACttt                              >  1:3094092/1‑70 (MQ=255)
                                          tgttgtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACtt                               <  1:3888863/69‑1 (MQ=255)
                                             tgtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACca                                                             >  1:4931922/1‑36 (MQ=255)
                                             tgtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACca                                                             >  1:3527272/1‑36 (MQ=255)
                                             tgtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACca                                                             >  1:1156736/1‑36 (MQ=255)
                                              gtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTAc                                 >  1:4123593/1‑63 (MQ=255)
                                              gtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTAc                                 >  1:3194175/1‑63 (MQ=255)
                                              gtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACtt                         >  1:1963810/1‑71 (MQ=255)
                                               tAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACttt                        <  1:314447/71‑1 (MQ=255)
                                               tAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACttt                        >  1:3661286/1‑71 (MQ=255)
                                               tAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACtt                         <  1:649078/70‑1 (MQ=255)
                                                aaaCTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGa                                                 >  1:1004588/1‑45 (MQ=255)
                                                 aaCTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGc                                                      >  1:5261197/1‑39 (MQ=255)
                                                 aaCTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTc                       <  1:2271148/70‑1 (MQ=255)
                                                     tCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGt                   >  1:3377800/1‑70 (MQ=255)
                                                      cagcagCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGtta                 >  1:772416/1‑71 (MQ=255)
                                                      cagcagCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGt                   <  1:3587214/69‑1 (MQ=255)
                                                      cagcagCGTTTGATGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGtta                 >  1:3994858/1‑71 (MQ=255)
                                                       agcagcGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGttat                <  1:375657/71‑1 (MQ=255)
                                                       agcagcGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGt                   <  1:2394222/68‑1 (MQ=255)
                                                          agcGTTTGCTGCGACTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGTTATTaa             >  1:1496321/1‑71 (MQ=255)
                                                           gcGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGTTATTAAc            >  1:296204/1‑71 (MQ=255)
                                                           gcGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACATTAACTTTCATGTTATTAAc            >  1:115510/1‑71 (MQ=255)
                                                             gTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGTTATTAAccc          >  1:1424613/1‑71 (MQ=255)
                                                                tgctgcGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGTTATTAACCctct       >  1:372088/1‑71 (MQ=255)
                                                                  ctgcGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGTTATTAACCCTCTGt     <  1:1352748/71‑1 (MQ=255)
                                                                  ctgcGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGTTATTAACCCTCTGt     <  1:5150770/71‑1 (MQ=255)
                                                                  ctgcGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGTTATTAACCCTCTGt     <  1:1788677/71‑1 (MQ=255)
                                                                    gcgcCTGCTACCAGCAGAGCTGGGACCAGTAGGGACAGTACTTTAACTTTCATGTTATTAACCCTCTGtt    <  1:3564345/70‑1 (MQ=255)
                                                                    gcgcCTGCTACCAGCAGAGCTGGGACCAGTAGGGACAGTACTTTAACTTTCATGTTATTAACCCTCTGtt    <  1:3557689/70‑1 (MQ=255)
                                                                     cgcCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGTTATTAACCCTCTGTtat  <  1:5333586/71‑1 (MQ=255)
                                                                     cgcCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGTTATTAACCCTCTGTta   >  1:4495536/1‑70 (MQ=255)
                                                                     cgcCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGTTATTAACCCTCTGTta   >  1:512514/1‑70 (MQ=255)
                                                                     cgcCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGTTATTAACCCTCTGTta   <  1:90997/70‑1 (MQ=255)
                                                                     cgcCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGTTATTAACCCTCTGTta   >  1:4810845/1‑70 (MQ=255)
                                                                     cgcCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGTTATTAACCCTCTGTta   >  1:4489252/1‑70 (MQ=255)
                                                                      |                                                                     
GGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGTTATTAACCCTCTGTTAT  >  minE/1440181‑1440320

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: