Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,745,317 C→T 7.6% P550P (CCG→CCA)  recJ ← recombination protein RecJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,745,3170CT 7.6% 149.4 / 3.6 53P550P (CCG→CCA) recJrecombination protein RecJ
Reads supporting (aligned to +/- strand):  ref base C (24/25);  new base T (2/2);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.89e-01

CGCATTTACTTTATCACTTTTAGGATGAATAAGCACGACAATATTTTCATCATTTTCTGTAAAAGGTAAACTTTTTCGCTTACTTCTATAATCTAATATTTGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATCCTAATATATTAATTGGCTGTTCATCTTGTAACTCAGGTTCTAAGTGTCCATTTTGCCAAAATAGCGCGGCGATATTACTTTCACCAAGTG  >  CP000730/1745179‑1745447
                                                                                                                                          |                                                                                                                                  
cGCATTTACTTTATCACTTTTAGGATGAATAAGCACGACAATATTTTCATCATTTTCTGTAAAAGGTAAACTTTTTCGCTTACTTCTATAATCTAATATTTGCTGTTCATTCACCGCAATATCTTGAATAATTATTTGCgg                                                                                                                                  <  2:473409/141‑1 (MQ=255)
   aTTTACTTTATCACTTTTAGGATGAATAAGCACGACAATATTTTCATCATTTTCTGTAAAAGGTAAACTTTTTCGCTTACTTCTATAATCTAATATTTGCTGTTCATTCATCGCAATATCTTGAATGATTATTTGCGGtta                                                                                                                               >  1:122449/1‑139 (MQ=255)
     ttACTTTATCACTTTTAGGATGAATAAGCACGACAATATTTTCATCATTTTCTGTAAAAGGTAAACTTTTTAGCTTACTTCTATAATCTAATATTTGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGCGAtgatt                                                                                                                             >  2:283894/1‑141 (MQ=255)
      tACTTTATCACTTTTAGGATGAATAAGCACGACAATATTTTCATCATTTTCTGTAAAAGGTAAACTTTTTCGCTTACTTCTATAATCTAATATTTGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGCGGtgattg                                                                                                                            >  2:119231/1‑141 (MQ=255)
           tattgcTTTTAGGATGAATAAGCACGACAATATTTTCATCATTTTCTGTAAAAGGTAAACTTTTTCGCTTACTTCTATAATCTAATATTTGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGCGGTGATTGATTAc                                                                                                                       <  2:329801/136‑1 (MQ=255)
                cTTTTAGGATGAATAAGCACGACAATATTTTCATCATTTTCTGTAAAAGGTAAACTTTTTCGCTTACTTCTATAATCTAATATTTGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGCGGTGATTGATTAccattc                                                                                                                  >  1:34385/1‑141 (MQ=255)
                cTTTTAGGATGAATAAGCACGACAATATTTTCATCATTTTCTGTAAAAGGTAAACTTTTTCGCTTACTTCTATAATCTAATATTTGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGCGGTGATTGATTAccattc                                                                                                                  >  1:385579/1‑141 (MQ=255)
                       gATGAATAAGCACGACGATATTTTCATCATTTTCTGTAAAAGGTAAACTTTTTCGCTTACTTCTATAATCTAATATTAGCTGTTCATTCATCGCAATCTCTAGAATAATTATTTGCGGTGTTTGATTACCCTTCCATTCGt                                                                                                           >  2:251588/1‑141 (MQ=255)
                       gATGAATAAGCACGACAATATTTTCATCATTTTCTGTAAAAGGTAAACTTTTTCGCTTACTTCTATAATCTAATATTTGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGt                                                                                                           >  1:606268/1‑141 (MQ=255)
                       gATGAATAAGCACGACAATATTTTCATCATTTTCTGTAAAAGGTAAACTTTTTCGCTTACTTCTATAATCTAATATTTGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGt                                                                                                           <  1:606476/141‑1 (MQ=255)
                            aTAAGCACGACAATATTTTCATCATTTTCTGTAAACGGTACCCTTTTTAGCTTTGTGCTGGAAGCTAATATTTGCTGTTCATTCATCGCAATATCGTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGAtt                                                                                                      <  2:140419/141‑1 (MQ=255)
                              aaGCACGACAATCTTTTCATCATTTTCTGTAAAAGGTAAACTTTTTCGCTTACTTCTATAATCTAATATTTGCTGTTCATTCTTCGAAATATCTTAAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTg                                                                                                    >  2:29412/1‑141 (MQ=255)
                              aaGCACGACAATATTTTCATCATTTTCTGTAAAAGGTAAACTTTTTCGCTTACTTCTATAATCTAATATTTGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTg                                                                                                    <  1:277876/141‑1 (MQ=255)
                               aGCACGACAATATTTTCATCATTTTCTGTAAAAGGTAAACTTTTTCGCTTACTTCTATAATCTAATATTTGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGa                                                                                                   <  1:582872/141‑1 (MQ=255)
                                 cACGACAATATTTTCATCATTTTCTGTAAAAGGTAAACTTTTTCGCTTACTTCTATAATCTAATATTTGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAAc                                                                                                 <  1:385186/141‑1 (MQ=255)
                                  aCGACAATATTTTCATCATTTTCTGTAAAAGGTAAACTTTTTCGCTTACTTCTATAATCTAATATTTGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAACa                                                                                                <  1:552691/141‑1 (MQ=255)
                                   cGACAATATTTTCATCATTTTCTGTAAAAGGTAAACTTTTTCGCTTACTTCTATAATCTAATATTTGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAACAg                                                                                               <  2:34385/141‑1 (MQ=255)
                                      aaaTATTTTCATCATTTTCTGTAAAAGGTAAACTTTTTCGCTTACTTCTATAATCTAATATTTGCTGTTCATTCATCGCAATATCTTGACTCATTATTTGCGGCGATTGAGTACCATCCCATTCGTTGATTTGAACAGATc                                                                                            <  1:91789/140‑1 (MQ=255)
                                       aaTGTTTTCATCATTTTCTGTAAAAGGTAAACTTTTTCGCTTACTTCTATAATCTAATATTTGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATcc                                                                                           <  1:283582/141‑1 (MQ=255)
                                       aaTATTTTCATCATTTTCTGTAAAAGGTAAACTTTTTCGCTTACTTCTATAATCTAATATTTGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATcc                                                                                           <  2:149118/141‑1 (MQ=255)
                                         tatTTTCATCATTTTCTGTAAAAGGTAAACTTTTTCGCTTACTTCTATAATCTAATATTTGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATCCTa                                                                                         >  1:419833/1‑141 (MQ=255)
                                                atcatTTTCTGTAAAAGGTAAACTTTTTCGCTTACTTCTATGATCTAATATTTGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATCCTAATATAtt                                                                                  <  2:91012/141‑1 (MQ=255)
                                                atcatTTTCTGTAAAAGGTAAACTTTTTCGCTTACTTCTATAATCTAATATTTGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATCCTAATATAtt                                                                                  <  1:268131/141‑1 (MQ=255)
                                                         tGTAAAAGGTAAACTTTTTCGCTTACTTCTATAATCTAATATTTGCTGTTCATTCATCGCAATCTCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATCCTAATATATTAATTGGCTg                                                                         <  1:450820/141‑1 (MQ=255)
                                                          gTAAAAGTTAAACTTTTTCGCTTACTTCTATAATCTAACATTAACTGTTCATTCATCGCAATATCTTGAATAATTATTGGCGGTGATTGATTACCATTCCACTCGTTGATTTGAACAGATCCTAATAGATTAATTGGCTGt                                                                        <  1:51990/141‑1 (MQ=255)
                                                                   aaaCTTTTTCGCTTACTTCTATAATCTAATATTTGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGTGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATCCTAATATATTAATTGGCTGTTCATCTTGt                                                               >  1:548681/1‑141 (MQ=255)
                                                                   aaaCTTTTTCGCTTACTTCTATAATCTAATATTTGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGTGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATCCTAATATATTAATTGGCTGTTCATCTTGt                                                               >  1:548705/1‑141 (MQ=255)
                                                                      cTTTTTCGCTTACTTCTATAATCTAATATTTGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGTGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATCCTAATATATTAATTGGCTGTTCATCTTGTAAc                                                            <  2:548681/141‑1 (MQ=255)
                                                                      cTTTTTCGCTTACTTCTATAATCTAATATTTGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGTGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATCCTAATATATTAATTGGCTGTTCATCTTGTAAc                                                            <  2:548705/141‑1 (MQ=255)
                                                                      cTTTTTCGCTTACTTCTATAATCTAATATTTGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTGCATTCGTTGATTTGAACAGATCCTAATATATTAATTGGCTGTTCATATTGTAAc                                                            >  2:70957/1‑141 (MQ=255)
                                                                       tttttCGCTTACTTCTATAATCTAATATTTGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATGCGTTGATTTGAACAGATCCTAATATGGCAAGTGGCTGTTCATCTTGTAACt                                                           <  2:474037/141‑1 (MQ=255)
                                                                        tttcCGATTACTTTTCTAATCATATATTTGTTGTTCATTCAGTGCAATTTTTTGAATAATTATTTGCGGTGATTTATTACCATTCCATTCGTTGATTTGAACAGATCCTAATATATTAATTGGCTGTTCATCTTGTAACTc                                                          <  2:129550/141‑1 (MQ=255)
                                                                          ttCGCTTACTTCTATAATCTAATATTTGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATCCTAATATATTAATTGGCTGTTCATCTTGTAACTCAg                                                        <  1:505805/141‑1 (MQ=255)
                                                                            cGCTTACTTCTATAATCTAATATTTGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATCCTAATATATTAATTGGCTGTTCATCTTGTAACTCAGGt                                                      >  2:30716/1‑141 (MQ=255)
                                                                                         aaTCTAATATTTGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATCCTAATATATTAATTGGCTGTTCATCTCGTAACTCAGGTTCTACGTGTCCAt                                         >  2:544206/1‑141 (MQ=255)
                                                                                             tAATATTTGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATCCTAATATATTAATTGGCTGTTCATCTTGTAACTCAGGTTCTAAGTGTCCATTTTg                                     >  1:136074/1‑141 (MQ=255)
                                                                                               attgttGCTGTTCATTCAGCGCAATATGTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATTCTAATATATTAATTGGCTGTTCATCTTGTAACTCAGGTTCTAAGTGTCCATTTTGcc                                   <  2:362730/137‑1 (MQ=255)
                                                                                                 atTTGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAAAAGATCCTAATATATTAATTGGCTGTTCATCTTGTAACGCAAGTTCTAAGTGTCCATCTTGCCaa                                 >  2:393890/1‑141 (MQ=255)
                                                                                                  tttGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATCCTAATATATTAATTGGCTGTTCATCTTGTCACTCAGGTTCTAAGTGTCCATTTTGCCaaa                                >  2:553899/1‑141 (MQ=255)
                                                                                                  tttGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATCCTAATATATTAATTGGCTGTTCATCTTGTAACTCAGGTTCTAAGTGTCCATTTTGCCaaa                                >  1:539821/1‑141 (MQ=255)
                                                                                                      cTGTTCATTCATCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATCCTAATATATTAATTGGCTGTTCATCTTGTAACTCAGGTTCTAAGTGTCCATTTTGCCAAAATAg                            >  1:526380/1‑141 (MQ=255)
                                                                                                       tGTTCTGTCATCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATCCTAATATATTAATTGGCTGTTCATCTTGTAACTCAGGTTCTAAGTGTCCATTTTGCCAAAATAgc                           <  2:666430/141‑1 (MQ=255)
                                                                                                       tGTTCATTCATCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATCCTAATATATTAATTGGCTGTTCATCTTGTAACTCAGGTTCTAAGTGTCCATTTTGCCAAAATAgc                           <  2:174650/141‑1 (MQ=255)
                                                                                                         ttcattcaTCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATCCTAATATATTAATTGGCTGTTCATCTTGTAACTCAGGTTCTAAGTGTCCATTTTGCCAAAATAgcgc                         <  1:543614/141‑1 (MQ=255)
                                                                                                           cattcatCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATCCTAATATATTAATTAGCTGTTCATCTTGTAACTCAGGTTCTAAGTGTCCATTTTGCCAAAATAGCgcgg                       >  2:119422/1‑141 (MQ=255)
                                                                                                               catCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATCCTAATATATTAATTGGCTGTTCATCTTGTAACTCAGGTTCTAAGTGTCCATTTTGCCAAAATAGCGCGGCGat                   >  1:247370/1‑141 (MQ=255)
                                                                                                               catCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATCCTAATATATTAATTGGCTGTTCAACTTGTAACTCAGTTTCTAAGGGTCCATTTTGCCAAAATAGCGCGGCGat                   >  1:452451/1‑141 (MQ=255)
                                                                                                                 tCGCAATATCTTGAATAATTATTTGCGGTGGTTGATTACCATTCCATTCGTTGCTTTGAACAGATCCTAATATATTAATTGGCTGTTCATCTTGTAACTCAGGTTCTAAGTGTCCATTTTGCCAAAATAGCGCGGCGatat                 >  1:625781/1‑141 (MQ=255)
                                                                                                                 tCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATCCTAATATATTAATTGGCTGTTCATCTTGTAACTCAGGTTCTAAGTGTCCATTTTGCCAAAATAGCGCGGCGatat                 >  1:234441/1‑141 (MQ=255)
                                                                                                                 tCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATCCTAATATATTAATTGGCTGTTCATCTTGTAACTCAGGTTCTAAGTGTCCATTTTGCCAAAATAGCGCGGCGatat                 >  2:218883/1‑141 (MQ=255)
                                                                                                                 tCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATCCTAATATATTAATTGGCTGTTCATCTTGTAACTCAGGTTCTAAGTGTCCATTTTGCCAAAATAGCGCGGCGatat                 >  1:234411/1‑141 (MQ=255)
                                                                                                                    cAATATCCTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATCCTAATATATTAATTGGCTGTTCATCTTGTAACTCAGGTTCTAAGTGTCCATTTTGCCAAAATAGCGCGGCGATATTAc              <  2:181838/141‑1 (MQ=255)
                                                                                                                                taatTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATCCTAATATATTAATTGGCTGTTCATCTTGTAACTCAGGTTCTAAGTGTCCATTTTGCCAAAATAGCGCGGCGATATTACTTTCACCAAgtg  <  2:483868/141‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                  
CGCATTTACTTTATCACTTTTAGGATGAATAAGCACGACAATATTTTCATCATTTTCTGTAAAAGGTAAACTTTTTCGCTTACTTCTATAATCTAATATTTGCTGTTCATTCATCGCAATATCTTGAATAATTATTTGCGGTGATTGATTACCATTCCATTCGTTGATTTGAACAGATCCTAATATATTAATTGGCTGTTCATCTTGTAACTCAGGTTCTAAGTGTCCATTTTGCCAAAATAGCGCGGCGATATTACTTTCACCAAGTG  >  CP000730/1745179‑1745447

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