Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 455,055 A→C 16.1% A10A (GCT→GCG ybgH ← predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported455,0550AC16.1% 158.1 / 31.9 81A10A (GCT→GCGybgHpredicted transporter
Reads supporting (aligned to +/- strand):  ref base A (34/34);  new base C (5/8);  total (39/42)
Fisher's exact test for biased strand distribution p-value = 5.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

ATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATAGTT  >  Exported/454987‑455122
                                                                    |                                                                   
aTCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAg                                                                    <  1:1708427/70‑1 (MQ=255)
  cagcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgc                                                                 >  1:4318729/1‑71 (MQ=255)
  cagcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgc                                                                 >  1:2903474/1‑71 (MQ=255)
   agcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgcg                                                                >  1:626790/1‑71 (MQ=255)
   agcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgcg                                                                >  1:1768082/1‑71 (MQ=255)
     taaggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATcgcgcgcg                                                              >  1:4027102/4‑71 (MQ=255)
     cagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGcgcgcg                                                              >  1:925877/1‑71 (MQ=255)
      agGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCgg                                                             >  1:2968215/1‑71 (MQ=255)
      agGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCgg                                                             >  1:1929373/1‑71 (MQ=255)
      aaggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCgg                                                             >  1:1049272/3‑71 (MQ=255)
      aaggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCgg                                                             >  1:7948/3‑71 (MQ=255)
         gCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtg                                                          <  1:3150876/71‑1 (MQ=255)
          cACGCATGCCGTAAAAGCTAAAATATTTCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                                                         <  1:1217570/71‑1 (MQ=255)
          cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAg                                                                    >  1:957106/1‑60 (MQ=255)
          cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                                                         <  1:4503315/71‑1 (MQ=255)
          cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                                                         <  1:2786885/71‑1 (MQ=255)
          cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                                                         <  1:4509645/71‑1 (MQ=255)
          cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtg                                                          <  1:1933962/70‑1 (MQ=255)
          cACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCtgt                                                         <  1:3778408/71‑1 (MQ=255)
            cGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGa                                                       <  1:70855/71‑1 (MQ=255)
            cGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGa                                                       <  1:294826/71‑1 (MQ=255)
             gCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGag                                                      >  1:2119683/1‑70 (MQ=255)
              cATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGTTGTGATg                                                     >  1:3128724/1‑71 (MQ=255)
              cATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGagg                                                     <  1:121246/71‑3 (MQ=255)
                tGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGagg                                                     <  1:1998365/69‑3 (MQ=255)
                  ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgtg                                                 >  1:442478/1‑71 (MQ=255)
                     tAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCtgt                                                         <  1:5459242/60‑1 (MQ=255)
                       aaaaCTGAAGTACTCCCAGATTTGTATCGCCACGATATAGTATATCGCGCGCGGCTGTGag                                                      >  1:3834926/1‑60 (MQ=255)
                       aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGttttc                                            <  1:4056179/70‑3 (MQ=255)
                       aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGttttc                                            <  1:1798406/70‑3 (MQ=255)
                       aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGttttc                                            <  1:3652356/70‑3 (MQ=255)
                       aaaGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAg                                                                    >  1:1839905/1‑47 (MQ=255)
                       aaaGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTAt                                            >  1:3473302/1‑71 (MQ=255)
                          gCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGc                                                            <  1:2532398/52‑1 (MQ=255)
                            tAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                                                         <  1:515501/53‑1 (MQ=255)
                             aaaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTa                                             >  1:4483838/1‑64 (MQ=255)
                              aaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTa                                             >  1:5443618/1‑63 (MQ=255)
                              aaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCataat                                     <  1:4807160/71‑1 (MQ=255)
                                 tatTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGtt                                               <  1:1532654/58‑1 (MQ=255)
                                   ttCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTc                                <  1:4033631/71‑1 (MQ=255)
                                      ccAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTcc                               >  1:4394748/1‑69 (MQ=255)
                                      ccAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCcac                             <  1:5599699/71‑1 (MQ=255)
                                      ccAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCca                              >  1:3020669/1‑70 (MQ=255)
                                      ccAGATTTGCAGAGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCca                              >  1:1175136/1‑70 (MQ=255)
                                       cAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCcaca                            <  1:1510223/71‑1 (MQ=255)
                                        aGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGa                                                       >  1:2925206/1‑43 (MQ=255)
                                        aGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGa                                                       >  1:3796198/1‑43 (MQ=255)
                                        aGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACaa                           >  1:4977867/1‑71 (MQ=255)
                                         gATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAAt                          >  1:5553338/1‑71 (MQ=255)
                                          aTTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAAtt                         <  1:5153035/71‑1 (MQ=255)
                                          aTTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAAtt                         >  1:1924153/1‑71 (MQ=255)
                                          aTTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAAtt                         <  1:3458865/71‑1 (MQ=255)
                                            ttGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGa                       >  1:5031569/1‑71 (MQ=255)
                                            ttGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGa                       >  1:2460373/1‑71 (MQ=255)
                                            ttGCAGAGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGa                       >  1:737611/1‑71 (MQ=255)
                                              gCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGa                       <  1:5020296/69‑1 (MQ=255)
                                               cAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACgg                    >  1:3548632/1‑71 (MQ=255)
                                                aGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGc                   <  1:1001850/71‑1 (MQ=255)
                                                 gcgcAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGc                   <  1:4602369/70‑1 (MQ=255)
                                                  cgcCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGttttc                                            >  1:2435057/1‑41 (MQ=12)
                                                  cgcAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCaa                 >  1:4018329/1‑71 (MQ=255)
                                                      aCGACATAGTAAATAGCGCGCGGCTGTGATGCGTGtt                                               <  1:940477/37‑1 (MQ=255)
                                                      aCGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTcc                               <  1:2095013/53‑1 (MQ=255)
                                                      aCGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTcc                               <  1:1744881/53‑1 (MQ=255)
                                                       cGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAAtt                         <  1:3388961/58‑1 (MQ=255)
                                                       cGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAAtt                         <  1:5311807/58‑1 (MQ=255)
                                                        gACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCcac                             >  1:1080846/1‑53 (MQ=255)
                                                        gACATAATAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCAc           <  1:1200868/71‑1 (MQ=255)
                                                         aCATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTc             <  1:4935546/68‑1 (MQ=255)
                                                          cATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGt         <  1:4355098/71‑1 (MQ=255)
                                                           aTAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTa        <  1:5211152/71‑1 (MQ=255)
                                                           aTAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTa        <  1:1986489/71‑1 (MQ=255)
                                                              gTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCataat                                     <  1:923398/39‑1 (MQ=255)
                                                              gTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTaa       >  1:5243582/1‑69 (MQ=255)
                                                              gTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTaa       >  1:3092661/1‑69 (MQ=255)
                                                               tAAATAGCGCGCGGCTGTGATGCGTGTTTATTCAtaataa                                   >  1:512561/1‑40 (MQ=255)
                                                               tAAATAGCGCGCGGCTGTGATGCGTGTTTATTCAtaataa                                   >  1:389757/1‑40 (MQ=255)
                                                               tAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTaa       <  1:3794278/68‑1 (MQ=255)
                                                               tAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTaa       <  1:4625531/68‑1 (MQ=255)
                                                               tAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTaa       <  1:1133868/68‑1 (MQ=255)
                                                                aaaTAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATAGt   >  1:2178223/1‑71 (MQ=255)
                                                                 aaTAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATAGtt  >  1:457941/1‑71 (MQ=255)
                                                                    aGCGCGCTGCTGTGATGCGTGTTTATTCATAATAATTCCACAAtt                         <  1:3459574/45‑1 (MQ=255)
                                                                    aGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATa     >  1:5512767/1‑65 (MQ=255)
                                                                    aGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATa     >  1:3374281/1‑65 (MQ=255)
                                                                    aGCGCGCGGCTGTGATGCGTGTTAATTCATAATAATTCCACAAtt                         >  1:3020230/1‑45 (MQ=255)
                                                                    |                                                                   
ATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATAGTT  >  Exported/454987‑455122

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: