Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 546,499 A→G 89.8% intergenic (+11/‑6) rplB → / → rpsS 50S ribosomal protein L2/30S ribosomal protein S19

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB143546,4990AG89.8% 113.9 / 4.3 40intergenic (+11/‑6)rplB/rpsS50S ribosomal protein L2/30S ribosomal protein S19
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base G (24/11);  total (26/14)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.67e-01

GTGGTGGTGAAGGTCGTACCTCCGGTGGTCGTCATCCGGTATCGCCATGGGGCTTCCCAACCAAGGGTGCTAAAACCCGTGGTAATAAGCGTACCGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATACGACAGTGCCACGTTCTCTGAAAAAAGGTCCTTTTATCGATCTTCACCTGTTGAAGAAGGTCGAAGTGGCGGTGGAGAAGAACGATCGCAAGCCAGTTAAAACCTGGTCGCGCCGTTCGATGATCCTGCCA  >  ACB143/546364‑546630
                                                                                                                                       |                                                                                                                                   
gtggtggtgAAGGTCGTACCTCCGGTGGGCGTCATCCGGTATCGCCATGGGGCTTCCCAACCAAGGGTGCTAAAACCCGTGGTAATAAGCGTACCGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATa                                                                                                                                     <  1:195609/136‑1 (MQ=255)
   gtggtgAAGGTCGTACCTCCGGTGGTCGTCATCCGGTATCGCCATGGGGCTTCCCAACCAAGGGTGCTAAAACCCGTGGTAATAAGCGTACCGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATGCGa                                                                                                                                  >  4:280412/1‑136 (MQ=255)
      gtgAAGGTCGTACCTCCGGTGGTCGTCATCCGGTATCGCCATGGGGCTTCCCAACCAAGGGTGCTAAAACCCGTGGTAATAAGCGTACCGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATGCGACAg                                                                                                                               >  1:170853/1‑136 (MQ=255)
      gtgAAGGTCGTACCTCCGGTGGTCGTCATCCGGTATCGCCATGGGGCTTCCCAACCAAGGGTGCTAAAACCCGTGGTAATAAGCGTACCGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATGCGACAg                                                                                                                               >  4:251195/1‑136 (MQ=255)
              cGTACCTCCGGTGGTCGTCATCCGGTATCGCCATGGGGCTTCCCAACCAAGGGTGCTAAAACCCGTGGTAATAAGCGTACCGACAACATGATCGTCCGTCGTCGCAAGTATCTAGAGGGATGCGACAGTGCCACGt                                                                                                                       >  2:140250/1‑136 (MQ=255)
                              gtcATCCGGTATCGCCATGGGGCTTCCCAACCAAGGGTGCTAAAACCCGTGGTAATAAGCGTACCGAGAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATGCGACAGTGCCACGTTCTCTGAAAAAAGGTc                                                                                                       >  2:258786/1‑136 (MQ=255)
                                  tCCGGTATCGCCATGGGGCTTCCCAACCAAGGGTGCTAAAACCCGTGGTAATAAGCGTACCGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATGCGACAGTGCCACGTTCTCTGAAAAAAGGTCCttt                                                                                                   <  8:154941/136‑1 (MQ=255)
                                  tCCGGTATCGCCATGGGGCTTCCCAACCAAGGGTGCTAAAACCCGTGGTAATAAGCGTACCGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATGCGACAGTGCCACGTTCTCTGAAAAAAGGTCCttt                                                                                                   <  8:229704/136‑1 (MQ=255)
                                  tCCGGTATCGCCATGGGGCTTCCCAACCAAGGGTGCTAAAACCCGTGGTAATAAGCGTACCGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATGCGACAGTGCCACGTTCTCTGAAAAAAGGTCCttt                                                                                                   <  3:280412/136‑1 (MQ=255)
                                   ccGGTATCGCCATGGGGCTTCCCAACCAAGGGTGCTAAAACCCGTGGTAATAAGCGTACCGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATGCGACAGTGCCACGTTCTCTGAAAAAAGGTCCttt                                                                                                   >  1:50242/1‑135 (MQ=255)
                                         tCGCCATGGGGCTTCCCAACCAAGGGTGCTAAAACCCGTGGTAATAAGCGTACCGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATGCGACAGTGCCACGTTCTCTGAAAAAAGGTCCTTTTATCGAt                                                                                            >  6:73567/1‑136 (MQ=255)
                                            ccATGGGGCTTCCCAACCAAGGGTGCTAAAACCCGTGGTAATAAGCGTACCGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATACGACAGTGCCACGTTCTCTGAAAAAAGGTCCTTTTATCGATCtt                                                                                         <  3:105693/136‑1 (MQ=255)
                                                    cTTCCCAACCAAGGGTGCTAAAACCCGTGGTAATAAGCGTACCGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATGCGACAGTGCCACGTTCTCTGAAAAAAGGTCCTTTTATCGATCTTCACCTGtt                                                                                 >  3:297256/1‑136 (MQ=255)
                                                         caaccaaGGGTGCTAAAACCCGTGGTAATAAGCGTACCGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATGCGACAGTGCCACGTTCTCTGAAAAAAGGTCCTTTTATCGATCTTCACCTGTTgaag                                                                             >  7:8735/1‑135 (MQ=255)
                                                         caaccaaGGGTGCTAAAACCCGTGGTAATAAGCGTACCGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATGCGACAGTGCCACGTTCTCTGAAAAAAGGTCCTTTTATCGATCTTCACCTGTTg                                                                                >  5:286234/1‑132 (MQ=255)
                                                          aaccaaGGGTGCTAAAACCCGTGGTAATAAGCGTACCGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATGCGACAGTGCCACGTTCTCTGAAAAAAGGTCCTTTTATCGATCTTCACCTGTTgaag                                                                             >  8:17008/1‑134 (MQ=255)
                                                          aaccaaGGGTGCTAAAACCCGTGGTAATAAGCGTACCGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATGCGACAGTGCCACGTTCTCTGAAAAAAGGTCCTTTTATCGATCTTCACCTGTTgaag                                                                             >  3:65561/1‑134 (MQ=255)
                                                          aaccaaGGGTGCTAAAACCCGTGGTAATAAGCGTACCGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATACGACAGTGCCACGTTCTCTGAAAAAAGGTCCTTTTATCGATCTTCACCTGTTgaag                                                                             >  7:92691/1‑134 (MQ=255)
                                                          aaccaaGGCTGCTAAAACCCGTGGTAATAAGCGTACCGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATGCGACAGTGCCACGTTCTCTGAAAAAAGGTCCTTTTATCGATCtt                                                                                         >  4:229711/1‑122 (MQ=255)
                                                          aaccaaGGCTGCTAAAACCCGTGGTAATAAGCGTACCGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATGCGACAGTGCCACGTTCTCTGAAAAAAGGTCCTTTTATCGATCtt                                                                                         <  3:229711/122‑1 (MQ=255)
                                                               aGGGTGCTAAAACCCGTGGTAATAAGCGTACCGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATGCGACAGTGCCACGTTCTCTGAAAAAAGGTCCTTTTATCGATCTTCACCTGTTGAAGAAGGTc                                                                       >  1:275093/1‑135 (MQ=255)
                                                                   tGCTAAAACCCGTGGTAATAAGCGTACCGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATGCGACAGTGCCACGTTCTCTGAAAAAAGGTCCTTTTATCGATCTTCACATGTTGAAGAAGGTCGAAGt                                                                  >  1:175553/1‑136 (MQ=255)
                                                                     cTAAAACCCGTGGTAATAAGCGTACCGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATGCGACAGTGCCACGTTCTCTGAAAAAAGGTCCTTTTATCGATCTTCACCTGTTGAAGAAGGTCGAAGTg                                                                 <  7:79320/135‑1 (MQ=255)
                                                                        aaaCCCGTGGTAATAAGCGTACCGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATGCGACAGTGCCACGTTCTCTGAAAAAAGGTCCTTTTATCGATCTTCACCTGTTGAAGAAGGTCGAAGTGGCgg                                                             >  8:73512/1‑136 (MQ=255)
                                                                                ggTAATAAGCGTACCGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATGCGACAGTGCCACGTTCTCTGAAAAAAGGTCCTTTTATCGATCTTCACCTGTTGAAGAAGGTCGAAgtggaggtggag                                                        >  3:240243/1‑133 (MQ=255)
                                                                                ggTAATAAGCGTACCGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATGCGACAGTGCCACGTTCTCTGAAAAAAGGTCCTTTTATCGATCTTCACCTGTTGAAGAAGGTCGAAGTGGCGGTGgag                                                        >  3:244297/1‑133 (MQ=255)
                                                                                  taataaGCGTACCGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATGCGACAGTGCCACGTTCTCTGAAAAAAGTTCATTTTATCGATCTTCACATTTTGAAGATGGTCTATTGGGCGGTGGAgaagaa                                                   >  3:121759/1‑136 (MQ=255)
                                                                                          gTACCGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATGCGACAGTGCCACGTTCTCTGAAAAAAGGTCCTTTTATCGATCTTCACCTGTTGAAGAAGGTCGAAGTGGCGGTGGAGAAGAACGATCGCa                                           <  8:218889/136‑1 (MQ=255)
                                                                                          gTACCGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATGCGACAGTGCCACGTTCTCTGAAAAAAGGTCCTTTTATCGATCTTCACCTGTTGAAGAAGGTCGAAGTGGCGGTGGAGAAGAACGATCGCa                                           <  7:92180/136‑1 (MQ=255)
                                                                                          gTACCGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATACGACAGTGCCACGTTCTCTGAAAAAAGGTCCTTTTATCGATCTTCACCTGTTGAAGAAGGTCGAAGTGGCGGTGGAGAAGAACGATCGCa                                           <  3:162642/136‑1 (MQ=255)
                                                                                              cGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATGCGACAGTGCCACGTTCTCTGAAAAAAGGTCCTTTTATCGATCTTCACCTGTTGAAGAAGGTCGAAGTGGCGGTGGAGAAGAACGATCGCAAGc                                        >  1:87128/1‑135 (MQ=255)
                                                                                                   acaTGATCGTCCGTCGTCGCAAGTAACTAGAGGGATGCGACAGTGCCACGTtc                                                                                                                     >  7:265568/1‑53 (MQ=255)
                                                                                                   acaTGATCGTCCGTCGTCGCAAGTAACTAGAGGGATGCGACAGTGCCACGTtc                                                                                                                     <  8:265568/53‑1 (MQ=255)
                                                                                                       gATCGTCCGTCGTCGCAAGTAACTAGAGGGATGCGACAGTGCCACGTTCTCTGAAAAAAGGTCCTTTTATCGATCTTCACCTGTTGAAGAAGGTCGAAGTGGCGGTGGAGAAGAACGATCGCAAGCCAGTTaaa                                <  2:136063/134‑1 (MQ=255)
                                                                                                       gATCGTCCGTCGTCGCAAGTAACTAGAGGGATGCGACAGTGCCACGTTCTCTGAAAAAAGGTCCTTTTATCGATCTTCACCTGTTGAAGAAGGTCGAAGTGGCGGTGGAGAAGAACGATCGCAAGCCAGTTaaa                                >  1:136063/1‑134 (MQ=255)
                                                                                                             ccgtcGTCGCAAGTAACTAGAGGGATGCGACAGTGCCACGTTCTCTGAAAAAAGGTCCTTTTATCGATCTTCACCTGTTGAAGAAGGTCGAAGTGGCGGTGGAGAAGAACGATCGCAAGCCAGTTAAAACCTGGTc                        <  6:162985/136‑1 (MQ=255)
                                                                                                              cgtcgtCGCAATTAACTAGAGGGATTCGACAGTGCCACGTTCTCTGAAAAAAGGTCCTTTTATCGATCTTCACCTGTTGAAGAAGGTCGAAGTGGCGGTGGAGAAGAACGATCGCAAGCCAGTTAAAACCTGGTc                        <  3:197033/135‑1 (MQ=255)
                                                                                                               gtcgtcGCAAGTAACTAGAGGGATGCGACAGTGCCACGTTCTCTGaaa                                                                                                              <  6:64868/48‑1 (MQ=255)
                                                                                                               gtcgtcGCAAGTAACTAGAGGGATGCGACAGTGCCACGTTCTCTGaaa                                                                                                              >  5:64868/1‑48 (MQ=255)
                                                                                                               gtcgtcGCAAGTAACTAGAGGGATGCGACAGTGCCACGTTCTCTGAAAAAAGGTCCTTTTATCGATCTTCACCTGTTGAAGAAGGTCGAAGTGGCGGTGGAGAAGAACGATCGCAAGCCAGTTAAAACCTGGTcgc                      >  6:210085/1‑136 (MQ=255)
                                                                                                                                    gATACGACAGTGCCACGTTCTCTGAAAAAAGGTCCTTTTATCGATCTTCACCTGTTGAAGAAGGTCGAAGTGGCGGTGGAGAAGAACGATCGCAAGCCAGTAAAAACCTGGTCGAGCCGTTCGATGATCCTGCca  >  7:193909/1‑135 (MQ=255)
                                                                                                                                       |                                                                                                                                   
GTGGTGGTGAAGGTCGTACCTCCGGTGGTCGTCATCCGGTATCGCCATGGGGCTTCCCAACCAAGGGTGCTAAAACCCGTGGTAATAAGCGTACCGACAACATGATCGTCCGTCGTCGCAAGTAACTAGAGGGATACGACAGTGCCACGTTCTCTGAAAAAAGGTCCTTTTATCGATCTTCACCTGTTGAAGAAGGTCGAAGTGGCGGTGGAGAAGAACGATCGCAAGCCAGTTAAAACCTGGTCGCGCCGTTCGATGATCCTGCCA  >  ACB143/546364‑546630

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: