Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,258,449 A→G 16.3% intergenic (+120/+45) 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,4490AG16.3% 135.6 / 16.0 55intergenic (+120/+45)yecP/torZpredicted methyltransferase/trimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base A (23/23);  new base G (3/6);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGAT  >  Exported/1258382‑1258516
                                                                   |                                                                   
aCCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGt                                                                  <  1:4203399/71‑1 (MQ=255)
aCCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGt                                                                  <  1:2336628/71‑1 (MQ=255)
 ccTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGtt                                                                 >  1:1082953/1‑71 (MQ=255)
    tGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGttt                                                                >  1:1692566/1‑69 (MQ=255)
      cTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCAt                                                            <  1:3921158/71‑1 (MQ=255)
      cTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCAt                                                            <  1:3633435/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)
         aTTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATaa                                                          >  1:4601760/1‑70 (MQ=255)
            cAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACg                                                                   >  1:1236330/1‑58 (MQ=255)
            cAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCt                                                       <  1:943738/70‑1 (MQ=255)
             aaTATATTTAATGGGTTATGTAGACCTGATAAGCATAGCGCATCAGGCTGCTTTGCGt                                                                  >  1:3009687/1‑58 (MQ=9)
             aaTATATTTAATGGGTTATGTAGACCTGATAAGCATAGCGCATCAGGCTGCTTTGCGt                                                                  >  1:1525217/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:4712131/39‑1 (MQ=11)
                               tGTAGACCTGATAAGCATAGCGCATCAGGCTGCTTTGCg                                                                   <  1:1362744/39‑1 (MQ=11)
                               tGTAGACCTGATAAGCATAGCGCATCAGGCTGCTTTGCg                                                                   <  1:2256842/39‑1 (MQ=11)
                               tGTAGACCTGATAAGCATAGCGCATCAGGCTGCTTTGCg                                                                   <  1:2115693/39‑1 (MQ=11)
                                gTAGACCTGATAAGCATAGCGCATCAGGCTGCTTTGCGttt                                                                <  1:420357/41‑1 (MQ=14)
                                gTAGACCTGATAAGCATAGCGCATCAGGCTGCTTTGCGttt                                                                <  1:5264968/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:5626581/1‑35 (MQ=255)
                                         aTAAGCATAGCGCATCAGGCAGCTTTACGTTTGCa                                                             >  1:2680155/1‑35 (MQ=255)
                                         aTAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTc                         >  1:5153246/1‑71 (MQ=255)
                                            aGCATAGCGCATCAGGCAGCTTTACGTTTGCATaa                                                          >  1:1141085/1‑35 (MQ=255)
                                            aGCATAGCGCATCAGGCAGCTTTACGTTTGCATaa                                                          >  1:440664/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:2140568/1‑58 (MQ=255)
                                                  gcgcATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCa                             >  1:1511959/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)
                                                                tttACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGAt  >  1:502412/1‑71 (MQ=255)
                                                                   |                                                                   
ACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGAT  >  Exported/1258382‑1258516

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: