Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 401,596 A→G 10.3% K47K (AAA→AAG USA300HOU_0370 → possible transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730401,5960AG10.3% 116.6 / 3.2 47K47K (AAA→AAGUSA300HOU_0370possible transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (19/23);  new base G (3/2);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGTAATCGATTGAAAGAATTACGAGCACGAGATGGCTTAAACCAAACGCAACTTGCTAAACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGTGGACTTGTAGGGGGTATCATAGGTGGAATTTTAGGTAATTT  >  CP000730/401457‑401735
                                                                                                                                           |                                                                                                                                           
tgCGTAATCGATTGAAAGAATTACGAGCACGAGATGGCTTAAACCAAACGCAACTTGCTAAACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATaaaaa                                                                                                                                            >  2:171917/1‑141 (MQ=255)
 gCGTAATCGATTGAAAGAATTACGAGCACGAGATGGCTTAAACCAAACGCAACTTGCTAAACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAAt                                                                                                                                           >  2:167699/1‑141 (MQ=255)
 gCGTAATCGATTGAAAGAATTACGAGCACGAGATGGCTTAAACCAAACGCAACTTGCTAAACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAAt                                                                                                                                           >  2:408292/1‑141 (MQ=255)
     aaTCGATTGAAAGAATTACGAGCACGAGATGGCTTAAACCAAACGCAACTTGCTAAACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACACTTTTATGCCATCAGTATTAACGGCAATAAAAATTGCt                                                                                                                                       >  2:260481/1‑141 (MQ=255)
       tCGATTGAAAGAATTACGAGCACGAGATGGCTTAAACCAAACGCAACTTGCTAAACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCg                                                                                                                                     <  2:60011/141‑1 (MQ=255)
                   gttACGAGCACGAGATGGCTTAAATCAAACGCAGCTTGCCAAACTAGCTGGTGTTTCAAGACAAACCATTTCGCTAATTGAGCGAAACAATTTTATGCCATCACTATTAACGGCAATAAAGATTGCTCGCATTTTCAATGa                                                                                                                         >  2:413761/2‑141 (MQ=255)
                       cGAGCACGAGATGGCTTAAACCAAACGCAACTTGCTAAACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACg                                                                                                                     >  1:223408/1‑141 (MQ=255)
                       cGAGCACGAGATGGCTTAAACCAAACGCAACTTGCTAAACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACg                                                                                                                     >  1:54759/1‑141 (MQ=255)
                           cACGAGATGGCTTAAACCAAACGCAACTTGCTAAACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGTCAATAAAAATTGCTCGCATTTTCAATGAAACGGTgg                                                                                                                 >  2:308343/1‑141 (MQ=255)
                               agaTGGCTTAAACCAAACGCAACTTGCTAAACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAAc                                                                                                             >  2:203681/1‑141 (MQ=255)
                                 aTGGCTTAAACCAAACGCAACTTGCTAAACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTg                                                                                                           <  1:374419/141‑1 (MQ=255)
                                       tAAACCAAACGCAACTTGCTAAACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTtta                                                                                                     <  1:279326/141‑1 (MQ=255)
                                       tAAACCAAACGCAACTTGCTAAACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTAGAAACTGTTTtta                                                                                                     >  2:62302/1‑141 (MQ=255)
                                        aaaCCAAACGCAACTTGCTAAACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTttat                                                                                                    <  2:103806/141‑1 (MQ=255)
                                          atcAAACGCAGCTTGCCAAACTAGCTGGTGTTTCAAGACAAACCATTTCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAGATTGCTCGCATTTTCAATGAAACGGTGGAAACGGTTTttatta                                                                                                  <  1:39702/139‑1 (MQ=255)
                                           ccAAACGCAACTTGCTAAACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTtattat                                                                                                 <  1:103861/141‑1 (MQ=255)
                                              aaCGCAACTTGCTAAACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTga                                                                                              <  2:261695/141‑1 (MQ=255)
                                                cGCAACTTGCTACACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTgagg                                                                                            <  2:366581/141‑1 (MQ=255)
                                                    aCTTGCTAAACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAggagga                                                                                        >  1:188210/1‑141 (MQ=255)
                                                        gCTAAACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTTATCGATACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAACCGGGGGAAACTGTTTTTATTATTGAGGAGGATGAg                                                                                    >  2:386808/1‑141 (MQ=255)
                                                        gCTAAACAAGCGGGCGTTCCAAGACAAACCATATCGCTAATTGATCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAg                                                                                    >  2:134054/1‑141 (MQ=255)
                                                        gCTAAACAAGCGGGCGTTACAAGACAAACCATATCGCTACTTGAGCGAAACAATTTTATGCCATCGGTATTTACGACAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAg                                                                                    <  1:386641/141‑1 (MQ=255)
                                                            aaCATGCGGGCCTACAAAAACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAGCGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCAt                                                                                <  2:116309/141‑1 (MQ=255)
                                                              cAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGa                                                                              <  2:59261/141‑1 (MQ=255)
                                                                  cGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAAt                                                                          >  2:246349/1‑141 (MQ=255)
                                                                   gggCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATa                                                                         >  2:106301/1‑141 (MQ=255)
                                                                        tttCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGAATGAGGCATGAAAATACTAAg                                                                    >  1:241879/1‑141 (MQ=255)
                                                                               aCAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATc                                                             >  1:306128/1‑141 (MQ=255)
                                                                                  aaCCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGa                                                          <  1:292928/141‑1 (MQ=255)
                                                                                    ccATATCGCTAATTGAGCGCAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAACCTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGata                                                        <  1:7697/141‑1 (MQ=255)
                                                                                      atatCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATc                                                      <  2:185775/141‑1 (MQ=255)
                                                                                       tatCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCt                                                     <  1:188822/141‑1 (MQ=255)
                                                                                        atCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGatatccgatatatctt                                                    <  1:356409/141‑13 (MQ=255)
                                                                                        atCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCtt                                                    <  1:167919/141‑1 (MQ=255)
                                                                                           gctcaattaGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTa                                                 >  2:286487/9‑141 (MQ=255)
                                                                                           gCTAATTGAGAGAAACAATTTTATGCCATCAGTATAAAAAAAAAAAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTa                                                 <  2:260655/141‑1 (MQ=255)
                                                                                            cTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTAc                                                <  1:167699/141‑1 (MQ=255)
                                                                                            cTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTAc                                                <  1:408292/141‑1 (MQ=255)
                                                                                                     cGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGTGGAc                                       <  2:371862/141‑1 (MQ=255)
                                                                                                         aCAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGCGGACTTGt                                   >  1:61711/1‑141 (MQ=255)
                                                                                                               ttATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGTGGACTTGTAggggg                             >  2:80985/1‑141 (MQ=255)
                                                                                                                tATGCCATCAGTATTAACGGCAATAAAGATTGCTCGCATTTTCAATGAAACGGTGGAAACGGTTTTTATTATTGAGGAGGATGAAGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGTGGAATTGTAGGGGGt                            >  2:337179/1‑141 (MQ=255)
                                                                                                                         aGTATTAACGGCAATAAAGATTGCTCGCATTTTCAATGAAACGGTGGAAACGGTTTTTATTATTGAGGAGGATGAAGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGTGGAATTGTAGGGGGTATCATTGGt                   <  1:419243/141‑1 (MQ=255)
                                                                                                                         aGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGTGGACTTGTAGGGGGTATCATAGGt                   <  1:25801/141‑1 (MQ=255)
                                                                                                                         aGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGTGGACTTGTAGGGGGTATCATAGGt                   <  1:25278/141‑1 (MQ=255)
                                                                                                                           tATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGTGGACTTGTAGGGGGTATCATAGGTgg                 <  2:256725/141‑1 (MQ=255)
                                                                                                                               aaCGGCAATAAAGATTGCTCGCATTTTCAATGAAACGGTGGAAACGGTTTTTATTATTGAGGAGGATGAAGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGTGGAATTGTAGGGGGTATCATTGGTGGAAtt             >  1:132/1‑141 (MQ=255)
                                                                                                                                      aTAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGTGGACTTGTAGGGGGTATCATAGGTGGAATTTTAGGTa      >  1:438089/1‑141 (MQ=255)
                                                                                                                                          aaaTTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGTGGACTTGTAGGGGGTATCATAGGTGGAATTTTAGGTAAttt  >  1:318580/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                           
TGCGTAATCGATTGAAAGAATTACGAGCACGAGATGGCTTAAACCAAACGCAACTTGCTAAACAAGCGGGCGTTTCAAGACAAACCATATCGCTAATTGAGCGAAACAATTTTATGCCATCAGTATTAACGGCAATAAAAATTGCTCGCATTTTCAATGAAACGGTGGAAACTGTTTTTATTATTGAGGAGGATGAGGCATGAAAATACTAAGATATATCGGATATCTTTTACTGGGTGGACTTGTAGGGGGTATCATAGGTGGAATTTTAGGTAATTT  >  CP000730/401457‑401735

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