Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 550,760 G→A 50.0% G33S (GGT→AGT)  rpsJ → 30S ribosomal protein S10

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947550,7600GA50.0% ‑3.6 / 28.0 24G33S (GGT→AGT) rpsJ30S ribosomal protein S10
Reads supporting (aligned to +/- strand):  ref base G (6/6);  new base A (7/5);  total (13/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAGGCTTTCGACCATCGCCTGATCGACCAATCCACCCAGGAAATCGTGGAAACCGCGAAACGTACTGGTGCACAAGTGCGTGGTCCAATTCCACTGCCTACCCGCAAAGAGCGTTTCACCGTTCTGGTCTCCCCGCACGTCAACAAAGACGCGCGTGACCAGTACGAGATTCGCACTCATAAGCGTGTTCTGGACATCGTCCAGCCAACGGA  >  NC_002947/550611‑550905
                                                                                                                                                     |                                                                                                                                                 
cctTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAGGCTTTCGACCATCGCCTGATCGACCAATCCACCCAGGAAATCGTGGAAACCGCGAAACGTACTg                                                                                                                                                   <  2:414406/150‑1 (MQ=255)
          ggCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAGGCTTTCGACCATCGCCTGATCGACCAATCCACCCAGGAAATCGTGGAAACCGCGAAACGTACTGGTGCACAAgc                                                                                                                                         <  2:387785/150‑2 (MQ=255)
            cGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAGGCTTTCGACCATCGCCTGATCGACCAATCCACCCAGGAAATCGTGGAAACCGCGAAACGTACTGGTGCACAAGTGc                                                                                                                                       <  2:438967/150‑1 (MQ=255)
                gTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAGGCTTTCGACCATCGCCTGATCGACCAATCCACCCAGGAAATCGTGGAAACCGCGAAACGTACTAGTGCACAAGTGCGTgg                                                                                                                                   >  2:22903/1‑150 (MQ=255)
                            ggAACAGCCTGGAGTCTGAAATCCAAAATCGAATCAGCAAATCCGTATCAGGTTGAAGGCTTTCGACCATCGCCTGATCGACCAATCCACCCAGGAAATCGTGGAAACCGCGAAACGTACTGGTGCACAAGTGCGTGGTCCAATTCCACt                                                                                                                       >  1:385316/1‑150 (MQ=255)
                              aaCAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAGGCTTTCGACCATCGCCTGATCGACCAATCCACCCAGGAAATCGTGGAAACCGCGAAACGTACTGGTGCACAAGTGCGTGGTCCAATTCCACTGc                                                                                                                     >  2:151109/1‑150 (MQ=255)
                                      ggAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAGGCTTTCGACCATCGCCTGATCGACCAATCCACCCAGGAAATCGTGGAAACCGCGAAACGTACTGGTGCACAAGTGCGTGGTCCAATTCCACTGCCTACCCGc                                                                                                             <  2:18267/150‑1 (MQ=255)
                                         gTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAGGCTTTCGACCATCGCCTGATCGACCAATCCACCCAGGAAATCGTGGAAACCGCGAAACGTACTAGTGCACAAGTGCGTGGTCCAATTCCACTGCCTACCCGCaaa                                                                                                          >  2:25618/1‑150 (MQ=255)
                                               aaTCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAGGCTTTCGACCATCGCCTGATCGACCAATCCACCCAGGAAATCGTGGAAACCGCGAAACGTACTAGTGCACAAGTGCGTGGTCCAATTCCACTGCCTACCCGCAAAGAGCGt                                                                                                    <  1:203760/150‑1 (MQ=255)
                                                 tCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAGGCTTTCGACCATCGCCTGATCGACCAATCCACCCAGGAAATCGTGGAAACCGCGAAACGTACTGGTGCACAAGTGCGTGGTCCAATTCCACTGCCTACCCGCAAAGAGCGttt                                                                                                  >  1:491709/1‑150 (MQ=255)
                                                           aaTCAGCAAATCCGTATCAGGTTGAAGGCTTTCGACCATCGCCTGATCGACCAATCCACCCAGGAAATCGTGGAAACCGCGAAACGTACTGGTGCACAAGTGCGTGGTCCAATTCCACTGCCTACCCGCAAAGAGCGTTTCACCGTTCTg                                                                                        <  2:491709/150‑1 (MQ=255)
                                                             tCAGCAAATCCGTATCAGGTTGAAGGCTTTCGACCATCGCCTGATCGACGAATCCACCCAGGAAATCGTGGAAACCGCGAAACGTACTGGTGCACAAGTGCGTGGTCCAATTCCACTGCCTACCCGCAAAGAGCGTTTCACCGTTCTGGt                                                                                      <  2:385316/150‑1 (MQ=255)
                                                                aggaagccGTATCAGGTTGAAGGCTTTCGACCATCGCCTGATCGACCAATCCACCCAGGAAATCGTGGAAACCGCGAAACGTACTAGTGCACAAGTGCGTGGACAAATTCCACTGCCTAACAGCAAAGGGCGTTTCACCGTTCTGGgctc                                                                                   >  1:443149/7‑150 (MQ=255)
                                                                                    aGGCTTTCGACCATCGCCTGATCGACCAATCCACCCAGGAAATCGTGGAAACCGCGAAACGTACTAGTGCACAAGTGCGTGGTCCAATTCCACTGCCTACCCGCAAAGAGCGTTTCACCGTTCTGGTCTCCCCGCACGTCAACAAAGAcg                                                               >  1:535519/1‑150 (MQ=255)
                                                                                        tttCGACCATCGCCTGATCGACCAATCCACCCAGGAAATCGTGGAAACCGCGAAACGTACTGGTGCACAAGTGCGTGGTCCAATTCCACTGCCTACCCGCAAAGAGCGTTTCACCGTTCTGGTCTCCCCGCACGTCAACAAAGAcgcgcg                                                           >  2:540220/1‑150 (MQ=255)
                                                                                            gACCATCGCCTGATCGACCAATCCACCCAGGAAATCGTGGAAACCGCGAAACGTACTAGTGCACAAGTGCGTGGTCCAATTCCACTGCCTACCCGCAAACAGCGGTTCACCGTTCTGGTCTCCCCGCAc                                                                            >  2:438525/1‑129 (MQ=255)
                                                                                            gACCATCGCCTGATCGACCAATCCACCCAGGAAATCGTGGAAACCGCGAAACGTACTAGTGCACAAGTGCGTGGTCCAATTCCACTGCCTACCCGCAAACACCGGTTCACCGTTCTGGTCTCCCCGCAc                                                                            <  1:438525/129‑1 (MQ=255)
                                                                                                            aCCAATCCACCCAGGAAATCGTGGAAACCGCGAAACGTACTGGTGCACAAGTGCGTGGTCCAATTCCACTGCCTACCCGCAAAGAGCGTTTCACCGTTCTGGTCTCCCCGCACGTCAACAAAGACGCGCGTGACCAGTACGAGATTCGCa                                       >  1:264493/1‑150 (MQ=255)
                                                                                                                   cACCCAGGAAATCGTGGAAACCGCGAAACGTACTAGTGCACAAGTGCGTGGTCCAATTCCACTGCCTACCCGCAAAGAGCGTTTCACCGTTCTGGTCTCCCCGCACGTCAACAAAGACGCGCGTGACCAGTACGAGATTCGCACTCATaa                                <  2:157475/150‑1 (MQ=255)
                                                                                                                        aGGAAATCGTGGAAACCGCGAAACGTACTGGTGCACAAGTGCGTGGTCCAATTCCACTGCCTACCCGCAAAGAGCGTTTCACCGTTCTGGTCTCCCCGCACGTCAACAAAGACGCGCGTGACCAGTACGAGATTCGCACTCATAAGCgtg                           <  2:426648/150‑1 (MQ=255)
                                                                                                                        aGGAAATCGTGGAAACCGCGAAACGTACTAGTGCACAAGTGCGTGGTCCAATTCCACTGCCTACCCGCAAAGAGCGTTTCACCGTTCTGGTCTCCCCGCACGTCAACAAAGACGCGCGTGACCAGTACGAGATTCGCACTCATAAGCgtg                           >  1:459136/1‑150 (MQ=255)
                                                                                                                                   gAAACCGCGAAACGTACTAGTGCACAAGTGCGTGGTCCAATTCCACTGCCTACCCGCAAAGAGCGTTTCACCGTTCTGGTCTCCCCGCACGTCAACAAAGACGCGCGTGACCAGTACGAGATTCGCACTCATAAGCGTGTTCTGGACATc                <  1:22903/150‑1 (MQ=255)
                                                                                                                                        cgcgAAACGTACTAGTGCACAAGTGCGTGGTCCAATTCCACTGCCTACCCGCAAAGAGCGTTTCACCGTTCTGGTCTCCCCGCACGTCAACAAAGACGCGCGTGACCAGTACGAGATTCGCACTCATAAGCGTGTTCTGGACATCGTccc           >  2:113734/1‑149 (MQ=255)
                                                                                                                                              aCGTACTAGTGCACAAGTGCGTGGTCCAATTCCACTGCCTACCCGCAAAGAGCGTTTCACCGTTCTGGTCTCCCCGCACGTCAACAAAGACGCGCGTGACCAGTACGAGATTCGCACTCATAAGCGTGTTCTGGACATCGTCCAGCCAAc     <  1:25618/150‑1 (MQ=255)
                                                                                                                                                 tACTGGTGCACAAGTGCGTGGTCCAATTCCACTGCCTACCCGCAAAGAGCGTTTCACCGTTCTGGTCTCCCCGCACGTCAACAAAGACGCGCGTGACCAGTACGAGATTCGCACTCATAAGCGTGTTCTGGACATCGTCCAGCCAACGGa  >  2:67158/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAGGCTTTCGACCATCGCCTGATCGACCAATCCACCCAGGAAATCGTGGAAACCGCGAAACGTACTGGTGCACAAGTGCGTGGTCCAATTCCACTGCCTACCCGCAAAGAGCGTTTCACCGTTCTGGTCTCCCCGCACGTCAACAAAGACGCGCGTGACCAGTACGAGATTCGCACTCATAAGCGTGTTCTGGACATCGTCCAGCCAACGGA  >  NC_002947/550611‑550905

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: