Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 550,552 C→G 12.0% intergenic (+30/‑112) tufB → / → rpsJ elongation factor Tu‑B/30S ribosomal protein S10

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947550,5520CG12.0% 108.4 / 8.7 42intergenic (+30/‑112)tufB/rpsJelongation factor Tu‑B/30S ribosomal protein S10
Reads supporting (aligned to +/- strand):  ref base C (19/18);  new base G (3/2);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.37e-01

AACATTCAGATGACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAGGCTTT  >  NC_002947/550403‑550701
                                                                                                                                                     |                                                                                                                                                     
aacaTTCAGATGACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTc                                                                                                                                                       >  2:20944/1‑150 (MQ=255)
aacaTTCAGATGACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTc                                                                                                                                                       >  2:143758/1‑150 (MQ=255)
  caTTCAGATGACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAg                                                                                                                                                     >  2:307359/1‑150 (MQ=255)
      cAGATGACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCggg                                                                                                                                                 <  1:307359/150‑1 (MQ=255)
             cTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTGAGCGGGGGttttt                                                                                                                                          <  2:220029/150‑1 (MQ=255)
             cTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGttttt                                                                                                                                          <  1:162496/150‑1 (MQ=255)
                      tGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTggg                                                                                                                                 <  1:20944/150‑1 (MQ=255)
                        aTCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGtt                                                                                                                               >  1:137029/1‑150 (MQ=255)
                                    gCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTgg                                                                                                                   <  1:143758/150‑1 (MQ=255)
                                       aTGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTggggg                                                                                                                <  2:137029/150‑1 (MQ=255)
                                                  tCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTGAGCGGGGGtttt                                                                                                                                           >  1:227272/1‑112 (MQ=255)
                                                  tCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTGAGCGGGGGtttt                                                                                                                                           <  2:227272/112‑1 (MQ=255)
                                                            gCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGcctcc                                                                                           >  1:162032/1‑150 (MQ=255)
                                                              tATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTCATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTggg                                                                                                                                 <  2:460105/110‑1 (MQ=255)
                                                              tATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTggg                                                                                                                                 >  1:460105/1‑110 (MQ=255)
                                                                 ccGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAACCCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTa                                                                                      <  1:306176/150‑1 (MQ=255)
                                                                  cGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCgg                                                                                                                                                  >  1:421021/1‑89 (MQ=255)
                                                                  cGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCgg                                                                                                                                                  <  2:421021/89‑1 (MQ=255)
                                                                      aaGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCggg                                                                                 >  1:611061/1‑150 (MQ=255)
                                                                         gCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGtttttt                                                                                                                                         <  2:640460/91‑1 (MQ=255)
                                                                         gCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGtttttt                                                                                                                                         <  2:640438/91‑1 (MQ=255)
                                                                         gCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGtttttt                                                                                                                                         >  1:640438/1‑91 (MQ=255)
                                                                         gCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGtttttt                                                                                                                                         >  1:640460/1‑91 (MQ=255)
                                                                                            cGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGt                                                                       <  2:40547/138‑1 (MQ=255)
                                                                                            cGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGt                                                                       >  1:40547/1‑138 (MQ=255)
                                                                                                       ccAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCt                                                >  2:332346/1‑150 (MQ=255)
                                                                                                           aaTCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGaaa                                            <  2:685462/150‑1 (MQ=255)
                                                                                                           aaTCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGaaa                                            >  1:622850/1‑150 (MQ=255)
                                                                                                           aaTCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGaaa                                            <  2:162032/150‑1 (MQ=255)
                                                                                                              ggtcGAGTGATTACTCGACCGATTGAAAAAACCCCCGCGCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAAtac                                         >  2:92031/3‑148 (MQ=255)
                                                                                                               gtcGAGTGATTACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCa                                        >  1:628007/2‑150 (MQ=255)
                                                                                                                  gAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTc                                                                                                            >  2:224442/1‑79 (MQ=255)
                                                                                                                  gAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTc                                                                                                            <  1:224442/79‑1 (MQ=255)
                                                                                                                  gAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATg                                     <  1:332346/150‑1 (MQ=255)
                                                                                                                   aGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTGTTTATTGGGTTGACACTAATTGTGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCTTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGc                                    >  2:651037/1‑150 (MQ=255)
                                                                                                                       aTCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCaaaa                                >  1:393363/1‑150 (MQ=255)
                                                                                                                                   aTTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCg                    >  1:777211/1‑150 (MQ=255)
                                                                                                                                   aTTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCg                    >  1:74300/1‑150 (MQ=255)
                                                                                                                                      gAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTAt                 <  2:74300/150‑1 (MQ=255)
                                                                                                                                              ccccGCTGAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGa         >  1:785977/1‑150 (MQ=255)
                                                                                                                                               cccGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCa               >  2:323053/1‑143 (MQ=255)
                                                                                                                                               cccGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCAGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCa               <  1:323053/143‑1 (MQ=255)
                                                                                                                                                ccGCTGAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAg       >  1:394728/1‑150 (MQ=255)
                                                                                                                                                ccGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAg       <  2:628007/150‑1 (MQ=255)
                                                                                                                                                     cAGCGGGGGTTTTTTTTATTTGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAGGCttt  <  1:457274/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AACATTCAGATGACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAGGCTTT  >  NC_002947/550403‑550701

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