Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,864,289 G→C 7.0% S185S (TCC→TCG mtlR ← DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,864,2890GC 7.0% 210.7 / 5.1 71S185S (TCC→TCGmtlRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (38/28);  new base C (2/3);  total (40/31)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

GGCCGAAGCAGGTAGAGGGGCAGTAAGCAGACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCGCACAGCCTGTTGATAGCGCTGAATTTG  >  AP009048/3864219‑3864354
                                                                      |                                                                 
ggCCGAAGCAGGTAGAGGGGCAGTAAGCAGACTTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTg                                                                   <  1:1248760/71‑1 (MQ=255)
 gCCGAAGCAGGTAGAGGGGCAGTAAGCAGACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTg                                                                   >  1:619431/1‑70 (MQ=255)
   cGAAGCAGGTAGAGGGGCAGTAAGCAGACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGaa                                                                <  1:3871603/71‑1 (MQ=255)
    gAAGCAGGTAGAGGGGCAGTAAGCAGACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGaaa                                                               <  1:1685423/71‑1 (MQ=255)
     aaGCAGGTAGAGGGGCAGTAAGCAGACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAAt                                                              >  1:1057622/1‑71 (MQ=255)
     aaGCAGGTAGAGGGGCAGTAAGCAGACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAAt                                                              >  1:339987/1‑71 (MQ=255)
      aGCAGGTAGAGGGGCAGTAAGCAGACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAAt                                                              >  1:932681/1‑70 (MQ=255)
        cAGGTAGAGGGGCAGTAAGCAGACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAg                                                           <  1:798164/71‑1 (MQ=255)
        cAGGTAGAGGGGCAGTAAGCAGACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAg                                                           <  1:4981206/71‑1 (MQ=255)
             agagGGGCAGTAAGCAGACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCtcag                                                      >  1:805301/1‑71 (MQ=255)
             agagGGGCAGTAAGCAGACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCtcag                                                      >  1:4015245/1‑71 (MQ=255)
               agGGGCAGTAAGCAGACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCtcagt                                                     >  1:3318741/1‑70 (MQ=255)
                   gCAGTAAGCAGACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGa                                                                 <  1:129368/54‑1 (MQ=255)
                    cAGTAAGCAGACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGa                                               <  1:2755227/71‑1 (MQ=255)
                    cAGTAAGCAGACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGa                                               <  1:4453672/71‑1 (MQ=255)
                    cAGTAAGCAGACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGa                                               <  1:1957853/71‑1 (MQ=255)
                    cAGTAAGCAGACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGa                                               <  1:4060674/71‑1 (MQ=255)
                      gTAAGCAGACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGa                                               >  1:4373337/1‑69 (MQ=255)
                       tAAGCAGACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGAAAg                                            <  1:786948/71‑1 (MQ=255)
                        aaGCAGACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGc                                                          >  1:3730013/1‑56 (MQ=255)
                         aGCAGACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGa                                               >  1:1766024/1‑66 (MQ=255)
                         aGCAGACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGa                                          >  1:3940401/1‑71 (MQ=255)
                         aGCAGACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGa                                          >  1:4736572/1‑71 (MQ=255)
                         aGCAGACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTCCAAAATCAGc                                                          >  1:3560880/1‑55 (MQ=37)
                         aGCAGACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTCCAAAATCAGc                                                          >  1:1024181/1‑55 (MQ=37)
                          gCAGACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTCCAAAATCa                                                            <  1:1574789/52‑1 (MQ=37)
                             gACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCtcagtc                                                    <  1:1089161/57‑1 (MQ=255)
                             gACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCAt                                      >  1:993254/1‑71 (MQ=255)
                             gACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCAt                                      >  1:2438483/1‑71 (MQ=255)
                             gACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCAt                                      >  1:1314826/1‑71 (MQ=255)
                                  tGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCAttgtt                                  >  1:1060230/1‑70 (MQ=255)
                                  tGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCAttgtt                                  >  1:84736/1‑70 (MQ=255)
                                      tACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGa                                <  1:425457/68‑1 (MQ=255)
                                      tACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCg                             >  1:2313886/1‑71 (MQ=255)
                                       aCTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGaa                                              >  1:2173953/1‑53 (MQ=255)
                                       aCTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGa                                               >  1:3410272/1‑52 (MQ=255)
                                       aCTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCAttgttg                                 >  1:4066434/1‑66 (MQ=255)
                                       aCTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCg                             >  1:2192971/1‑70 (MQ=255)
                                       aCTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCg                             >  1:4109333/1‑70 (MQ=255)
                                          tttGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGa                                               <  1:930048/49‑1 (MQ=255)
                                          tttGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCGca                           >  1:1326581/1‑69 (MQ=255)
                                           ttGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCGcac                          <  1:2012717/69‑1 (MQ=255)
                                           ttGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCGCACAg                        >  1:1376303/1‑71 (MQ=255)
                                           ttGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCGCACAg                        >  1:1452843/1‑71 (MQ=255)
                                              aaaGGCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGAc                                         >  1:2088397/1‑51 (MQ=255)
                                               aaGGCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCGCACAGCCTg                    >  1:3801019/1‑71 (MQ=255)
                                                  gCTTTTTTTAAGCTGATTTTGGAAATCAGCtcagtcag                                                  <  1:1171897/38‑1 (MQ=255)
                                                  gCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCGCACAGCCTGTTg                 <  1:4083739/71‑1 (MQ=255)
                                                     ttttttAAGCTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCGcaca                         >  1:2195712/1‑60 (MQ=255)
                                                     ttttttAAGCTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCGCACAGCCTGTTGa                >  1:2401539/1‑69 (MQ=255)
                                                     ttttttAAGCTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCGCACAGCCTGTTGa                >  1:2988837/1‑69 (MQ=255)
                                                     ttttttAAGCTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCGCACAGCCTGTTGa                >  1:2553281/1‑69 (MQ=255)
                                                      tttttAAGCTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCg                             <  1:168219/55‑1 (MQ=255)
                                                          tAAGCTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGAt                               >  1:3673666/1‑49 (MQ=255)
                                                          tAAGCTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGAt                               >  1:1378079/1‑49 (MQ=255)
                                                           aaGCTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGa                                >  1:83550/1‑47 (MQ=255)
                                                              cTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCa                                       <  1:2470422/37‑1 (MQ=255)
                                                              cTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCGCACAGCCTGTTGATAgcg           <  1:1368498/65‑1 (MQ=255)
                                                              cTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCGCACAGCCTGTTGATAGCGCTGAAt     <  1:4260996/71‑1 (MQ=255)
                                                               tGATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCGCACAg                        <  1:232379/51‑1 (MQ=255)
                                                                gATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCGCGCAGCCTGTTGATAGCGCTGAAttt   >  1:4324860/1‑71 (MQ=255)
                                                                gATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCGCACAGCCTGTTg                 >  1:4505317/1‑57 (MQ=255)
                                                                gATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCGCACAGCCTGTTGATAGCGCTGAAt     <  1:1524511/69‑1 (MQ=255)
                                                                 aTTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCGCACAg                        >  1:1630105/1‑49 (MQ=255)
                                                                 aTTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCGCACAGCCTGTTGAt               <  1:1279280/58‑1 (MQ=255)
                                                                 aTTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCGCACAGCCTGTTGATAGCGCTGAAttt   <  1:2542662/70‑1 (MQ=255)
                                                                 aTTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCGCACAGCCTGTTGATAGCGCTGAAttt   <  1:3779454/70‑1 (MQ=255)
                                                                 aTTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCGCACAGCCTGTTGATAGCGCTGAAtt    <  1:2787139/69‑1 (MQ=255)
                                                                 aTTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCGCACAGCCTGTTGATAGCGCTGAAt     <  1:4701073/68‑1 (MQ=255)
                                                                 aTTTTCGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCGCACAGCCTGTTGATAGCGCTGAAttt   <  1:2341083/70‑1 (MQ=255)
                                                                 aTTTTCGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCGCACAGCCTGTTGATAGCGCTGAAttt   <  1:2341081/70‑1 (MQ=255)
                                                                  ttttGGAAATCAGCTCAGTCAGTGAAAGGACCAttgt                                   <  1:624689/37‑1 (MQ=255)
                                                                  ttttGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCGCACAGCCTGTTGATAGCGCTGAATttg  >  1:4460851/1‑70 (MQ=255)
                                                                      |                                                                 
GGCCGAAGCAGGTAGAGGGGCAGTAAGCAGACGTTGCTTACTTTTGAAAGGCTTTTTTTAAGCTGATTTTGGAAATCAGCTCAGTCAGTGAAAGGACCATTGTTGATCGCACAGCCTGTTGATAGCGCTGAATTTG  >  AP009048/3864219‑3864354

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: