Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 593,430 T→A 12.2% intergenic (+22/‑58) rpoC → / → USA300HOU_0538 DNA‑directed RNA polymerase subunit beta'/possible RNA‑binding ribosomal protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730593,4300TA12.2% 86.2 / 5.2 33intergenic (+22/‑58)rpoC/USA300HOU_0538DNA‑directed RNA polymerase subunit beta'/possible RNA‑binding ribosomal protein
Reads supporting (aligned to +/- strand):  ref base T (16/13);  new base A (2/2);  total (18/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

TATTGGTAAGTTAATTCCAGCTGGTACTGGTATGAGACGTTATAGCGACGTAAAATACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAAC‑AGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATGCAAGCAGAACTTTGGAGGATAAATTATTGTCTAAGGAAAAAGTTGCA  >  CP000730/593294‑593565
                                                                                                                                         |                                                                                                                                       
tattGGTAAGTTAATTCCAGCTGGTACTGGTATGAGACGTTATAGCGACGTAAAATACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAAC‑AGAGGCTAATGCtt                                                                                                                                     <  2:253657/141‑1 (MQ=255)
tattGGTAAGTTAATTCCAGCTGGTACTGGTATGAGACGTTATAGCGACATAAAAAAAAAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAAC‑AGAGGCTAATGCtt                                                                                                                                     <  2:143394/141‑1 (MQ=255)
 atggaaaaGTTAATTCAAGTTGGTACTGGTAGGAGATGTTATAGAGAAGAAAAAAAAAAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAAGAGGCTAAAGcat                                                                                                                                     <  2:206832/135‑3 (MQ=255)
     gTAAGTTAATTCCAGCTGGTACTGGTATGAGACGTTATAGCGACGTAAAATACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAAC‑AGAGGCTAATGCTTTAGcc                                                                                                                                <  2:56033/141‑1 (MQ=255)
                ccAGCTGGTACTGGTATGAGACGTTATAGCGACGTAAAATACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAAC‑AGAGGCTAATGCTTTAGCCTCTTGTTAttt                                                                                                                     <  1:183982/141‑1 (MQ=255)
                   gCTGGTACTGGTATGAGACGTTATAGCGACGTAAAATACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAAC‑AGAGGCTAATGCTTTAGCCTCTTGTTATTTTTa                                                                                                                  >  1:192710/1‑141 (MQ=255)
                   gCTGGTACTGGTATGAGACGTTATAGCGACGTAAAATACGAAAAAACAGCTAAACCAATTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAAC‑AGAGGCTAATGCTTTAGCCTCTTGTTATTTTTa                                                                                                                  >  1:172323/1‑141 (MQ=255)
                           tggtaTGAGACGTTATAGCGACGTAAAATACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAAC‑AGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAatt                                                                                                          >  1:83069/1‑141 (MQ=255)
                             gtaTGAGACGTTATAGCGACGTAAAATACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAAC‑AGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAattat                                                                                                        <  2:333474/141‑1 (MQ=255)
                                         ataGCGACGTAAAATACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAAC‑AGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGt                                                                                            >  2:166762/1‑141 (MQ=255)
                                            gCGACGTAAAATACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAAC‑AGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGa                                                                                         >  2:37551/1‑141 (MQ=255)
                                              gACGTAAAATACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAAC‑AGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACg                                                                                       <  1:125417/141‑1 (MQ=255)
                                               aCGTAAAATACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAAC‑AGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGa                                                                                      >  2:100061/1‑141 (MQ=255)
                                                    aaaTACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAAC‑AGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATtct                                                                                 <  1:129087/141‑1 (MQ=255)
                                                      aTACGAAAAAACAGCTAACCCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAAC‑AGAGGCTAATGCTTTAGGCTCTTGTTATTTTTATGTAAATTATTTGCTTTAATGTTGACGAATTctct                                                                               >  2:285303/1‑141 (MQ=255)
                                                              aaaCAGCTAAACCAGTTGCAGTATTTTGCTCTCAATAGATAGTTATCGAAGAACAATTTTATTGC‑GGAGTCGAGTGCTTTATCCTCTTTTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAAt                                                                       <  1:47955/141‑1 (MQ=255)
                                                                   gCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAAC‑AGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTaa                                                                  >  1:307183/1‑141 (MQ=255)
                                                                    cTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAAC‑AGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAat                                                                 >  1:310734/1‑141 (MQ=255)
                                                                               gCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAAGAGGCTAAAGCATTAGCctc                                                                                                                              <  1:89473/70‑1 (MQ=25)
                                                                               gCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAAGAGGCTAAAGCATTAGCctc                                                                                                                              >  2:89473/1‑70 (MQ=25)
                                                                                     gTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAAC‑AGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATGc                                                >  2:132597/1‑141 (MQ=255)
                                                                                      ttGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAAC‑AGAGGCTAATCCTTTAGCCGCATGTTATTCTTATGTAATTTATTTGTTTTAATGTTGACGAATTCTCATGTTGCATGTTAATATATTAAAGGTTGATGCa                                               >  2:145482/1‑141 (MQ=255)
                                                                                                                tAACAAGTATATAAC‑AGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATGCAAGCAGAACTTTGGAGGATAAattatt                     >  2:163004/1‑141 (MQ=255)
                                                                                                                 aaCAAGTATATAAC‑AGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATGCAAGCAGAACTTTGGAGGATAAATTATTg                    >  2:240059/1‑141 (MQ=255)
                                                                                                                 aaCAAGTATATAAC‑AGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTAGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATGCAAGCAGAACTTTGGAGGATAAATTATTg                    >  1:162973/1‑141 (MQ=255)
                                                                                                                  aCAAGTATATAAC‑AGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATGCAAGCAGAACTTTGGAGGATAAATTATTGt                   >  1:155936/1‑141 (MQ=255)
                                                                                                                    aaGTATATAAC‑AGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATGCAAGCAGAACTTTGGAGGATAAATTATTGTCt                 <  1:74613/141‑1 (MQ=255)
                                                                                                                      gTATATAAC‑AGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGTTTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATGCAAGCAGAACTGTGGAGGATAAATTATTGTTTaa               >  1:173256/1‑141 (MQ=255)
                                                                                                                       tatttaGC‑AGAGGCTAATGGTTTGGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTGCATGTTAGTATATTAAAGGTTGGGGCAAGCAGAACTTTGGAGGATGCATTATTGTCTAAg              <  2:125190/141‑1 (MQ=255)
                                                                                                                        atataAC‑AGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATGCAAGCAGAACTTTGGAGGATAAATTATTGTCTAAgg             <  2:201149/141‑1 (MQ=255)
                                                                                                                             acaagagGCTAAAGCATGAGCCTCTTGTTATATTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATGCAAGCAGAACTTTGGAGGATAAATTATTGACTagtggata         >  2:94376/4‑134 (MQ=255)
                                                                                                                             aC‑AGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATGCAAGCAGAACTTTGGAGCATAAATTATTGTCTAAGGaaaaa        <  2:157986/141‑1 (MQ=255)
                                                                                                                                    tcTAATGCTTTAGCCTGTTGTTATTTTTATTTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATGCAAGCAGAACTTTGGAGGATAAATTATTGTCTAAGGAAAAAGTTGCa  <  2:108851/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                       
TATTGGTAAGTTAATTCCAGCTGGTACTGGTATGAGACGTTATAGCGACGTAAAATACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAAC‑AGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATGCAAGCAGAACTTTGGAGGATAAATTATTGTCTAAGGAAAAAGTTGCA  >  CP000730/593294‑593565

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