Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 584,121 G→T intergenic (+205/‑67) rplA → / → rplJ ribosomal protein L1/ribosomal protein L10

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730584,1210GT98.1% 162.9 / NA 52intergenic (+205/‑67)rplA/rplJribosomal protein L1/ribosomal protein L10
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (31/20);  total (31/21)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCGTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTATCAAATTCAGTT  >  CP000730/583990‑584256
                                                                                                                                   |                                                                                                                                       
tCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTGTTTGCCTTGGGTAGaa                                                                                                                                <  1:817938/141‑1 (MQ=255)
   aaGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGt                                                                                                                             <  1:153534/141‑1 (MQ=255)
    aGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAACATTCACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTg                                                                                                                            >  1:1354738/1‑141 (MQ=255)
    aGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTg                                                                                                                            >  1:726647/1‑141 (MQ=255)
    aGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTg                                                                                                                            >  1:1135501/1‑141 (MQ=255)
    aGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTg                                                                                                                            >  1:1229126/1‑141 (MQ=255)
         tatTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCtttt                                                                                                                       <  1:1206992/141‑1 (MQ=255)
         tatTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCtttt                                                                                                                       >  1:942151/1‑141 (MQ=255)
           ttACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCtttttt                                                                                                                     >  1:1051522/1‑141 (MQ=255)
            tACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGTTGAGCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCttttttt                                                                                                                    >  1:1293982/1‑141 (MQ=255)
            tACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTGGGGAAGAAAGGGCttttttt                                                                                                                    >  1:423729/1‑141 (MQ=255)
                   tgtgATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAAtt                                                                                                             <  1:1252195/141‑1 (MQ=255)
                          attTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCTAGGCTAGAATTGACTCGGACGTCATCATGTATGATCTTTCACGCACTATTTGCCTAGGGTAGAAAGTGCTCTTTTTATTAATTTTaaaaa                                                                                                      <  1:1028846/141‑1 (MQ=255)
                             ttACCTAAGACAGTGGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAg                                                                                                   >  1:898068/1‑141 (MQ=255)
                             ttACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAg                                                                                                   >  1:483396/1‑141 (MQ=255)
                             ttACCTAAGACAGTAGGAGGTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGCACGTGATGATCTATGATCGTTCAAGCACTTTGTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAg                                                                                                   >  1:729609/1‑141 (MQ=255)
                                    aGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCaaa                                                                                            >  1:1335449/1‑141 (MQ=255)
                                        aGTAGGAGTTATTTATAACTTAAAATTTGTCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAAtt                                                                                        >  1:888922/1‑141 (MQ=255)
                                                    ttatAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAACATTTAAATGGAGgtg                                                                            >  1:1260220/1‑141 (MQ=255)
                                                            tAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAAGTTAAATGGAGGTGTCGGaagg                                                                    >  1:48066/1‑139 (MQ=255)
                                                            tAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAAGTTGTAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATg                                                                    >  1:106518/1‑141 (MQ=255)
                                                              aaaTTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTc                                                                  <  1:458578/141‑1 (MQ=255)
                                                               aaTTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAACTGGAGGTGTCTGAATGTCt                                                                 <  1:821189/141‑1 (MQ=255)
                                                               aaTTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCt                                                                 <  1:197988/141‑1 (MQ=255)
                                                                       cTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCAt                                                         >  1:582004/1‑141 (MQ=255)
                                                                         gCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTg                                                       >  1:1119496/1‑141 (MQ=255)
                                                                                cTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTaa                                                >  1:823263/1‑141 (MQ=255)
                                                                                 tAAAATTGACTTGAACGGGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTaaa                                               <  1:1235861/141‑1 (MQ=255)
                                                                                        gACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAaacaac                                        >  1:199368/1‑141 (MQ=255)
                                                                                             gAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGtt                                   <  1:1238556/141‑1 (MQ=255)
                                                                                                   gatgatCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGaa                             >  1:102259/1‑141 (MQ=255)
                                                                                                   gatgatCTATGATCTTCCAAGCACTTTTTGCCTTGGGTGGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGaa                             >  1:764025/1‑141 (MQ=255)
                                                                                                    atgatctatgatctTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGATGTAATTGAAGCTAAAAAACAACTAGTTGATGaaa                            >  1:519542/1‑141 (MQ=255)
                                                                                                     tgatctatgatctTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACACCTAGTTGATGAAAt                           <  1:734167/141‑1 (MQ=255)
                                                                                                      gatctatgatctTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTCAAATGGAGGAGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTAGATGAAAtt                          >  1:1123838/1‑141 (MQ=255)
                                                                                                       atctatgatctTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTg                         <  1:800993/141‑1 (MQ=255)
                                                                                                                    cAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTATc            <  1:449449/141‑1 (MQ=255)
                                                                                                                    cAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTATc            >  1:1129826/1‑141 (MQ=255)
                                                                                                                    cAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTATc            <  1:24093/141‑1 (MQ=255)
                                                                                                                    cAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCAACAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTATc            <  1:1141173/141‑1 (MQ=255)
                                                                                                                     aaGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTATCa           >  1:671565/1‑141 (MQ=255)
                                                                                                                      aGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTATCaa          >  1:467541/1‑141 (MQ=255)
                                                                                                                      aGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTATCaa          <  1:24946/141‑1 (MQ=255)
                                                                                                                      aGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTATCaa          >  1:886209/1‑141 (MQ=255)
                                                                                                                      aGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTATCaa          <  1:151565/141‑1 (MQ=255)
                                                                                                                      aGCACTCTTTGACGTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTACAAAACAACTAGTTGATGAAATTGCTGAGGTACTATCaa          <  1:1229033/141‑1 (MQ=255)
                                                                                                                         aCTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTATCAAAtt       <  1:1190682/141‑1 (MQ=255)
                                                                                                                         aCTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTATCAAAtt       <  1:892567/141‑1 (MQ=255)
                                                                                                                         aCTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTAGTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTATCTAAtt       <  1:936356/141‑1 (MQ=255)
                                                                                                                              ttGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTCTCATTGAAGGAATAAAACAACGAGTTGATGAAATTGCTGAGGTACTATCAAATTCAGtt  >  1:445746/1‑141 (MQ=255)
                                                                                                                              ttGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTGTCAAATTCAGtt  >  1:62243/1‑141 (MQ=255)
                                                                                                                              ttGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTATCAAATTCAGtt  >  1:1134017/1‑141 (MQ=255)
                                                                                                                                   |                                                                                                                                       
TCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCGTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTATCAAATTCAGTT  >  CP000730/583990‑584256

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: