Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,479,750 C→A 11.8% Q353K (CAA→AAA)  ybdR → Zn‑dependent oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,479,7500CA11.8% 94.3 / 2.4 34Q353K (CAA→AAA) ybdRZn‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (14/16);  new base A (2/2);  total (16/18)
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.32e-01

AGCGTGCCCGGCGTCTATGCAGGGTTCATTCACGGCTTCATGTTCGGTGATGCCTTCGACAAAGGCCTGACCTTCAAGATGGGTCAGACCCATGTTCAAAAATACCTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTA  >  NC_002947/4479606‑4479895
                                                                                                                                                |                                                                                                                                                 
aGCGTGCCCGGCGTCTATGCAGGGTTCATTCACGGCTTCATGTTCGGTGATGCCTTCGACAAAGGCCTGACCTTCAAGATGGGTCAGACCCATGTTCAAAAATACCTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCg                                                                                                                                              <  1:255610/150‑1 (MQ=255)
 gCGTGCCCGGCGTCTATGCAGGGTTCATTCACGGCTTCATGTTCGGTGATGCCTTCGACAAAGGCCTGACCTTCAAGATGGTTCAGACCCATGTTCAAAAATACCTTCCCGAACTGCTGGAGCACATCGAAGCAGTGCGCCTGCAACCgg                                                                                                                                             >  2:248665/1‑150 (MQ=255)
                 tGCAGGGTTCATTCACGGCTTCATGTTCGGTGATGCCTTCGACAAAGGCCTGACCTTCAAGATGGGTCAGACCCATGTTCAAAAATACCTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCa                                                                                                                             >  2:113048/1‑150 (MQ=255)
                    aGGGTTCATTCACGGCTTCATGTTCGGTGATGCCTTCGACAAAGGCCTGACCTTCAAGATGGGTCAGACCCATGTTCAAAAATACCTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCg                                                                                                                          >  2:440947/1‑150 (MQ=255)
                               aCGGCTTCATGTTCGGTGATGCCTTCGACAAAGGCCTGACCTTCAAGATGGGTCAGACCCATGTTCAAAAATACCTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTag                                                                                                               >  2:170442/1‑150 (MQ=255)
                                     tCATGTTCGGTGATGCCTTCGACAAAGGCCTGACCTTCAAGATGGGTCAGACCCATGTTCAAAAATACCTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAgg                                                                                                         <  1:603144/150‑1 (MQ=255)
                                        tGTTCGGTGATGCCTTCGACAAAGGCCTGACCTTCAAGATGGGTCAGACCCATGTTCAAAAATACCTTCCCGAACTGCTTGAGCACATCGAAGCAGGGAGCCGTCCACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAggctg                                                                                                      >  2:146057/1‑148 (MQ=255)
                                                    ccTTCGACAAAGGCCTGACCTTCAAGATGGGTCAGACCCATGTTCAAAAATACCTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACa                                                                                          >  1:169030/1‑150 (MQ=255)
                                                      ttCGACAAAGGCCTGACCTTCAAGATGGGTCAGACCCATGTTCAAAAATACCTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGaa                                                                                                                                           >  1:791313/1‑99 (MQ=255)
                                                      ttCGACAAAGGCCTGACCTTCAAGATGGGTCAGACCCATGTTCAAAAATACCTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGaa                                                                                                                                           <  2:791313/99‑1 (MQ=255)
                                                      ttCGACAAAGGCCTGACCTTCAAGATGGGTCAGACCCATGTTCAAAAATACCTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAggcgg                                                                                                      >  2:527587/1‑136 (MQ=255)
                                                      ttCGACAAAGGCCTGACCTTCAAGATGGGTCAGACCCATGTTCAAAAATACCTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAggcgg                                                                                                      <  1:527587/136‑1 (MQ=255)
                                                        cGACAAAGGCCTGACCTTCAAGATGGGTCAGACCCATGTTCAAAAATACCTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAAt                                                                                      <  1:514366/150‑1 (MQ=255)
                                                         gACAAAGGCCTGACCTTCAAGATGGGTCAGACCCATGTTCAAAAATACCTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGAAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATg                                                                                     >  1:82090/1‑150 (MQ=255)
                                                         gACAAAGGCCTGACCTTCAAGATGGGTCAGACCCATGTTCAAAAATACCTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGAAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATg                                                                                     >  1:82141/1‑150 (MQ=255)
                                                            aaaGGCCTGACCTTCAAGATGGGTCAGACCCATGTTCAAAAATACCTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCCGACCTGATCGTTATTCACCTCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTc                                                                                  <  2:689839/150‑1 (MQ=255)
                                                              aGGCCTGACCTTCAAGATGGGTCAGACCCATGTTCAAAAATACCTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGa                                                                                <  1:170442/150‑1 (MQ=255)
                                                               ggCCTGACCTTCAAGATGGGTCAGACCCATGTTCAAAAATACCTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGAAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGAt                                                                               <  2:82141/150‑1 (MQ=255)
                                                               ggCCTGACCTTCAAGATGGGTCAGACCCATGTTCAAAAATACCTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGAAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGAt                                                                               <  2:82090/150‑1 (MQ=255)
                                                                       cTTCAAGATGGGTCAGACCCATGTTCAAAAATACCTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCa                                                                       >  2:757028/1‑150 (MQ=255)
                                                                              aTGGGTCAGACCCATGTTCAAAAATACCTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAAc                                                                <  2:601096/150‑1 (MQ=255)
                                                                               tGGGTCAGACCCATGTTCAAAAATACCTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACt                                                               <  1:440947/150‑1 (MQ=255)
                                                                                           aTGTTCAAAAATACCTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCa                                                   >  1:615136/1‑150 (MQ=255)
                                                                                              ttCAAAAATACCTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGtt                                                                                   >  1:475403/1‑115 (MQ=255)
                                                                                              ttCAAAAATACCTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGtt                                                                                   <  2:475403/115‑1 (MQ=255)
                                                                                               tCAAAAATACCTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCt                                               <  1:248665/150‑1 (MQ=255)
                                                                                                 aaaaaTACCTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCg                                             <  1:113048/150‑1 (MQ=255)
                                                                                                      tACCTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCa                                        <  1:757028/150‑1 (MQ=255)
                                                                                                               gAACTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCa                               >  1:770073/1‑150 (MQ=255)
                                                                                                                        gAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGt                                            <  2:162187/128‑1 (MQ=255)
                                                                                                                        gAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGt                                            >  1:162187/1‑128 (MQ=255)
                                                                                                                                           gcCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAg                                                       >  2:384187/1‑98 (MQ=255)
                                                                                                                                           gcCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAg                                                       <  1:384187/98‑1 (MQ=255)
                                                                                                                                            ccTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTa  <  2:169030/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
AGCGTGCCCGGCGTCTATGCAGGGTTCATTCACGGCTTCATGTTCGGTGATGCCTTCGACAAAGGCCTGACCTTCAAGATGGGTCAGACCCATGTTCAAAAATACCTTCCCGAACTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTA  >  NC_002947/4479606‑4479895

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