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,1210GT95.7% 143.5 / ‑0.6 46intergenic (+205/‑67)rplA/rplJribosomal protein L1/ribosomal protein L10
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (30/14);  total (31/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 = 9.31e-01

ATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCGTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTATCAAATTCAGTTTCAACAGTA  >  CP000730/583973‑584265
                                                                                                                                                    |                                                                                                                                                  
atTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTg                                                                                                                                                  >  2:83513/1‑151 (MQ=255)
  tGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTggg                                                                                                                                                >  1:601644/1‑151 (MQ=255)
   gACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCGTGGGt                                                                                                                                               <  2:642298/151‑1 (MQ=255)
          aGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGt                                                                                                                                        <  2:965124/151‑1 (MQ=255)
             aaTTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCt                                                                                                                                     <  2:432115/151‑1 (MQ=255)
                        tataTTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTaa                                                                                                                          >  1:815794/1‑151 (MQ=255)
                             tACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTa                                                                                                                    >  2:23877/1‑151 (MQ=255)
                                  aTTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTaaaaaa                                                                                                               >  1:993432/1‑151 (MQ=255)
                                     gtgATTATTTTACCTAAGACAGTAGGAGTTATTTAGAACTTAAAATTTTGCCTGCCTAGGCTAAAATTTACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAGGGGCTTTTTTTATTAATTTTAAAAAAAGc                                                                                                            >  2:635927/1‑151 (MQ=255)
                                              ttACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTAAAATTAATAAAAAAAGCActttctac                                                                                                   <  1:468520/151‑8 (MQ=38)
                                                aCCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAAttt                                                                                                 >  1:751415/1‑151 (MQ=255)
                                                         aGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAgg                                                                                        <  1:302345/151‑1 (MQ=255)
                                                                gTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTTTTTTTTTTTTAAAGTGCTTTTTTTATTAATTTTAAAAAAATCATCAATTATTTAAATGGAGGTGTCTGa                                                                                 >  2:466054/1‑151 (MQ=255)
                                                                  tatttatAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTTTTTTTTTTTTAAAGTGCTTTTTTTATTCATTTTAATAAATGCACCAAAAATTTAAATGGTGGTGTCTGCAt                                                                               >  2:211796/1‑151 (MQ=255)
                                                                  tatttatAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAAt                                                                               >  1:853071/1‑151 (MQ=255)
                                                                    tttatAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGt                                                                             >  2:613999/1‑151 (MQ=255)
                                                                        taACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTTCTTTTTTTTTAAAGTGCTTTTTTTATTAATTTTAAAATAAGCACCAAAAATTTAAATGGAGGTGTCTGCCTGTCTGc                                                                         >  2:394297/1‑151 (MQ=255)
                                                                          aCTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTa                                                                       >  2:217920/1‑151 (MQ=255)
                                                                           cTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTAt                                                                      >  2:710087/1‑151 (MQ=255)
                                                                            ttAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATc                                                                     <  1:905288/151‑1 (MQ=255)
                                                                               aaaTTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCAtt                                                                  <  1:23877/151‑1 (MQ=255)
                                                                                aaTTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTGGGGTAGAAAGTGCTTTTTTTATTAATTTTACAAAAAGCACCAAAAATTCAAATTGAGGGGTCTGAATGTCTGCTATCAttt                                                                 >  2:621747/1‑150 (MQ=255)
                                                                                aaTTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCGTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGTAGGTGTCTGAATGTCTGCTATCATTg                                                                 >  2:344821/1‑151 (MQ=255)
                                                                                 aTTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGa                                                                >  1:353050/1‑151 (MQ=255)
                                                                                  tttATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGaa                                                               >  2:380343/1‑151 (MQ=255)
                                                                                         tGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTaaaa                                                       >  2:585815/1‑151 (MQ=255)
                                                                                          gCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTTTTTTTTTTAAAAAAAGCACCAAAAATTTAAATGTAGTTGTCTGAATGTCTGCTATCATTGAAGCTaaaaa                                                      >  2:363101/1‑151 (MQ=255)
                                                                                             gAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAaaca                                                   <  2:993432/151‑1 (MQ=255)
                                                                                               ggCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAaacaac                                                 <  2:663756/151‑1 (MQ=255)
                                                                                               ggCTAAAATGGAATTGAACGTGATGATATATGATCTTTCAAGCAATTTTTGCATTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAATTGGAGGTGTCTGATTTTCTTCTATCATTAAAGCTAAAAaacaac                                                 <  2:96465/151‑1 (MQ=255)
                                                                                                 cTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAAACAACt                                                >  1:212863/1‑151 (MQ=255)
                                                                                                  tAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAg                                              >  1:868162/1‑151 (MQ=255)
                                                                                                           cTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGaaa                                     >  1:613158/1‑151 (MQ=255)
                                                                                                               aaCGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGc                                 <  2:650608/151‑1 (MQ=255)
                                                                                                                          ctatgatctTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTATTTTATTAATTTTAAAAAAAGCACAAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAAATAGATGATGAAATTGCTGAGGTACTAt                      >  1:496676/1‑151 (MQ=255)
                                                                                                                           tatgatctTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTATc                     >  1:608753/1‑151 (MQ=255)
                                                                                                                              gatctTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCaaaaa                                                                                                    >  2:1020356/1‑70 (MQ=255)
                                                                                                                              gatctTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCaaaaa                                                                                                    <  1:1020356/70‑1 (MQ=255)
                                                                                                                                  tttCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTATTTTTATAAAAAGCACTAAAATTTTAAATTGTGGTGTTTGAATGTCTGCTATCATTTAAGCTAAAAAACAACTAGTTGATTAAATTGCTGAGGTACTATCAAATTCa              >  2:282695/1‑151 (MQ=255)
                                                                                                                                  tttCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGc                                 >  1:208630/1‑132 (MQ=255)
                                                                                                                                  tttCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGc                                 <  2:208630/132‑1 (MQ=255)
                                                                                                                                     cAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTATCAAATTCAGtt           >  1:393810/1‑151 (MQ=255)
                                                                                                                                     cAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTATCAAATTCAGtt           <  2:815794/151‑1 (MQ=255)
                                                                                                                                       aGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTATCAAATTCAGTTTc         >  1:746255/1‑151 (MQ=255)
                                                                                                                                          aCTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTATCAAATTCAGTTTCAAc      >  1:889380/1‑151 (MQ=255)
                                                                                                                                              tttGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTATCAAATTCAGTTTCAACAGTa  <  1:411156/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCGTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTATCAAATTCAGTTTCAACAGTA  >  CP000730/583973‑584265

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: