Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,440,236 A→G 7.6% A22A (GCT→GCC ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,440,2360AG 7.6% 217.7 / 7.7 79A22A (GCT→GCCompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base A (40/33);  new base G (4/2);  total (44/35)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGTTAT  >  minE/1440171‑1440306
                                                                 |                                                                      
aaaTAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCagcagc                                                                   >  1:3032291/1‑71 (MQ=255)
  aTAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGt                                                                 >  1:1576358/1‑71 (MQ=255)
     gTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTtg                                                              <  1:3793346/71‑1 (MQ=255)
         aGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTtgctgc                                                          <  1:3771150/71‑1 (MQ=255)
         aGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGtgcag                                                           >  1:233014/1‑68 (MQ=12)
             cGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCt                                                      >  1:3762868/1‑71 (MQ=255)
             cGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCt                                                      >  1:3326044/1‑71 (MQ=255)
             cGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCt                                                      >  1:593970/1‑71 (MQ=255)
             cGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCt                                                      <  1:1553432/71‑1 (MQ=255)
              gTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCg                                                                  <  1:3423355/58‑1 (MQ=255)
               tCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTtgct                                                            <  1:322023/63‑1 (MQ=255)
               tCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTtgct                                                            <  1:814031/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:3331113/1‑71 (MQ=255)
                 tACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTa                                                  >  1:3587972/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:3788287/1‑70 (MQ=255)
                      tACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcag                                              >  1:636137/1‑70 (MQ=255)
                      tACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcag                                              >  1:858912/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)
                           tACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCggctgca                                                                  <  1:2687171/45‑2 (MQ=25)
                           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:1635220/1‑70 (MQ=255)
                                     aTTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCag                               >  1:3850430/1‑70 (MQ=255)
                                     aTTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCag                               <  1:3001916/70‑1 (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:2343046/70‑1 (MQ=255)
                                                           aaCTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTc         <  1:3179039/70‑1 (MQ=255)
                                                           aaCTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTc         <  1:2702851/70‑1 (MQ=255)
                                                           aaCTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCa        >  1:683238/1‑71 (MQ=255)
                                                           aaCTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCa        <  1:3863909/71‑1 (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)
                                                                 |                                                                      
AAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGTTAT  >  minE/1440171‑1440306

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: