Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,395,628 A→T 10.4% intergenic (‑193/+50) dacB ← / ← PP_2099 D‑alanyl‑D‑alanine carboxypeptidase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,395,6280AT10.4% 137.5 / 7.5 50intergenic (‑193/+50)dacB/PP_2099D‑alanyl‑D‑alanine carboxypeptidase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (23/21);  new base T (2/3);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.03e-01

GAGTGAATACCCTGCGGCGAGGGAAAAGGGCTTTGGGGCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTG‑ACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATCC‑AGCAGCGGGCTGACTTCGATCTTGGTCAGTTCGCTGCGGCCTTT  >  NC_002947/2395481‑2395764
                                                                                                                                                   |                                                                                                                                          
gAGTGAATACCCTGCGGCGAGGGAAAAGGGCTTTGGGGCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATg                                                                                                                                          >  1:251227/1‑150 (MQ=255)
  gTGAATACCCTGCGGCGAGGGAAAAGGGCTTTGGGGCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATggg                                                                                                                                        <  2:54223/150‑1 (MQ=255)
     aaTACCCTGCGGCGAGGGAAAAGGGCTTTGGGGCTGTAAGAGATGATGGCCCCACGACGTAAACGGTACGCATTATGCCCCAAGCGCAGCCGCCATGCGTCGCGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCCCCCCAGTCGGGTg                                                                                                                                     <  1:843977/150‑1 (MQ=255)
        aCCCTGCGGCGAGGGAAAAGGGCTTTGGGGCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCCATTCGGGTGGcc                                                                                                                                  <  1:598264/150‑1 (MQ=255)
          ccTGCGGCGAGGGAAAAGGGCTTTGGGGCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCCTCTCCAATAAAAAGGCCACCCCATTCGGGTGGCCtt                                                                                                                                <  1:202991/150‑1 (MQ=255)
          ccTGCGGCGAGGGAAAAGGGCTTTGGGGCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCTCCAATAAAAAGGCCACCCCATTCGGGTGGCCtt                                                                                                                                >  2:202991/1‑150 (MQ=255)
            tGCGGCGAGGGAAAAGGGCTTTGGGGCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCtttt                                                                                                                              <  2:656424/150‑1 (MQ=255)
                  gAGGGAAAAGGGCTTTGGGGCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCaa                                                                                                                        <  1:509933/150‑1 (MQ=255)
                  gAGGGAAAAGGGCTTTGGGGCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCCATTCGGGTGGCCtttttattgg                                                                                                                        <  1:645110/150‑7 (MQ=255)
                    gggAAAAGGGCTTTGGGGCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATc                                                                                                                      >  2:744438/1‑150 (MQ=255)
                     ggAAAAGGGCTTTGGGGCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCa                                                                                                                     >  2:770874/1‑150 (MQ=255)
                          aGGGCTTTGGGGCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCtt                                                                                                                >  2:344959/1‑150 (MQ=255)
                           gggCTTTGGGGCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCttt                                                                                                               >  1:54489/1‑150 (MQ=255)
                             gCTTTGGGGCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCa                                                                                                             >  2:31565/1‑150 (MQ=255)
                               tttGGGGCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGCTGCGTCATTCCCCACGGAAAACGATCGCCTATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGc                                                                                                           <  1:744438/150‑1 (MQ=255)
                                 tGGGGCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGccc                                                                                                         <  2:54489/150‑1 (MQ=255)
                                  ggggCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGAGAGCGGTTGGCAGCAATCAACCTTTCAGCCCg                                                                                                        >  1:641444/1‑150 (MQ=255)
                                         tAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTG‑a                                                                                                >  2:742880/1‑150 (MQ=255)
                                           agagATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCt                                                                                                   >  1:704486/1‑146 (MQ=255)
                                           agagATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCt                                                                                                   <  2:704486/146‑1 (MQ=255)
                                               atgatgTCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTG‑ACAAGTa                                                                                          >  1:676159/1‑150 (MQ=255)
                                                     gCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTG‑ACAAGTATCACTc                                                                                    >  1:836825/1‑150 (MQ=255)
                                                        ccACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAACCAAACTTTAAGACCGAAGCGT‑ACGAATAGCACGCAGc                                                                                 >  2:398837/1‑150 (MQ=255)
                                                             acgTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTG‑ACAAGTATCACTCAGCCttgt                                                                            <  1:742880/150‑1 (MQ=255)
                                                             acgTAAACAGTACGCATTATGCCCCAAGAGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTGAGCCACAATCTG‑ACAAGTATCACTCAGCCttgt                                                                            <  1:627371/150‑1 (MQ=255)
                                                                tAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTG‑ACAAGTATCACTCAGCCTTGTTGa                                                                         <  2:251227/150‑1 (MQ=255)
                                                                  aaCAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTG‑ACAAGTATCACTCAGCCTTGTTGATa                                                                       >  1:783883/1‑150 (MQ=255)
                                                                    cAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTG‑ACAAGTATCACTCAGCCTTGTTGATAgg                                                                     <  1:42226/150‑1 (MQ=255)
                                                                     aGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTG‑ACAAGTATCACTCAGCCTTGTTGATAGGc                                                                    <  2:133340/150‑1 (MQ=255)
                                                                       tACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTG‑ACAAGTATCACTCAGCCTTGTTGATAGGCtt                                                                  >  2:699772/1‑150 (MQ=255)
                                                                           cATTATGGCCCAAGCCCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGACAATAAAAAGGCCACCCGAATGGGGTGGCCGTTTCAGCAATCAACCTTTCAGCCCGAAGCTG‑ACAAGTATCACTCAGCCTTGTTGATAGGCTTTTcc                                                              <  2:87602/150‑1 (MQ=255)
                                                                             ttATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTG‑ACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCgg                                                            <  2:783883/150‑1 (MQ=255)
                                                                                    ccAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTG‑ACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAg                                                     <  1:508023/150‑1 (MQ=255)
                                                                                       aGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCGGCCCGAAGCTG‑ACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCa                                                  <  1:770874/150‑1 (MQ=255)
                                                                                         cgcAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTG‑ACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATc                                                <  1:398837/150‑1 (MQ=255)
                                                                                                    tGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTG‑ACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTATGTGG‑ATCCGAGCAGCGGGc                                    <  2:676159/150‑1 (MQ=255)
                                                                                                        tCACGTTCCACGGAAAAGCATCGCCAATAAAAAGGGCACCCGAATGAGGTGACCTTTTCAGCAAACAACCTTTCAGCCCGAAGCTG‑ACAAGTATCACTCAGCCTTGTTGAGAGGCTTTTCCGGCTACCAGACATCC‑AGCAGTGGGCTGAc                                >  2:361380/1‑150 (MQ=255)
                                                                                                            gTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTG‑ACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATCC‑AGCAGCGGGCTGACTTCg                            >  2:799808/1‑150 (MQ=255)
                                                                                                            gTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCACGAAGCTG‑ACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATCC‑AGCAGCGGGCTGACTTCg                            >  2:332987/1‑150 (MQ=255)
                                                                                                             ttCGACGGAAAAGCATCGCCAATAAAAAGGCAACCCGAACGGGGACGACTTTTAAGCAACCCACCATTGAGCCC‑AAACTGCACAGGTGTCACTCAGCCTTGTGCATAGGCTTTTCCGGGTACCAGGCATCC‑AGCAGCGGGCTGCCTTCgg                           >  2:313251/1‑149 (MQ=255)
                                                                                                               cGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTG‑TCAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATCC‑AGCAGCGGGCTGACTTCGATc                         <  1:699772/150‑1 (MQ=255)
                                                                                                                gACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTG‑ACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATc                                                >  2:194867/1‑127 (MQ=255)
                                                                                                                gACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTG‑ACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATc                                                <  1:194867/127‑1 (MQ=255)
                                                                                                                  cGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTG‑ACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATCC‑AGCAGCGGGCTGACTTCGATCTTg                      <  2:836825/150‑1 (MQ=255)
                                                                                                                   ggAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTG‑ACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATCC‑AGCAGCGGGCTGACTTCGATCTTgg                     >  1:377744/1‑150 (MQ=255)
                                                                                                                               ccAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTG‑ACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATCC‑AGCAGCGGGCTGACTTCGATCTTGGTCAGTTCGCTGc         >  2:572714/1‑150 (MQ=255)
                                                                                                                                  aTAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTG‑ACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATCC‑AGCAGCGGGCTGACTTCGATCTTGGTCAGTTCGCTGCGGc      >  1:144527/1‑150 (MQ=255)
                                                                                                                                      aaaGGCCACCCCATTCGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTG‑ACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACAAGGCATCC‑AGCAGGGGGTAGACTGCGGTCTTGGGCATTGCTCTGCGGCCttt  >  2:8035/1‑150 (MQ=255)
                                                                                                                                             aCCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTG‑ACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATCC‑AGCAGCGGGCTGACTTCGATCTTg                      >  2:268073/1‑123 (MQ=255)
                                                                                                                                             aCCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTG‑ACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATCC‑AGCAGCGGGCTGACTTCGATCTTg                      <  1:268073/123‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
GAGTGAATACCCTGCGGCGAGGGAAAAGGGCTTTGGGGCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTG‑ACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATCC‑AGCAGCGGGCTGACTTCGATCTTGGTCAGTTCGCTGCGGCCTTT  >  NC_002947/2395481‑2395764

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: ‑