Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,440,236 A→G 11.7% 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
*Exported1,440,2360AG11.7% 210.1 / 16.3 85A22A (GCT→GCCompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base A (46/29);  new base G (4/6);  total (50/35)
Fisher's exact test for biased strand distribution p-value = 3.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGTTAT  >  Exported/1440167‑1440306
                                                                     |                                                                      
agagAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTcag                                                                       >  1:2106970/1‑71 (MQ=255)
agagAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTcag                                                                       >  1:4428153/1‑71 (MQ=255)
   gaAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTcagcag                                                                    >  1:1661661/1‑71 (MQ=255)
   aaaaTAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTcggctg                                                                    >  1:1430497/2‑69 (MQ=9)
        aGTGCAGGCCGTCTACTTTTCCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGttt                                                               >  1:1521748/1‑71 (MQ=255)
        aGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTca                                                                        >  1:4284700/1‑62 (MQ=255)
         gTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTtg                                                              <  1:3563717/71‑1 (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)
                   tCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTg                                                     <  1:3581500/70‑1 (MQ=255)
                     tACTTTTCCGTACAAATCCAGTTTATTACCATCTTTATTATAGACTTCAGCCGCATTTGCTGCGCCAGc                                                    <  1:3016962/69‑1 (MQ=2)
                     tACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTca                                                                        >  1:5107266/1‑49 (MQ=255)
                     tACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGttt                                                               <  1:3279653/58‑1 (MQ=255)
                     tACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTa                                                  >  1:1528568/1‑71 (MQ=255)
                      aCTTTTCCGTACAAATCCAGTTTATTACCATCTTTATTATAGACTTCAGCCGCATTTGCTGCGcc                                                       >  1:470653/1‑65 (MQ=2)
                      aCTTTTCCGTACAAATCCAGTTTATTACCATCTTTATTATAGACTTCAGCCGCATTTGCTGCGcc                                                       >  1:165270/1‑65 (MQ=2)
                      aCTTTTCCGTACAAATCCAGTTTATTACCATCTTTATTATAGACTTCAGCCGCATTTGCTGCGcc                                                       >  1:1285980/1‑65 (MQ=2)
                      aCTTTTCCGTACAAATCCAGTTTATTACCATCTTTATTATAGACTTCAGCCGCATTTGCTGCGcc                                                       >  1:2294653/1‑65 (MQ=2)
                       cTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTa                                                  <  1:4598928/69‑1 (MQ=255)
                        tttACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTgcgc                                                        <  1:934179/62‑1 (MQ=255)
                        tttACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTgcgc                                                        <  1:4836368/62‑1 (MQ=255)
                        tttACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTAc                                                 >  1:4481721/1‑69 (MQ=255)
                          tACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcag                                              >  1:4389517/1‑70 (MQ=255)
                          tACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcag                                              >  1:520484/1‑70 (MQ=255)
                          tACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcag                                              >  1:1937160/1‑70 (MQ=255)
                          tACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcag                                              >  1:1621446/1‑70 (MQ=255)
                          tACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcag                                              >  1:496206/1‑70 (MQ=255)
                            ccatACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGTGCAGCACCTGCGGCgagcag                                           <  1:3618084/68‑1 (MQ=255)
                               tACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTcagca                                                                     <  1:5490912/42‑1 (MQ=255)
                               tACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTcagca                                                                     <  1:1616235/42‑1 (MQ=255)
                               tACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGt                                                                 >  1:4747511/1‑46 (MQ=255)
                               tACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTtgct                                                            <  1:4509361/51‑1 (MQ=255)
                               tACAGATCTAATTTGTTCCGTCTTTGTTGTAAACTTcagca                                                                     <  1:4550452/41‑1 (MQ=39)
                                aCAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcag                                              >  1:2338551/1‑64 (MQ=255)
                                aCAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGTGCAGCACCTGcggc                                                 <  1:1619161/61‑4 (MQ=255)
                                aCAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGTGCAGCACCTGCGGCgagcag                                           <  1:235588/67‑1 (MQ=255)
                                 cAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTg                                      <  1:1884941/71‑1 (MQ=255)
                                  aGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGttt                                                               <  1:3575862/45‑1 (MQ=255)
                                  aGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTgg                                     >  1:4972598/1‑71 (MQ=255)
                                  aGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTgg                                     >  1:1228954/1‑71 (MQ=255)
                                   gATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGtgcagc                                                          <  1:5338854/49‑1 (MQ=16)
                                    aTCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCTTTtgc                                                             >  1:4884715/1‑45 (MQ=255)
                                    aTCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTtgc                                                             >  1:1788181/1‑45 (MQ=255)
                                     tCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGtt                                                                >  1:4738667/1‑41 (MQ=255)
                                      cTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGTGCAGCACCTGCGGCGAgcagggcaggaat                                  >  1:4935129/1‑67 (MQ=255)
                                      cTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGTGCAGCACCTGCGGCGAgcagggcaggaat                                  >  1:2648995/1‑67 (MQ=255)
                                       tAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGt                                                                 >  1:439235/1‑38 (MQ=255)
                                       tAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGt                                                                 >  1:1495279/1‑38 (MQ=255)
                                       tAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGTGCAGCACCt                                                      <  1:1763009/49‑1 (MQ=12)
                                        aaTTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTtg                                                              >  1:5116743/1‑40 (MQ=255)
                                        aaTTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCag                               >  1:2153372/1‑71 (MQ=255)
                                         aTTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGAc                                  <  1:788333/67‑1 (MQ=255)
                                         aTTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGtgcagcac                                                        >  1:1627459/1‑43 (MQ=9)
                                         aTTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGtgcagca                                                         >  1:2573411/1‑43 (MQ=9)
                                          tttGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGAcc                                 <  1:5335710/67‑1 (MQ=255)
                                          tttGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCagga                             <  1:3937048/71‑1 (MQ=255)
                                          tttGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCagga                             <  1:204270/71‑1 (MQ=255)
                                          tttGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCag                               >  1:2103598/1‑69 (MQ=255)
                                            tgttgCCGTCTTTGTTGTAAACTTCAGCAGCGTTtg                                                              >  1:2438350/1‑36 (MQ=255)
                                            tgttgCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCaggagg                           >  1:1721441/1‑71 (MQ=255)
                                             gttgCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCaggagg                           <  1:988937/70‑1 (MQ=255)
                                             gttgCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCaggagg                           <  1:3760441/70‑1 (MQ=255)
                                                 ccGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAg                      <  1:75376/71‑1 (MQ=255)
                                                   gTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTa                    >  1:3093427/1‑71 (MQ=255)
                                                    tCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGCGCTGGGACCAGGAGGGACAGTAc                   >  1:2184546/1‑71 (MQ=255)
                                                    tCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGa                                   <  1:2971903/55‑1 (MQ=255)
                                                    tCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGAcc                                 <  1:3814144/57‑1 (MQ=255)
                                                    tCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTAc                   >  1:1652277/1‑71 (MQ=255)
                                                     cTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTa                                                  <  1:4306071/39‑1 (MQ=255)
                                                      tttGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCag                               <  1:4770638/57‑1 (MQ=255)
                                                      tttGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACtt                 >  1:5492746/1‑71 (MQ=255)
                                                      tttGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACtt                 >  1:4708676/1‑71 (MQ=255)
                                                      tttGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACtt                 >  1:3535897/1‑71 (MQ=255)
                                                      tttGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACtt                 >  1:2461263/1‑71 (MQ=255)
                                                        tgttgtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGATCTGGGACCAGGAGGGACAGTACttt                <  1:1410667/70‑1 (MQ=255)
                                                        tgttgtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACttt                >  1:4949594/1‑70 (MQ=255)
                                                        tgttgtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACttt                >  1:3094092/1‑70 (MQ=255)
                                                        tgttgtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACttt                <  1:169379/70‑1 (MQ=255)
                                                        tgttgtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACttt                >  1:3792894/1‑70 (MQ=255)
                                                        tgttgtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACtt                 <  1:3888863/69‑1 (MQ=255)
                                                           tgtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACca                                               >  1:1156736/1‑36 (MQ=255)
                                                           tgtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACca                                               >  1:3527272/1‑36 (MQ=255)
                                                           tgtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACca                                               >  1:4931922/1‑36 (MQ=255)
                                                            gtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTAc                   >  1:4123593/1‑63 (MQ=255)
                                                            gtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTAc                   >  1:3194175/1‑63 (MQ=255)
                                                            gtAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACtt           >  1:1963810/1‑71 (MQ=255)
                                                             tAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACttt          >  1:3661286/1‑71 (MQ=255)
                                                             tAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACttt          <  1:314447/71‑1 (MQ=255)
                                                             tAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACtt           <  1:649078/70‑1 (MQ=255)
                                                              aaaCTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGa                                   >  1:1004588/1‑45 (MQ=255)
                                                               aaCTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGc                                        >  1:5261197/1‑39 (MQ=255)
                                                               aaCTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTc         <  1:2271148/70‑1 (MQ=255)
                                                                   tCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGt     >  1:3377800/1‑70 (MQ=255)
                                                                    cagcagCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGtta   >  1:772416/1‑71 (MQ=255)
                                                                    cagcagCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGt     <  1:3587214/69‑1 (MQ=255)
                                                                    cagcagCGTTTGATGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGtta   >  1:3994858/1‑71 (MQ=255)
                                                                     agcagcGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGttat  <  1:375657/71‑1 (MQ=255)
                                                                     agcagcGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGt     <  1:2394222/68‑1 (MQ=255)
                                                                     |                                                                      
AGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAGCAGAGCTGGGACCAGGAGGGACAGTACTTTAACTTTCATGTTAT  >  Exported/1440167‑1440306

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: