Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,258,443 A→T 18.0% intergenic (+114/+51) yecP → / ← torZ predicted methyltransferase/trimethylamine N‑oxide reductase system III, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,258,4430AT18.0% 152.1 / 24.4 61intergenic (+114/+51)yecP/torZpredicted methyltransferase/trimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base A (25/25);  new base T (3/8);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

ACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAGCTGG  >  Exported/1258377‑1258511
                                                                  |                                                                    
aCAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCtt                                                                  >  1:4506173/1‑71 (MQ=255)
  aaaaCCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCttt                                                                 <  1:1847400/70‑1 (MQ=255)
     aCCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGt                                                             <  1:2336628/71‑1 (MQ=255)
     aCCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGt                                                             <  1:4203399/71‑1 (MQ=255)
      ccTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGtt                                                            >  1:1082953/1‑71 (MQ=255)
         tGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGttt                                                           >  1:1692566/1‑69 (MQ=255)
           cTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCAt                                                       <  1:3633435/71‑1 (MQ=255)
           cTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCAt                                                       <  1:3921158/71‑1 (MQ=255)
            tAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCAt                                                       >  1:2458677/1‑70 (MQ=255)
            tAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCAt                                                       >  1:3643659/1‑70 (MQ=255)
            aagttCAATATATTTAATGGGTTATGTAGACCTGATAAGCATAGCGCATCAGGCTGCTTTGCGttt                                                           >  1:2317224/4‑66 (MQ=9)
             aaTTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCtt                                                                  >  1:831128/1‑58 (MQ=255)
             aaTTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCtt                                                                  >  1:3091987/1‑58 (MQ=255)
              aTTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATaa                                                     >  1:4601760/1‑70 (MQ=255)
               ttCAATATATTTAATGGGTTATGTAGACCTGATAAGCATAGCGCATCAGGCTGCtt                                                                  <  1:1378907/56‑1 (MQ=9)
               ttCAATATATTTAATGGGTTATGTAGACCTGATAAGCATAGCGCATCAGGCTGCtt                                                                  <  1:5190587/56‑1 (MQ=9)
                tCAATATATTTAATGGGTTATGTAGACCTGATAAGCATAGCGCATCAGgctgct                                                                   <  1:2949309/54‑1 (MQ=9)
                 cAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACg                                                              >  1:1236330/1‑58 (MQ=255)
                 cAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCt                                                  <  1:943738/70‑1 (MQ=255)
                  aaTATATTTAATGGGTTATGTAGACCTGATAAGCATAGCGCATCAGGCTGCTTTGCGt                                                             >  1:1525217/1‑58 (MQ=9)
                  aaTATATTTAATGGGTTATGTAGACCTGATAAGCATAGCGCATCAGGCTGCTTTGCGt                                                             >  1:3009687/1‑58 (MQ=9)
                   atatatTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAg                                               <  1:357589/71‑1 (MQ=255)
                   atatatTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAg                                               <  1:2701218/71‑1 (MQ=255)
                    tatatTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCa                                                >  1:948591/1‑69 (MQ=255)
                       atTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAgcgc                                            <  1:3789462/70‑1 (MQ=255)
                        ttGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGccc                                          >  1:5067378/1‑71 (MQ=255)
                           cAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGtt                                       <  1:3583970/71‑1 (MQ=255)
                             gggACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTc                                     <  1:1467958/71‑1 (MQ=255)
                              ggACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTg                                                          <  1:2382131/49‑1 (MQ=255)
                                  tATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAgc                                              <  1:5310211/57‑1 (MQ=255)
                                   aTGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTa                                                                <  1:2472096/38‑1 (MQ=255)
                                    tGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGgcgc                              <  1:4297408/71‑1 (MQ=255)
                                    tGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGgc                                <  1:1599104/69‑1 (MQ=255)
                                    tGTAGACCTGATAAGCATAGCGCATCAGGCTGCTTTGCg                                                              <  1:2115693/39‑1 (MQ=11)
                                    tGTAGACCTGATAAGCATAGCGCATCAGGCTGCTTTGCg                                                              <  1:2256842/39‑1 (MQ=11)
                                    tGTAGACCTGATAAGCATAGCGCATCAGGCTGCTTTGCg                                                              <  1:4712131/39‑1 (MQ=11)
                                    tGTAGACCTGATAAGCATAGCGCATCAGGCTGCTTTGCg                                                              <  1:1362744/39‑1 (MQ=11)
                                     gTAGACCTGATAAGCATAGCGCATCAGGCTGCTTTGCGttt                                                           <  1:5264968/41‑1 (MQ=14)
                                     gTAGACCTGATAAGCATAGCGCATCAGGCTGCTTTGCGttt                                                           <  1:420357/41‑1 (MQ=14)
                                        ggCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTAtt                          <  1:1631156/71‑1 (MQ=255)
                                         gCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTc                                     <  1:3030715/59‑1 (MQ=255)
                                         gCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTc                                     <  1:3712108/59‑1 (MQ=255)
                                         gCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTc                         <  1:1320987/71‑1 (MQ=255)
                                           cTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCAc                       >  1:3640606/1‑71 (MQ=255)
                                            tGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTAtt                          >  1:183704/1‑67 (MQ=255)
                                            tGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTAtt                          >  1:3065350/1‑67 (MQ=255)
                                              aTAAGCATAGCGCATCAGGCAGCTTTACGTTTGCa                                                        >  1:2680155/1‑35 (MQ=255)
                                              aTAAGCATAGCGCATCAGGCAGCTTTACGTTTGCa                                                        >  1:5626581/1‑35 (MQ=255)
                                              aTAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTc                    >  1:5153246/1‑71 (MQ=255)
                                                 aGCATAGCGCATCAGGCAGCTTTACGTTTGCATaa                                                     >  1:440664/1‑35 (MQ=255)
                                                 aGCATAGCGCATCAGGCAGCTTTACGTTTGCATaa                                                     >  1:1141085/1‑35 (MQ=255)
                                                  gCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTAc                <  1:1335804/71‑1 (MQ=255)
                                                     tAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGt                     <  1:4147410/63‑1 (MQ=255)
                                                     tAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGcc             <  1:4514489/71‑1 (MQ=255)
                                                       gcgcATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCa                        >  1:1511959/1‑58 (MQ=255)
                                                       gcgcATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCa                        >  1:2140568/1‑58 (MQ=255)
                                                        cgcATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACg                      <  1:424829/59‑1 (MQ=255)
                                                          cATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAAc        <  1:3369029/71‑1 (MQ=255)
                                                              aGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAgct    >  1:4949797/1‑71 (MQ=255)
                                                              aGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAgct    >  1:1801471/1‑71 (MQ=255)
                                                              aGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAgct    >  1:1311196/1‑71 (MQ=255)
                                                                gCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAgctgg  >  1:4401814/1‑71 (MQ=255)
                                                                  |                                                                    
ACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAGCTGG  >  Exported/1258377‑1258511

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: