Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,440,244 T→G 6.9% N20H (AAC→CAC) ‡ ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,440,2440TG 6.9% 221.0 / 12.0 87N20H (AAC→CAC) ‡ompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base T (43/38);  new base G (4/2);  total (47/40)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

GTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGTTATTAACCCT  >  minE/1440176‑1440313
                                                                    |                                                                     
gTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTtg                                                                     <  1:3793346/71‑1 (MQ=255)
    aGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTtgctgc                                                                 <  1:3771150/71‑1 (MQ=255)
    aGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGtgcag                                                                  >  1:233014/1‑68 (MQ=12)
        cGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCt                                                             >  1:3326044/1‑71 (MQ=255)
        cGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCt                                                             <  1:1553432/71‑1 (MQ=255)
        cGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCt                                                             >  1:3762868/1‑71 (MQ=255)
        cGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCt                                                             >  1:593970/1‑71 (MQ=255)
          tCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTtgct                                                                   <  1:814031/63‑1 (MQ=255)
          tCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTtgct                                                                   <  1:322023/63‑1 (MQ=255)
          tCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGc                                                           <  1:4104165/71‑1 (MQ=255)
          acTACTTTTCCGTACAAATCCAGTTTATTACCATCTTTATTATAGACTTCAGCCGCATTTGCTGCGCCAGc                                                           <  1:2431350/70‑1 (MQ=2)
           caaCTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGTGCAGCACCTGcg                                                          >  1:3655588/3‑70 (MQ=12)
            tACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCt                                                          >  1:794748/1‑70 (MQ=255)
            tACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCt                                                          >  1:2471053/1‑70 (MQ=255)
            tACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTa                                                         >  1:3587972/1‑71 (MQ=255)
            tACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTa                                                         >  1:3331113/1‑71 (MQ=255)
            tACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCCCCTGCTa                                                         >  1:939324/1‑71 (MQ=255)
              cTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTtgc                                                                    >  1:980708/1‑58 (MQ=255)
              cTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTa                                                         <  1:3433959/69‑1 (MQ=255)
              cTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTAc                                                        >  1:1208978/1‑70 (MQ=255)
               tttACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGttt                                                                      >  1:1812468/1‑55 (MQ=255)
                ttACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTtgctgc                                                                 >  1:1136594/1‑59 (MQ=255)
                ttACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcag                                                     >  1:3663847/1‑71 (MQ=255)
                ttACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcag                                                     >  1:2607038/1‑71 (MQ=255)
                 tACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTAcc                                                       >  1:132367/1‑68 (MQ=255)
                 tACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcag                                                     >  1:858912/1‑70 (MQ=255)
                 tACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcag                                                     >  1:636137/1‑70 (MQ=255)
                 tACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcag                                                     >  1:3788287/1‑70 (MQ=255)
                   ccGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcagcag                                                  <  1:432744/71‑1 (MQ=255)
                   ccGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcagcag                                                  <  1:2371333/71‑1 (MQ=255)
                    cGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCaga                                                 <  1:1056624/71‑1 (MQ=255)
                    cGTACAAATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCaga                                                 <  1:3328512/71‑1 (MQ=255)
                     gTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCagag                                                <  1:2666157/71‑1 (MQ=255)
                     atACAGATCTAATTTGTTGCCGGCTTTATTATAAACTTCGGCTGCATGtgca                                                                   <  1:1315586/51‑2 (MQ=255)
                      tACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGTGCAGCACCTGCGGCgagcaggg                                                >  1:3501576/1‑68 (MQ=255)
                        cAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCt                                              <  1:1738901/70‑1 (MQ=255)
                         aGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcagca                                                   <  1:3897543/64‑1 (MQ=255)
                         aGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGTGCAGCACCTGCGGCGAgcagggcagg                                            <  1:709534/71‑1 (MQ=255)
                           aTCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGTGCAGCACCTGCGGCGAgcagggcaggaa                                          >  1:3958810/1‑69 (MQ=255)
                             cTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCagag                                                >  1:1837368/1‑63 (MQ=255)
                             cTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCagag                                                >  1:4043625/1‑63 (MQ=255)
                             cTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGAc                                         >  1:1722421/1‑70 (MQ=255)
                             cTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGAc                                         >  1:1211838/1‑70 (MQ=255)
                              tAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcagcag                                                  <  1:271597/60‑1 (MQ=255)
                               aaTTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTgg                                            >  1:891821/1‑65 (MQ=255)
                                aTTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTtgc                                                                    >  1:37051/1‑40 (MQ=255)
                                aTTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCagg                                     <  1:2844197/71‑1 (MQ=255)
                                aTTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCag                                      <  1:3001916/70‑1 (MQ=255)
                                aTTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCag                                      >  1:3850430/1‑70 (MQ=255)
                                aTTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCag                                      >  1:1635220/1‑70 (MQ=255)
                                aTTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCag                                      <  1:1220555/70‑1 (MQ=255)
                                 tttGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTtgctg                                                                  >  1:2487293/1‑41 (MQ=255)
                                 tttGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCagga                                    <  1:3297991/71‑1 (MQ=255)
                                 tttGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCag                                      >  1:3458204/1‑69 (MQ=255)
                                 tttGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCag                                      >  1:1221244/1‑69 (MQ=255)
                                   tgttgCCGTCTTTGTTGTAAACTTCAGCAGCGTTtgc                                                                    >  1:1305596/1‑37 (MQ=255)
                                   tgttgCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCag                                      >  1:2807663/1‑67 (MQ=255)
                                     ttgCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGa                                <  1:2438223/71‑1 (MQ=255)
                                        ccGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAg                             <  1:3334533/71‑1 (MQ=255)
                                           tCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCt                                                             <  1:1823416/36‑1 (MQ=255)
                                             tttGTTGTAAACTTGAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCagga                                    >  1:1808216/1‑59 (MQ=255)
                                             tttGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCaga                                                 >  1:1210149/1‑46 (MQ=255)
                                             tttGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCATAGCTGGGACCagga                                    >  1:590274/1‑59 (MQ=255)
                                             tttGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACtt                        >  1:260650/1‑71 (MQ=255)
                                               tgttgtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTa                           >  1:3604781/1‑66 (MQ=255)
                                               tgttgtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACttt                       <  1:2193015/70‑1 (MQ=255)
                                                gttgtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTa                      <  1:3488606/70‑1 (MQ=255)
                                                 ttgtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTaa                     <  1:1011212/70‑1 (MQ=255)
                                                  tgtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACt                   >  1:1116879/1‑71 (MQ=255)
                                                    tAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACt                         <  1:9456/63‑1 (MQ=255)
                                                    tAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACtt                  <  1:3962571/70‑1 (MQ=255)
                                                      aaCTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTc                <  1:3179039/70‑1 (MQ=255)
                                                      aaCTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTc                <  1:2702851/70‑1 (MQ=255)
                                                      aaCTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTc                <  1:2343046/70‑1 (MQ=255)
                                                      aaCTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCa               <  1:3863909/71‑1 (MQ=255)
                                                      aaCTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCa               >  1:683238/1‑71 (MQ=255)
                                                        cTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGa                                          >  1:2886290/1‑42 (MQ=255)
                                                           cagcagCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGtta          >  1:1469422/1‑71 (MQ=255)
                                                           cagcagCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGtta          >  1:2564176/1‑71 (MQ=255)
                                                           cagcagCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGtt           <  1:2795904/70‑1 (MQ=255)
                                                           cagcagCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACATTAACTTTCATGtt           <  1:19316/70‑1 (MQ=255)
                                                            agcagcGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACttt                 <  1:1659941/63‑1 (MQ=255)
                                                            agcagcGGTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGttat         <  1:4081949/71‑1 (MQ=255)
                                                             ggagcGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGttttt        <  1:3133954/69‑1 (MQ=255)
                                                             gcagcGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGttat         >  1:1691405/1‑70 (MQ=255)
                                                                 cGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGa                                <  1:989697/43‑1 (MQ=255)
                                                                  gTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGttatta       >  1:1514106/1‑67 (MQ=255)
                                                                  gTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGTTATTAAc     >  1:3586559/1‑69 (MQ=255)
                                                                   tttGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGTTATTAACCct  <  1:3210611/71‑1 (MQ=255)
                                                                   tttGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGTTATTAACCct  >  1:943163/1‑71 (MQ=255)
                                                                    |                                                                     
GTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGTTATTAACCCT  >  minE/1440176‑1440313

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