Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 236,569 C→T 25.4% A168V (GCC→GTC)  dnaQ → DNA polymerase III epsilon subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913236,5690CT25.4% 125.0 / 44.3 75A168V (GCC→GTC) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base C (24/32);  new base T (9/10);  total (33/42)
Fisher's exact test for biased strand distribution p-value = 7.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.89e-01

GGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGACACAACAGCAACAAGGTGAAGCAAC  >  NC_000913/236470‑236668
                                                                                                   |                                                                                                   
ggAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGc                                                                                                     <  2:971476/100‑1 (MQ=255)
ggAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGc                                                                                                     >  2:290075/1‑100 (MQ=255)
ggAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGc                                                                                                     <  1:929375/100‑1 (MQ=255)
  aaaaTGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGccc                                                                                                   <  2:2048468/100‑1 (MQ=255)
   aaaTGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCa                                                                                                  <  1:2432469/100‑1 (MQ=255)
   aaaTGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCa                                                                                                  <  2:2540237/100‑1 (MQ=255)
   aaaTGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCa                                                                                                  >  2:209462/1‑100 (MQ=255)
    aaTGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAg                                                                                                 >  1:671762/1‑100 (MQ=255)
     aTGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTTAACGAACGCTGCACGGGGCATTACTCGATGCCCAGa                                                                                                <  1:1834020/100‑1 (MQ=255)
      tGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCa                                                                                                  >  1:1115829/1‑97 (MQ=255)
           cccGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCtt                                                                                           >  2:999333/1‑99 (MQ=255)
            ccGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGc                                                                                         >  2:2669/1‑100 (MQ=255)
            ccGGTAAGCGCAACAGCCTCGATGCGTTATGTACTCGCGACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGc                                                                                         <  1:533288/100‑1 (MQ=255)
             cGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCg                                                                                        >  1:1673853/1‑100 (MQ=255)
             cGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCg                                                                                        <  1:1679715/100‑1 (MQ=255)
              ggTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCgg                                                                                       <  2:2602870/100‑1 (MQ=255)
              ggTAAGCGCAACAGCCTCGATGCGTTATGTGCTCCCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGc                                                                                         <  2:309001/98‑1 (MQ=255)
               gTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGa                                                                                      >  1:1454335/1‑100 (MQ=255)
                  aGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGt                                                                                   >  2:1490313/1‑100 (MQ=255)
                  aGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGt                                                                                   <  1:129339/100‑1 (MQ=255)
                   gcgcAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGtt                                                                                  <  1:1625457/100‑1 (MQ=255)
                    cgcAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGACCGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGttt                                                                                 <  2:414214/100‑1 (MQ=255)
                    cgcAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGttt                                                                                 <  1:690953/100‑1 (MQ=255)
                    cgcAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGtt                                                                                  <  2:1867707/99‑1 (MQ=255)
                    cgcAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGAATTACTCGATGCCCAGATCCTTGCGGAAGttt                                                                                 >  2:1435974/1‑100 (MQ=255)
                     gcAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGttt                                                                                 >  1:770332/1‑99 (MQ=255)
                      cAACAGCCTCGATGCGTTATGTGCTCGCTCCGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTTCTCGATGTCCAGATCCTTGCGGAAGTTTaa                                                                               <  1:999338/100‑2 (MQ=255)
                        aCAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCt                                                                             <  2:730979/100‑1 (MQ=255)
                           gCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGc                                                                          >  1:1261216/1‑100 (MQ=255)
                           gCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGc                                                                          >  2:1065034/1‑100 (MQ=255)
                           gCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGc                                                                          >  1:331268/1‑100 (MQ=255)
                            ccTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCg                                                                         <  2:85927/100‑1 (MQ=255)
                            ccTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCg                                                                         >  2:1056439/1‑100 (MQ=255)
                             cTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGa                                                                        <  1:389701/100‑1 (MQ=255)
                              tCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGAt                                                                       <  1:2456978/100‑1 (MQ=255)
                               cGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATg                                                                      >  1:2424163/1‑100 (MQ=255)
                                gATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGACATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATg                                                                      >  2:796557/1‑99 (MQ=255)
                                    cGTTATGTGCTCTCTACGAAATAGATAACAGTAAACAAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCgg                                                                 <  1:1082811/100‑1 (MQ=255)
                                    cGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCgg                                                                 >  1:1120348/1‑100 (MQ=255)
                                     gTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCggt                                                                >  2:146285/1‑100 (MQ=255)
                                     gTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCggt                                                                >  1:934210/1‑100 (MQ=255)
                                       tATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCggtgg                                                              <  1:2030787/100‑1 (MQ=255)
                                        aTGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCggtggt                                                             <  1:2539808/100‑1 (MQ=255)
                                         tgtgCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTc                                                            <  1:15720/100‑1 (MQ=255)
                                            gCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCaaa                                                         <  1:290078/100‑1 (MQ=255)
                                              tCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACg                                                       <  1:1651452/100‑1 (MQ=255)
                                                  tACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGa                                                   <  1:636639/100‑1 (MQ=255)
                                                   aCGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGa                                                   >  1:629118/1‑99 (MQ=255)
                                                    cGAAATAGATAACAGTAAACGAACGCTAAAAGGGGCATTACTCGATGTCCAGATCCTTGAGGAAGTTTATCTGACGATGACCAGTGGTCAAACGTCGATg                                                 >  2:2593686/1‑100 (MQ=255)
                                                           gATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTg                                          <  2:1747351/100‑1 (MQ=255)
                                                           gATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTg                                          <  1:1144778/100‑1 (MQ=255)
                                                            aTAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGc                                         <  2:2267132/100‑1 (MQ=255)
                                                             tAACAGTAAACGAACTCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTCCg                                        <  2:924656/100‑1 (MQ=255)
                                                             tAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCg                                        >  2:2045676/1‑100 (MQ=255)
                                                              aaCAGTAAACGGGTTCTGCACGGGGCATTACTCGATGTCCGGATCTTTGCGGAAGTTTATCTGGCGATGACCGGTGGTTAAACGTCGATGGCTTTTGCGa                                       <  2:1978292/100‑1 (MQ=255)
                                                               aCAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGAt                                      <  2:2478318/100‑1 (MQ=255)
                                                                cAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATg                                     >  2:2371928/1‑100 (MQ=255)
                                                                 aGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATgg                                    <  1:209464/100‑1 (MQ=255)
                                                                  gTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGa                                   >  1:2381904/1‑100 (MQ=255)
                                                                   tAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGaa                                  >  1:2451680/1‑100 (MQ=255)
                                                                          acgCTGCACGGGGCATTACTCGATGCCCAGATCCTGGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGgaga                             <  2:2446995/98‑1 (MQ=255)
                                                                           cgCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGac                          <  2:2245563/100‑1 (MQ=255)
                                                                             cTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGacac                        >  1:2102673/1‑100 (MQ=255)
                                                                              tGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGAcaca                       >  2:394463/1‑100 (MQ=255)
                                                                              tGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGAcaca                       <  1:304655/100‑1 (MQ=255)
                                                                              tGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGAcaca                       <  2:1804550/100‑1 (MQ=255)
                                                                               gCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGACacaa                      >  2:1475672/1‑100 (MQ=255)
                                                                                 aCGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGACacaaca                    <  2:970108/100‑1 (MQ=255)
                                                                                 aCGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGACacaaca                    >  2:2373286/1‑100 (MQ=255)
                                                                                  cGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGACACAACAg                   <  2:1770675/100‑1 (MQ=255)
                                                                                          tACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAGACGTCGATGGCTTTTGCGATGGAAGGAGTGACACAACAGCAACAAgg           >  2:2421720/1‑100 (MQ=255)
                                                                                          tACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGACACAACAGCAACAAgg           >  2:2421728/1‑100 (MQ=255)
                                                                                          tACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGACACAACAGCAACAAgg           <  1:1537925/100‑1 (MQ=255)
                                                                                          tACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGACACAACAGCAACAAgg           <  1:170959/100‑1 (MQ=255)
                                                                                           aCTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGACACAACAGCAACAAGGt          >  2:621441/1‑100 (MQ=255)
                                                                                           aCTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGACACAACAGCAACAAGGt          >  1:77038/1‑100 (MQ=255)
                                                                                               gATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGACACAACAGCAACAAGGTGAAg      <  2:251412/100‑1 (MQ=255)
                                                                                                aTGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGACACAACAGCAACAAGGTGAAGc     <  2:134371/100‑1 (MQ=255)
                                                                                                aTGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGACACAACAGCAACAAGGTGAAGc     <  1:1445155/100‑1 (MQ=255)
                                                                                                   cccAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGACACAACAGCAACAAGGTGAAGcaac  >  1:1828521/1‑100 (MQ=255)
                                                                                                   cccAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGACACAACAGCAACAAGGTGAAGcaac  >  1:1730093/1‑100 (MQ=255)
                                                                                                   cccAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGACACAACAGCAACAAGGTGAAGcaac  >  2:2045439/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGACACAACAGCAACAAGGTGAAGCAAC  >  NC_000913/236470‑236668

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: