Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,030,799 C→T 21.1% intergenic (+37/‑458) PP_4433 → / → dadA‑I amino acid MFS transporter/D‑amino acid dehydrogenase small subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,030,7990CT21.1% 66.5 / 15.7 38intergenic (+37/‑458)PP_4433/dadA‑Iamino acid MFS transporter/D‑amino acid dehydrogenase small subunit
Reads supporting (aligned to +/- strand):  ref base C (15/15);  new base T (3/5);  total (18/20)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-01

ACCTGCCTACTACATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTCGCTATCTCTTG  >  NC_002947/5030651‑5030944
                                                                                                                                                    |                                                                                                                                                 
aCCTGCCTACTACATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCa                                                                                                                                                  >  1:1156149/1‑150 (MQ=255)
    gCCTACTACATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTGATTAGCAc                                                                                                                                              >  1:308759/1‑150 (MQ=255)
     ccTACTACATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTGATTAGCAcc                                                                                                                                             <  2:308759/150‑1 (MQ=255)
       tactacATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTGATTAGCACCCa                                                                                                                                           <  2:147758/150‑1 (MQ=255)
       tactacATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTGATTAGCACCCa                                                                                                                                           <  2:689022/150‑1 (MQ=255)
        actacATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAg                                                                                                                                          <  1:118437/150‑1 (MQ=255)
        actacATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAg                                                                                                                                          >  1:376287/1‑150 (MQ=255)
          tacATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGc                                                                                                                                        >  2:947434/1‑150 (MQ=255)
          tacATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGc                                                                                                                                        <  2:376287/150‑1 (MQ=255)
                    ttGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAg                                                                                                                              >  1:836280/1‑150 (MQ=255)
                      gCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCt                                                                                                                            <  1:289771/150‑1 (MQ=255)
                       cATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTc                                                                                                                           <  2:836280/150‑1 (MQ=255)
                                               tgttCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTGATTAGCACCCAgg                                                                                                                                         <  1:459583/112‑1 (MQ=255)
                                               tgttCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTGATTAGCACCCAgg                                                                                                                                         >  2:459583/1‑112 (MQ=255)
                                                     ttGTGAACGTGACTTAACGTTGCTCGATTCGGTGGGATTGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGt                                                                                             <  2:561136/150‑1 (MQ=255)
                                                        tgATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGa                                                                                          >  1:596816/1‑150 (MQ=255)
                                                             gagaCTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATaa                                                                                     <  2:1156149/150‑1 (MQ=255)
                                                                 cTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAagcag                                                                                 >  1:941296/1‑150 (MQ=255)
                                                                           cTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTc                                                                                                                           >  1:758242/1‑98 (MQ=255)
                                                                           cTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTc                                                                                                                           <  2:758242/98‑1 (MQ=255)
                                                                            tCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCAc                                                                      >  2:363119/1‑150 (MQ=255)
                                                                            tCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCAc                                                                      >  1:970922/1‑150 (MQ=255)
                                                                              gATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGc                                                                                                          <  2:1115181/112‑1 (MQ=255)
                                                                              gATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGc                                                                                                          >  1:1115181/1‑112 (MQ=255)
                                                                                                   tGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTGATTAGCACCCAgg                                                                                                                                         >  2:292425/1‑60 (MQ=37)
                                                                                                   tGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTGATTAGCACCCAgg                                                                                                                                         <  1:292425/60‑1 (MQ=37)
                                                                                                       gCAAAATAACCTAATAGGCAAATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCACGCCAAAACGACTGGCTAGGGCGATATTAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACctc                                           >  1:745167/1‑150 (MQ=255)
                                                                                                          aaaTAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGAcc                                                                                                     <  1:816940/89‑1 (MQ=255)
                                                                                                          aaaTAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGAcc                                                                                                     >  2:816940/1‑89 (MQ=255)
                                                                                                                        gCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGa                          <  1:947434/150‑1 (MQ=255)
                                                                                                                             tGGTATCGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCAGTCTTGCGACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACa                     <  1:74693/150‑1 (MQ=255)
                                                                                                                             tGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACa                     <  2:941296/150‑1 (MQ=255)
                                                                                                                                     agCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTc             <  1:363119/150‑1 (MQ=255)
                                                                                                                                       ccTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAgt                   <  2:1021772/142‑1 (MQ=255)
                                                                                                                                       ccTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAgt                   >  1:1021772/1‑142 (MQ=255)
                                                                                                                                        cTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTc                                                       >  2:1045301/1‑105 (MQ=255)
                                                                                                                                        cTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTc                                                       <  1:1045301/105‑1 (MQ=255)
                                                                                                                                                tAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTCGCTATCTCTTg  >  2:911405/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ACCTGCCTACTACATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTCGCTATCTCTTG  >  NC_002947/5030651‑5030944

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