Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 138,172 A→G 13.5% N665N (AAT→AAC htrE ← predicted outer membrane usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE138,1720AG13.5% 84.7 / 8.1 37N665N (AAT→AAChtrEpredicted outer membrane usher protein
Reads supporting (aligned to +/- strand):  ref base A (15/17);  new base G (3/2);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

AAGTCAGTCCACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCATAA  >  minE/138102‑138226
                                                                      |                                                      
aaGTCAGTCCACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTa                                                        <  1:344376/71‑1 (MQ=255)
    cagtcCACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTg                                                    >  1:1122620/1‑71 (MQ=255)
         cACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGa                                                <  1:805729/70‑1 (MQ=255)
         cACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGa                                                <  1:1404032/70‑1 (MQ=255)
         cACCGCTATGCAATACAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGaa                                               <  1:2152791/70‑1 (MQ=255)
           ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAAt                                             <  1:109707/71‑1 (MQ=255)
            cGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATGGCTGTTTGCAGAAAtt                                            >  1:1159777/1‑71 (MQ=255)
                   cAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTg                                     >  1:1320777/1‑71 (MQ=255)
                      tACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTa                                                        <  1:128942/49‑1 (MQ=255)
                       aCAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCa                                  >  1:535032/1‑70 (MQ=255)
                       aCAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTAtt                                                      >  1:21342/1‑50 (MQ=255)
                       aCAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAg                                 >  1:860998/1‑71 (MQ=255)
                        cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGaa                                         >  1:503905/1‑62 (MQ=255)
                        cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTAtt                                                      >  1:214569/1‑49 (MQ=255)
                        cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTAtt                                                      >  1:1637024/1‑49 (MQ=255)
                         aaaaCCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGaa                                         >  1:771119/1‑61 (MQ=255)
                         aaaaCCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCg                               <  1:715777/71‑1 (MQ=255)
                          aaaCCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCg                               <  1:851341/70‑1 (MQ=255)
                           aaCCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGc                                    <  1:1984532/64‑1 (MQ=255)
                             ccaccGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAAc                                        >  1:589097/1‑58 (MQ=255)
                              caccGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTccc                          <  1:1313231/71‑1 (MQ=255)
                                         cTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCAt                       >  1:1629531/1‑63 (MQ=255)
                                         cTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCAt                       >  1:1614653/1‑63 (MQ=255)
                                         cTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACt                >  1:1259263/1‑70 (MQ=255)
                                          tGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCAt                       <  1:310656/62‑1 (MQ=255)
                                          tGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCAt                       <  1:489221/62‑1 (MQ=255)
                                            agagaAACTTGACGGCTGTTATCGCTATTTGCAGTAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATa            <  1:1588688/71‑1 (MQ=255)
                                            agagaAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCAt             <  1:1453082/70‑1 (MQ=255)
                                            agagaAACTTGACGGATGTTATCGCTATTTGCAGAAATTGAACCTGCCAGGGTTCCCCATGGTGACTCa              >  1:1209693/1‑69 (MQ=255)
                                            agagaAACTTACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATgg                     >  1:2130939/1‑61 (MQ=255)
                                               gaAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCt         <  1:1559575/71‑1 (MQ=255)
                                               gaAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCt         <  1:708396/71‑1 (MQ=255)
                                                aaaCTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGt                    <  1:16656/59‑1 (MQ=255)
                                                    ttGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGt                              >  1:649171/1‑45 (MQ=255)
                                                      gACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCAt    <  1:1858927/69‑1 (MQ=255)
                                                       aCGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCAAGCGTTCCCCATGGTGACTCATAGCTGGCATaa  <  1:911941/70‑1 (MQ=255)
                                                            tGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCAt    >  1:325755/1‑63 (MQ=255)
                                                              ttATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTccc                          <  1:653766/39‑1 (MQ=255)
                                                              ttATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCa              <  1:344289/51‑1 (MQ=255)
                                                                      aTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACt                <  1:587155/41‑1 (MQ=255)
                                                                      aTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACt                <  1:890348/41‑1 (MQ=255)
                                                                      |                                                      
AAGTCAGTCCACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCATAA  >  minE/138102‑138226

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: