Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,440,197 G→A 13.2% Y35Y (TAC→TAT ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,440,1970GA13.2% 234.7 / 21.0 91Y35Y (TAC→TATompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base G (45/34);  new base A (10/2);  total (55/36)
Fisher's exact test for biased strand distribution p-value = 1.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

GACGCATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAG  >  Exported/1440127‑1440262
                                                                      |                                                                 
gACGCATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCg                                                                   >  1:4437036/1‑71 (MQ=255)
 aCGCATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCg                                                                   >  1:4789128/1‑70 (MQ=255)
  cGCATGTATGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTa                                                                 <  1:4228482/71‑1 (MQ=255)
    cATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACa                                                               <  1:5384074/71‑1 (MQ=255)
      tGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCg                                                                   <  1:4823278/65‑1 (MQ=255)
      tGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGa                                                             <  1:1611298/71‑1 (MQ=255)
        tAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATc                                                           >  1:371044/1‑71 (MQ=255)
        tAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATc                                                           >  1:757451/1‑71 (MQ=255)
            tCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAAt                                                       <  1:106870/71‑1 (MQ=255)
             cTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGTCCGTCTACTTTACCGTACAGATCTAAtt                                                      >  1:46530/1‑71 (MQ=255)
                gTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCGGTCTACTTTACCGTACAGATCTAAttt                                                     <  1:2837429/69‑1 (MQ=255)
                gTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAAttt                                                     <  1:5245665/69‑1 (MQ=255)
                gTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATttgt                                                   <  1:4054538/71‑1 (MQ=255)
                 tCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATttgtt                                                  <  1:4415298/71‑1 (MQ=255)
                 tCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATttgtt                                                  <  1:3972393/71‑1 (MQ=255)
                    ccATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGcc                                               >  1:3906280/1‑71 (MQ=255)
                    ccATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGcc                                               <  1:624218/71‑1 (MQ=255)
                     cATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCg                                              >  1:689984/1‑71 (MQ=255)
                       tCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAg                                                              <  1:5521815/53‑1 (MQ=255)
                       tCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAg                                                              <  1:5167924/53‑1 (MQ=255)
                       tCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAg                                                              <  1:5394859/53‑1 (MQ=255)
                       tCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCAGTCTACTTTACCGTACAGATCTAAtt                                                      <  1:3704217/61‑1 (MQ=255)
                        cTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCt                                           >  1:2346539/1‑71 (MQ=255)
                        cTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCt                                           >  1:3663035/1‑71 (MQ=255)
                         tACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTa                                                                 >  1:4845125/1‑48 (MQ=255)
                         tACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTa                                                                 >  1:78174/1‑48 (MQ=255)
                         tACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCtt                                          >  1:431957/1‑71 (MQ=255)
                         tACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCtt                                          >  1:2625383/1‑71 (MQ=255)
                         tACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCtt                                          >  1:1934871/1‑71 (MQ=255)
                         tACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCtt                                          >  1:2906427/1‑71 (MQ=255)
                         tACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCtt                                          >  1:4393193/1‑71 (MQ=255)
                         tACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCtt                                          >  1:3347306/1‑71 (MQ=255)
                           cATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCTTCTttg                                        <  1:180053/71‑1 (MQ=255)
                           cATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCt                                           <  1:2752969/68‑1 (MQ=255)
                           cATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTttg                                        <  1:696502/71‑1 (MQ=255)
                           cATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTttg                                        <  1:3816773/71‑1 (MQ=255)
                            aTCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCg                                                                   <  1:4857748/43‑1 (MQ=255)
                            aTCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCg                                                                   <  1:914361/43‑1 (MQ=255)
                             tCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAAt                                                       <  1:780484/54‑1 (MQ=255)
                             tCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATttgt                                                   >  1:4338540/1‑58 (MQ=255)
                                 tgttgtCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTa                                   >  1:2978238/1‑70 (MQ=255)
                                 tgttgtCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTa                                   <  1:4519547/70‑1 (MQ=255)
                                  gttgtCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTaa                                  <  1:3901941/70‑1 (MQ=255)
                                   ttatCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATaa                                  >  1:4584766/4‑69 (MQ=2)
                                   ttatCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAAc                                >  1:5075339/4‑71 (MQ=9)
                                    tgtCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTttgtt                                      >  1:5530487/1‑64 (MQ=255)
                                     atCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATttgt                                                   >  1:1072543/2‑50 (MQ=2)
                                      tCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACtt                              >  1:323894/1‑70 (MQ=255)
                                      tCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACtt                              <  1:2968300/70‑1 (MQ=255)
                                        agagAAATAGTGCAGGCCGTCTACTTTACCGTACa                                                               >  1:3978740/1‑35 (MQ=255)
                                        agagAAATAGTGCAGGCCGTCTACTTTACCGTACa                                                               >  1:3062996/1‑35 (MQ=255)
                                        agagAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAAt                                                       <  1:4981965/43‑1 (MQ=255)
                                        agagAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATttgt                                                   >  1:1637962/1‑47 (MQ=255)
                                        agagAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATttgt                                                   >  1:63663/1‑47 (MQ=255)
                                        agagAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTcag                           >  1:2106970/1‑71 (MQ=255)
                                        agagAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTcag                           >  1:4428153/1‑71 (MQ=255)
                                          aaaaaTAATGCAGGCCATCAACTTTGCCATACAGATCTAAttt                                                     <  1:3358368/41‑1 (MQ=2)
                                           gaAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTcagcag                        >  1:1661661/1‑71 (MQ=255)
                                           aaaaTAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGt                                             >  1:5335882/2‑50 (MQ=14)
                                           aaaaTAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTcggctg                        >  1:1430497/2‑69 (MQ=9)
                                            aaaTAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTc                             >  1:3737092/1‑65 (MQ=255)
                                                aGTGCAGGCCGTCTACTTTTCCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGttt                   >  1:1521748/1‑71 (MQ=255)
                                                aGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTca                            >  1:4284700/1‑62 (MQ=255)
                                                 gTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTaaa                                 >  1:2672177/1‑56 (MQ=255)
                                                 gTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTtg                  <  1:3563717/71‑1 (MQ=255)
                                                 gTGCAGGCCGTCTACTTTAACGTACAGATCTAATTTGTTGCCg                                              >  1:1951569/1‑43 (MQ=255)
                                                    cAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTca                            <  1:4235278/58‑1 (MQ=255)
                                                     aGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTtgctgc              <  1:482148/71‑1 (MQ=255)
                                                      ggCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTgcg             <  1:767531/71‑1 (MQ=255)
                                                      ggCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGtgcagca             <  1:2463377/71‑2 (MQ=12)
                                                          gTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCtt                                          <  1:1041773/38‑1 (MQ=255)
                                                           tCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTg         <  1:3581500/70‑1 (MQ=255)
                                                             tACTTTTCCGTACAAATCCAGTTTATTACCATCTTTATTATAGACTTCAGCCGCATTTGCTGCGCCAGc        <  1:3016962/69‑1 (MQ=2)
                                                             tACTTTACCGTACAGATCTAATTTGTTGCCGTCTttgtt                                      >  1:605422/1‑39 (MQ=255)
                                                             tACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTca                            >  1:5107266/1‑49 (MQ=255)
                                                             tACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTc                             >  1:2609717/1‑48 (MQ=255)
                                                             tACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGttt                   <  1:3279653/58‑1 (MQ=255)
                                                             tACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTa      >  1:1528568/1‑71 (MQ=255)
                                                             aaCTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAAc                                >  1:5362326/2‑45 (MQ=17)
                                                             aaCTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAAc                                >  1:1064332/2‑45 (MQ=17)
                                                             aaCTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAAc                                >  1:4446776/2‑45 (MQ=17)
                                                             aaCTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAAc                                >  1:1861550/2‑45 (MQ=17)
                                                             aaCTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTcg                            >  1:2065895/2‑48 (MQ=17)
                                                              aCTTTTCCGTACAAATCCAGTTTATTACCATCTTTATTATAGACTTCAGCCGCATTTGCTGCGcc           >  1:470653/1‑65 (MQ=2)
                                                              aCTTTTCCGTACAAATCCAGTTTATTACCATCTTTATTATAGACTTCAGCCGCATTTGCTGCGcc           >  1:2294653/1‑65 (MQ=2)
                                                              aCTTTTCCGTACAAATCCAGTTTATTACCATCTTTATTATAGACTTCAGCCGCATTTGCTGCGcc           >  1:165270/1‑65 (MQ=2)
                                                              aCTTTTCCGTACAAATCCAGTTTATTACCATCTTTATTATAGACTTCAGCCGCATTTGCTGCGcc           >  1:1285980/1‑65 (MQ=2)
                                                               cTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTa      <  1:4598928/69‑1 (MQ=255)
                                                                tttACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTgcgc            <  1:4836368/62‑1 (MQ=255)
                                                                tttACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTgcgc            <  1:934179/62‑1 (MQ=255)
                                                                tttACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTAc     >  1:4481721/1‑69 (MQ=255)
                                                                  tACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcag  >  1:520484/1‑70 (MQ=255)
                                                                  tACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcag  >  1:496206/1‑70 (MQ=255)
                                                                  tACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcag  >  1:4389517/1‑70 (MQ=255)
                                                                  tACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcag  >  1:1937160/1‑70 (MQ=255)
                                                                  tACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcag  >  1:1621446/1‑70 (MQ=255)
                                                                      |                                                                 
GACGCATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAG  >  Exported/1440127‑1440262

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: