Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,440,242 G→A 7.7% N20N (AAC→AAT) ‡ ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,440,2420GA 7.7% 198.3 / 5.0 65N20N (AAC→AAT) ‡ompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base G (35/25);  new base A (2/3);  total (37/28)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

AATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGTTATTAACCC  >  minE/1440172‑1440312
                                                                      |                                                                      
aaTAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCg                                                                        <  1:1566089/71‑1 (MQ=255)
  tAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGtt                                                                      >  1:2540822/1‑71 (MQ=255)
    gtccAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGttt                                                                     <  1:2142652/67‑1 (MQ=255)
    gTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGttt                                                                     <  1:1046119/70‑1 (MQ=255)
    gTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGttt                                                                     <  1:3622522/70‑1 (MQ=255)
    gTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGttt                                                                     <  1:1923652/70‑1 (MQ=255)
    gTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTtg                                                                    <  1:1767317/71‑1 (MQ=255)
        aGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGtgcag                                                                 <  1:563832/70‑3 (MQ=12)
         ggCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTtg                                                                    <  1:476000/66‑1 (MQ=255)
            cGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCt                                                            >  1:2511712/1‑71 (MQ=255)
              tCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGc                                                          <  1:3700968/71‑1 (MQ=255)
               cTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGc                                                          >  1:3220835/1‑70 (MQ=255)
                tACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTtgctgc                                                                >  1:3232651/1‑63 (MQ=255)
                tACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTtgctgc                                                                >  1:1141268/1‑63 (MQ=255)
                tACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTtgctgc                                                                >  1:2443296/1‑63 (MQ=255)
                  cTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTtg                                                                    >  1:3695843/1‑57 (MQ=255)
                   tttACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGc                                                          >  1:1568550/1‑66 (MQ=255)
                     tACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACca                                                     <  1:3271827/69‑1 (MQ=255)
                        cGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcagca                                                  >  1:31472/1‑69 (MQ=255)
                          tACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGc                                              <  1:3335980/71‑1 (MQ=255)
                          tACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGc                                              <  1:3365432/71‑1 (MQ=255)
                           aCAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGtgcagc                                                                <  1:135513/52‑1 (MQ=16)
                           aCAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGTGCAGCACCTGCGGCgagcaggg                                               >  1:2582446/1‑67 (MQ=255)
                            cAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCt                                             <  1:3209872/70‑1 (MQ=255)
                            cAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTg                                            <  1:1699741/71‑1 (MQ=255)
                            cAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTg                                            >  1:154628/1‑71 (MQ=255)
                             aGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCt                                             <  1:1511783/69‑1 (MQ=255)
                             aGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGTGCAGCACCTGCGGCGAgcagggcagg                                           <  1:735141/71‑1 (MQ=255)
                               aTCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGa                                         >  1:3899087/1‑71 (MQ=255)
                               aTCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGa                                         >  1:1703399/1‑71 (MQ=255)
                                tCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGAc                                        >  1:2645839/1‑71 (MQ=255)
                                tCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGAc                                        >  1:1112491/1‑71 (MQ=255)
                                 cTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGAc                                        >  1:3276002/1‑70 (MQ=255)
                                 cTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGAc                                        >  1:2993910/1‑70 (MQ=255)
                                   aaTTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCag                                     >  1:2720678/1‑71 (MQ=255)
                                    aTTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcagcag                                                 >  1:1453407/1‑58 (MQ=255)
                                    aTTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGa                                         <  1:711962/66‑1 (MQ=255)
                                    aTTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCagg                                    <  1:669712/71‑1 (MQ=255)
                                    aTTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCag                                     >  1:333375/1‑70 (MQ=255)
                                    aTTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCag                                     >  1:2489047/1‑70 (MQ=255)
                                    aTTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGTGCAGCACCTg                                                           >  1:3165087/1‑48 (MQ=12)
                                     tttGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCagga                                   <  1:2817526/71‑1 (MQ=255)
                                     tttGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCagga                                   <  1:3751257/71‑1 (MQ=255)
                                       tgttgCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCaggagg                                 >  1:3042709/1‑71 (MQ=255)
                                         ttgCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTa                                                        >  1:2374682/1‑46 (MQ=255)
                                           gCCGTCTTTGTTGTAAACTTCAGCAGCGTTtgctgc                                                                >  1:2959110/1‑36 (MQ=255)
                                             cGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTgg                                           <  1:3276464/55‑1 (MQ=255)
                                               tCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTAc                         >  1:420652/1‑71 (MQ=255)
                                                   tgttgtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACttt                      <  1:349458/70‑1 (MQ=255)
                                                   tgttgtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACttt                      >  1:1298770/1‑70 (MQ=255)
                                                   tgttgtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACt                        <  1:3597181/68‑1 (MQ=255)
                                                    gttgtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTaa                    >  1:2999125/1‑71 (MQ=255)
                                                    gttgtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTaa                    >  1:277056/1‑71 (MQ=255)
                                                      tgtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTaa                    <  1:2823685/69‑1 (MQ=255)
                                                      tgtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACt                  >  1:881125/1‑71 (MQ=255)
                                                        tAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACt                  >  1:1787957/1‑69 (MQ=255)
                                                          aaCTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCa              <  1:228369/71‑1 (MQ=255)
                                                          aaCTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCa              >  1:1546298/1‑71 (MQ=255)
                                                          aaCTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCa              >  1:964902/1‑71 (MQ=255)
                                                           aCTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCAt             >  1:1114062/1‑71 (MQ=255)
                                                           aCTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCAt             >  1:97950/1‑71 (MQ=255)
                                                            cTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATg            >  1:1829322/1‑71 (MQ=255)
                                                              tCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGt           >  1:901488/1‑70 (MQ=255)
                                                                  cagcGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTaa                    <  1:3947422/57‑1 (MQ=255)
                                                                  cagcGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTaa                    <  1:2920811/57‑1 (MQ=255)
                                                                     cGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGTTATTAAcc   <  1:1786993/71‑1 (MQ=255)
                                                                      gTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGTTATTAAccc  <  1:998020/71‑1 (MQ=255)
                                                                      |                                                                      
AATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGTTATTAACCC  >  minE/1440172‑1440312

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: