Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,525,990 A→G 7.8% S651S (TCT→TCC)  dnaQ ← DNA‑directed DNA polymerase III epsilon subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,525,9900AG 7.8% 154.5 / 2.6 59S651S (TCT→TCC) dnaQDNA‑directed DNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (28/26);  new base G (3/2);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACATCACTCGGTATAAATACACTTGTATTTTTTGCACTTTGTAACGACGTGTTAACTTCAAATGTATTAAAGTGAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATACATCTTTTACAGCGTAATCTTTCACATATAAACGTATTGTTGATTTTTGAGATAAATTTTTAATCGAAATGAAACTAGTATG  >  CP000730/1525856‑1526121
                                                                                                                                      |                                                                                                                                   
aCATCACTCGGTATAAATACACTTGTATTTTTTGCACTTTGTAACGACGTGTTAACTTCAAATGTATTAAAGTGAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGAtata                                                                                                                               >  1:417173/1‑141 (MQ=255)
aCATCACTCGGTATAAATACACTTGTATTTTTTGCACTTTGTAACGACGTGTTAACTTCAAATGTATTAAAGTGAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGAtata                                                                                                                               <  2:339630/141‑1 (MQ=255)
  aTCACTTGGTATAAATACACTTGTATTTTTTGCACTTTGTAACGACGTGTTAACTTCAAATGTATTAAAGTGAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATaaa                                                                                                                             <  2:497491/141‑1 (MQ=255)
     aCTCAGTATTAATACACTTGTATTTTTTGCACTTTGTAACGACGTGTTAACTTCAAATGTATTAAAGTGAACATCTTGGTTGAATAACTGTTTAAAGGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAGAt                                                                                                                          <  2:172470/141‑1 (MQ=255)
      cTCGGTATAAATACACTTGTATTTTTTGCACTTTGTAACGACGTTTTAACTTCAAATGTATTAAAGTGAGCATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAATTTTTAAAGTTCCAAATAGaagagtt                                                                                                                         >  2:328509/1‑136 (MQ=255)
         ggTATAAATACACTTGTATTTTTTGCACTTTGTAACGACGTGTTAAATCCACATATGTTAAATTGAGCGGCTTTGTTGAATGACTGTCCGAATTCGCCAAACCTAGTATTAAATTTTAAATCACCAGATATAAAAAtcacc                                                                                                                      >  2:120233/1‑137 (MQ=255)
           tataAATACACTTGTATTTTTTGCACTTTGTAACGACGTGTTAACTTCAAATGTATTAAAGTGAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGa                                                                                                                    <  2:461363/141‑1 (MQ=255)
             taAATACACTTGTATTTTTTGCACTTTGTAACGACGTGTTAACTTCAAATGTATTAAAGTGAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGAtt                                                                                                                  <  2:495322/141‑1 (MQ=255)
              aaaTACACTTGTATTTTTTGCACTTTGTAACGACGTGTTAACTTCAAATGTATTAAAGTGAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGAttt                                                                                                                 <  1:40460/141‑1 (MQ=255)
                   cacTTGTATTTTTTGCACTTTGTAACGACGTGTTAACTTCAAATGTATTAAAGTGAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAAtt                                                                                                            <  1:460424/141‑1 (MQ=255)
                   cacTTGTATTTTTTGCACTTTGTAACGACGTGTTAACTTCAAATGTATTAAAGTGAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAAtt                                                                                                            >  2:134393/1‑141 (MQ=255)
                    acTTGTATTTTTTGCACTTTGTAACGACGTGTTAACTTCAAATGTATTAAAGTGAACATCTATGTTGAATAACTATTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAACTGAGCGATTTAAAttt                                                                                                           >  1:398286/1‑141 (MQ=255)
                         tATTTTTTGCACTTTGTAACGACGTGTTAACTTCAAATGTATTAAAGTGAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCa                                                                                                      >  1:20060/1‑141 (MQ=255)
                          aTTTTTTGCACTTTGTAACGACGTGTTAACTTCAAATGTATTAAGGTGAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCaa                                                                                                     >  1:323315/1‑141 (MQ=255)
                          aTTTTTTGCACTTTGTAACGACGTGTTAACTTCAAATGTATTAAAGTGAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCaa                                                                                                     >  1:561745/1‑141 (MQ=255)
                                 gCACTTTGTAACGACGTGTTAACTTCAAATGTATTAAAGTGAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGt                                                                                              >  1:359578/1‑141 (MQ=255)
                                 gCACTTTGTAACGACGTGTTAACTTCAAATGTATTAAAGTGAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGt                                                                                              <  1:514426/141‑1 (MQ=255)
                                 gCACTTTGTAACGACGTGTTAACTTCAAATGTATTAAAGTGAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGt                                                                                              <  1:229943/141‑1 (MQ=255)
                                  cACTTTGTAACGACGTGTTAACTTCAAATGTATTAAAGTGAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGtt                                                                                             <  1:198646/141‑1 (MQ=255)
                                      ttGTAACGACGTGTTAACTTCAAATGTATTAAAGTGAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGt                                                                                         >  1:288194/1‑141 (MQ=255)
                                           acgacgTGTTAACTTCAAATGTATTAAAGTGAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATa                                                                                    >  2:393579/1‑141 (MQ=255)
                                           acgacgTGTTAACTTCAAATGTATTAAAGTGAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTGAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATa                                                                                    <  1:311412/141‑1 (MQ=255)
                                             gacgTGTTAACTTCAAAAGATTTTGAGTGAACAGCTTTGTTGAATAACTGTGTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATACa                                                                                  >  2:405880/1‑141 (MQ=255)
                                               gggtttAACTTCAAATGTATTAAAGTGAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATACATc                                                                                <  2:569986/137‑1 (MQ=255)
                                                 tgtTAACTTCAAATGTATTAAAGTGAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATACATCtt                                                                              <  2:111095/141‑1 (MQ=255)
                                                    tAACTTCAAATGTATTAAAGTGAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAACAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATACATCTTTTa                                                                           <  2:20089/141‑1 (MQ=255)
                                                     aaCTTCAAATGTATTAAAGTGAAACACTTTGTTTAAGAACTGTTGAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATCCTGGTTTCGGTAATACATCTTTTAc                                                                          <  2:118284/141‑1 (MQ=255)
                                                          cAAATGTATTAAAGTGAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTTAATTTTTCCAATACTTGTTTCGTTAATACATCTTTTACAGCGt                                                                     <  1:536153/141‑1 (MQ=255)
                                                           aaaTGTATTAAAGTGAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCGGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATACATCTTTTACAGCGTa                                                                    >  2:273689/1‑141 (MQ=255)
                                                           aaaTGTATTAAAGTGAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCGGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATACATCTTTTACAGCGTa                                                                    >  2:273700/1‑141 (MQ=255)
                                                             aTGTATTAAAGGTAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTGAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGCTTCCGTACTTGATCTTTTACAGCGTAAt                                                                  <  2:549875/141‑1 (MQ=255)
                                                                 aTTAAAGTGAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATACATCTTTTACAGCGTAATCttt                                                              >  1:122507/1‑141 (MQ=255)
                                                                 aTTAAAGTGAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATACATCTTTTACAGCGTAATCttt                                                              >  1:140808/1‑141 (MQ=255)
                                                                    aaaGTGAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATACATCTTTTACAGCGTAATCTTTcac                                                           <  1:210292/141‑1 (MQ=255)
                                                                    aaaGTGAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATACATCTTTTACAGCGTAATCTTTcac                                                           <  1:566569/141‑1 (MQ=255)
                                                                          aaCATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCGGATATAAAAATAAGCGATTTAAATTTTTCAAATACTTGTTTCGTTAATACATCTTTTCCATCGTCATCTTTCCCATATaa                                                     >  2:527554/1‑141 (MQ=255)
                                                                                tttGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATACATCTTTTACAGCGTAATCTTTCACATATAAACGTAt                                               <  1:523516/141‑1 (MQ=255)
                                                                                 ttgtttAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATACATCTTTTACAGCGTAATCTTTCACATATAAACGTAtt                                              <  2:309950/141‑1 (MQ=255)
                                                                                  tgttgAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCATTAATACATCTTTTACAGCGTAATCTTTCACATATAAACGTAttt                                             >  1:515066/1‑140 (MQ=255)
                                                                                   gttgAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATACATCTTTTACAGCGTAATCTTTCCCATATAAACGTAttgt                                            >  2:219211/1‑141 (MQ=255)
                                                                                   gttgAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATACATCTTTTACAGCGTAATCTTTCACATATAAACGTAttgt                                            >  2:146793/1‑141 (MQ=255)
                                                                                            cTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTAGTTAATACATCTTTTACAGCGTAATCTTTCACATATAAACGTATTGTTGATTTTTg                                   >  1:501031/1‑141 (MQ=255)
                                                                                               tttAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATACATCTTTTACAGCGTAATCTTTCACATATAAACGTATTGTTGATTTTTgagc                                >  2:348537/1‑140 (MQ=255)
                                                                                                ttAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATACATCTTTTACAGCGTAATCTTTCACATATAAACGTATTGTTGATTTTTGAGAt                               >  2:364172/1‑141 (MQ=255)
                                                                                                    aaGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATACATCTTTTACAGCGTAATCTTTCACATATAAACGTATTGTTGATTTTTGAGATAAAt                           >  1:60336/1‑141 (MQ=255)
                                                                                                        ttCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATACATCTTTTACAGCGTAATCTTTCACATATAAACGTATTGTTGATTTTTGAGATAAAttttt                       >  2:18689/1‑141 (MQ=255)
                                                                                                         tCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATACATCTTTTACAGCGTAATCTTTCACATATAAACGTATTGTTGATTTTTGAGATAAATTTTTa                      >  2:46139/1‑141 (MQ=255)
                                                                                                         tCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATACATCTTTTACAGCGTAATCTTTCACATATAAACGTATTGTGGATTTTTGAGATAAATTTTTa                      >  1:12393/1‑141 (MQ=255)
                                                                                                          cAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATACATCTTTTACAGCGTAATCTTTCACATATAAACGTATTGTTGATTTTTGAGATAAATTTTTaa                     >  1:11620/1‑141 (MQ=255)
                                                                                                                attGATTAAATTTTTAAGTTCCAGATATAAAAATGAGCGATTTATTTTTTTCCTTTTCTTGTTTCGTTAATACATCTTTTACAGCGTAATCTTTCACATATAAACGTATTGTTGATTTTTGAGATAAATTTTTAATAGaaa               <  2:259343/139‑1 (MQ=255)
                                                                                                                aaTGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTATATTTTTCCAATACTTGTTTCGTTAATACATCTTTTACAGCGTAATCTTTCACATATAAACGTATTGTTGATTTTTGAGATAAATTTTTAATCGaaa               <  2:288194/141‑1 (MQ=255)
                                                                                                                 aTGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATACATCTTTTACAGCGTAATCTTTCACATATAAACGTATTGTTGATTTTTGAGATAAATTTTTAATCGAAAt              <  1:19546/141‑1 (MQ=255)
                                                                                                                 aTGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATACATCTTTTACAGCGTAATCTTTCACATATAAACGTATTGTTGATTTTTGAGATAAATTTTTAATCGAAAt              <  2:259738/141‑1 (MQ=255)
                                                                                                                 aTGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATACATCTTTTACAGCGTAATCTTTCACATATAAACGTATTGTTGATTTTTGAGATAAATTTTTAATCGAAAt              <  2:492019/141‑1 (MQ=255)
                                                                                                                 aTGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATACATCTTTTACAGCGTAATCTTTCACATATAAACGTATTGTTGATTTTTGAGATAAATTTTTAATCGAAAt              >  2:177618/1‑141 (MQ=255)
                                                                                                                   gATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATCCTTGTTTCGTTAATACATCTTTTACAGCGTAATCTTTCACATATAAACGTATTGTTGATTTTTGAGATAAATTTTTAATCGAAATGa            >  1:278646/1‑141 (MQ=255)
                                                                                                                    aTTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATACATCTTTTACAGCGTAATCTTTCACATATAAACGTATTGTTGATTTTTGAGATAAATTTTTAATTGAAATGaa           >  2:341449/1‑141 (MQ=255)
                                                                                                                             tAAAGTACCGGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATACATCTTTTACAGCGTAATCTTTCACATATAAACGTATTGTTGATTTTTGAGATAAATTTTTAATCGAAATGAAACTAGTATg  <  1:273700/141‑1 (MQ=255)
                                                                                                                             tAAAGTACCGGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATACATCTTTTACAGCGTAATCTTTCACATATAAACGTATTGTTGATTTTTGAGATAAATTTTTAATCGAAATGAAACTAGTATg  <  1:273689/141‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                   
ACATCACTCGGTATAAATACACTTGTATTTTTTGCACTTTGTAACGACGTGTTAACTTCAAATGTATTAAAGTGAACATCTTTGTTGAATAACTGTTTAAAAGCTTCAAACGAATGATTAAATTTTAAAGTACCAGATATAAAAATGAGCGATTTAAATTTTTCCAATACTTGTTTCGTTAATACATCTTTTACAGCGTAATCTTTCACATATAAACGTATTGTTGATTTTTGAGATAAATTTTTAATCGAAATGAAACTAGTATG  >  CP000730/1525856‑1526121

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