Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,024,308 T→A 12.5% intergenic (+123/+31) dapA‑II → / ← PP_2640 4‑hydroxy‑tetrahydrodipicolinate synthase/GNAT family acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,024,3080TA12.5% 100.3 / 6.4 40intergenic (+123/+31)dapA‑II/PP_26404‑hydroxy‑tetrahydrodipicolinate synthase/GNAT family acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (15/20);  new base A (3/2);  total (18/22)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

CTGTGTTCACTCGGCGCCTGAAGCGGCAACGGTTGATGGCGGGTTGACAAGCCGCCACACCCTTGGTTGGTTCGGGGTCGGCGGATGCCTCGCGAGCGGATTTCTGCCAGCAACGCGCTGAGACAGAGTATTGCCAGGGCCGCTTTGCGGCCCTGTCGCGACCCCAGGTTGCTCCTACGCTGGCTGGGTCGGCAGGCGCTTTTCCATCGGCTGATACACCAGCCCGGCAAACTCACCCACGTCGCCTGCCAGGGTGAAGCCGTAGCGTTGGTAGGCGGGTACAGAAG  >  NC_002947/3024165‑3024451
                                                                                                                                               |                                                                                                                                                
cTGTGTTCACTCGGCGCCTGAAGCGGCAACGGTTGATGGCGGGTTGACAAGCCGCCACACCCTTGGTTGGTTCGGGGTCGGCGGATGCCTCGCGAGTGGATTTCTGCCAGCAACGCGCTGAGACAGAGTATTGCCAGGGCCGCTTTGCg                                                                                                                                             <  5:110338/149‑1 (MQ=255)
   tgtTCACTCGGCGCCTGAAGCGGCAACGGTTGATGGCGGGTTGACAAGCCGCCACACCCTTGGTTGGTTCGGGGTCGGCGGATGCCTCGCGAGCGGATTTCTGCCAGCAACGCGCTGAGACAGAGTATTGCCAGGGCCGCTTTGCGGcc                                                                                                                                          <  4:119572/149‑1 (MQ=255)
     ttCACTCGGCGCCTGAAGCGGCAACGGTTGATGGCGGGTTGACAAGCCGCCACACCCTTGGTTGGTTCGGGGTCGGCGGATGCCTCGCGAGCGGATTTCTGCCAGCAACGCGCTGAGACAGAGTATTGCCAGGGCCGCTTTGCGGCCCt                                                                                                                                        <  5:8944/149‑1 (MQ=255)
        aCTCGGCGCCTGAAGCGGCAACGGTTGATGGCGGGTTGACAAGCCGCCACACCCTTGGTTGGTTCGGGGTCGGCGGATGCCTCGCGAGCGGATTTCTGCCAGCAACGCGCTGAGACAGAGTATTGACAGGGCCGCAAAGCGGCCCTgg                                                                                                                                      <  1:9549/148‑2 (MQ=255)
         cTCGGCGCCTGAAGCGGCAACGGTTGATGGCGGGTTGACAAGCCGCCACACCCTTGGTTGGTTCGGGGTCGGCGGATGCCTCGCGAGCGGATTTCTGCCAGCAACGCGCTGAGACAGAGTATTGCCAGGGCCGCTTTGCGGCCCTGTcg                                                                                                                                    <  3:252081/149‑1 (MQ=255)
            ggCGCGTGAAGCGGCATGGGTTGATGGCTGGTTGACAAGGCGCCACACCCTTTGTTGGTTCGTGGTCGGCGGATGCCTCGCGAGCGGATTTCTGCCAGCAACTCGCTGAGACAGAGTATTGCCAGGGCCGCTTTGCGGCCCTGTcgc                                                                                                                                   <  3:231577/147‑1 (MQ=255)
               gcCTGAAGCGGCAACGGTTGATGGCGGGTTGACAAGCCGCCACACCCTTGGTTGGTTCGGGGTCGGCGGATGCCTCGCGAGCGGATTTCTGCCAGCAACGCGCTGAGACAGAGTATTGCCAGGGCCGCAAAGCGGCCCt                                                                                                                                        <  4:175064/139‑1 (MQ=255)
               gcCTGAAGCGGCAACGGTTGATGGCGGGTTGACAAGCCGCCACACCCTTGGTTGGTTCGGGGTCGGCGGATGCCTCGCGAGCGGATTTCTGCCAGCAACGCGCTGAGACAGAGTATTGCCAGGGCCGCAAAGCGGCCCt                                                                                                                                        >  3:175064/1‑139 (MQ=255)
                                    tGGCGGGTTGACAAGCCGCCACACCCTTGGTTGGTTCGGGGTCGGCGGATGCCTCGCGAGCGGATTTCTGCCAGCAACGCGCTGAGACAGAGTATTGCCAGGGCCGCTTTGCGGCCCTGTCGCGACCCCAGGTTGATCCTACgctggct                                                                                                         >  5:247898/1‑149 (MQ=255)
                                      gCGGGTTGACAAGCCGCCACACCCTTGGTTGGTTCGGGGTCGGCGGATGCCTCGCGAGCGGATTTCTGCCAGCAACGCGCTGAGACAGAGTATTGCCAGGGCCGCTTTGCGGCCCTGTCGCGACCCCAGGTTGCTCCTACgctggctg                                                                                                        <  4:153851/148‑1 (MQ=255)
                                      gCGGGTTGACAAGCCGCCACACCCTTGGTTGGTTCGGGGTCGGCGGATGCCTCGCGAGCGGATTTCTGCCAGCAACGCGCTGAGACAGAGTATTGCCAGGGCCGCTTTGCGGCCCTGTCGCGACCCCAGGTTGCTCCTACGctggctgg                                                                                                       <  6:175833/149‑1 (MQ=255)
                                        gggTTGACAAGCCGCCACACCCTTGGTTGGTTCGGGGTCGGCGGATGCCTCGCGAGCGGATTTCTGCCAGCAACGCGCTGAGACAGAGTATTGCCAGGGCCGCTTTGCGGCCCTGTcg                                                                                                                                    >  2:92684/1‑118 (MQ=255)
                                        gggTTGACAAGCCGCCACACCCTTGGTTGGTTCGGGGTCGGCGGATGCCTCGCGAGCGGATTTCTGCCAGCAACGCGCTGAGACAGAGTATTGCCAGGGCCGCTTTGCGGCCCTGTCGCGACCCCAGGTTGCTCCTACGCTGGCTGGGt                                                                                                     <  7:154876/149‑1 (MQ=255)
                                         ggTTGACAAGCCGCCACACCCTTGGTTGGTTCGGGGTCGGCGGATGCCTCGCGAGCGGATTTCTGCCAGCAACGCGCTGAGACAGAGTATTGCCAGGGCCGCTTTGCGGCCCTGTCGCGACCCCAGGTTGCTCCTACGCTGGATGGGTc                                                                                                    >  2:198480/1‑149 (MQ=255)
                                           ttGACAAGCCGCCACACCCTTGGTTGGTTCGGGGTCGGCGGATGCCTCGCGAGCGGATTTCTGCCAGCAACGCGCTGAGACAGAGTATTGCCAGGGCCGCTTTGCGGCCCTGTCGCGACCCCAGGTTGCTCCTACGCTGGCTGGGTCgg                                                                                                  >  3:239361/1‑149 (MQ=255)
                                                aaGCCGCCACACCCTTGGTTGGTTCGGGGTCGGCGGATGCCTCGCGAGCGGATTTCTGCCAGCAACGCGCTGAGACAGAGTATTGCCAGGGCCGCTTTGCGGCCCTGTCGCGACCCCAGGTTGCTCCTACGCTGGCTGGGTCGGCAGgc                                                                                             >  6:255002/1‑149 (MQ=255)
                                                      ccACACCCTTGGTTGGTTCGGGGTCGGCGGATGCCTCGCGAGCGGATTTCTGCCAGCAACGCGCTGAGACAGAGTATTGCCAGGGCCGCTTTGCGGCCCTGTCGCGACCCCAGGTTGCTCCTACGCTGGCTGGGTCGGCAGGCGCtttt                                                                                       <  4:97406/149‑1 (MQ=255)
                                                       cactgccTTTGTTGGTTCGGGGGCTGCGTATGCCTCGCGAGGGGATTTATGCCAGCAAGTCGCTGAGACAGAGTATTGCGAGTGCCGCTTTGCGGCCCTGTCGCGAGCCCAGGTTGATCCTACGCTTGCTGGGTCGGCAGGCGCtttt                                                                                       <  6:247898/143‑1 (MQ=255)
                                                        acacCCTTGGTTGGTTCGGGGTCGGCGGATGCCTCGCGAGCGGATTTCTGCCAGCAACGCGCTGAGACAGAGTATTGCCAGGGCCGCTTTGCGGCCCTGTCGCGACCCCAGGTTGCTCCTACGCTGGCTGGGTCGGCAGGCGCTTTTcc                                                                                     <  5:255002/149‑1 (MQ=255)
                                                              ttggttggTTCGGGGTCGGCGGATGCCTCGCGAGCGGATTTCTGCCAGCAACGCGCTGAGACAGAGTATTGCCAGGGCCGCTTTGCGGCCCTGTCGCGACCCCAGGTTGCTCCTACGCTGGCTGGGTCGGCAGGCGCTTTTCCATCGGc                                                                               >  8:1433/1‑149 (MQ=255)
                                                                                 cggATGCCTCGCGAGCGGATTTCTGCCAGCAACGCGCTGAGACAGAGTATTGCCAGGGCCGCTTTGCGGCCCTGTCGCGACCCCAGGTTGCTCCTACGCTGGCTGGGTCGGCAGGCGCTTTTCCATCGGCTGATACACCAGCCCGGCa                                                            >  7:182869/1‑148 (MQ=255)
                                                                                 cggATGCCTCGCGAGCGGATTTCTGCCAGCAACGCGCTGAGACAGAGTATTGCCAGGGCCGCTTTGCGGCCCTGTCGCGACCCCAGGTTGCTCCTACGCTGGCTGGGTCGGCAGGCGCTTTTCCATCGGCTGATACACCAGCCCGGCa                                                            >  2:168656/1‑148 (MQ=255)
                                                                                 cggATGCCTCGCGAGCGGATTTCTGCCAGCAACGCGCTGAGACAGAGTATTGCCAGGGCCGCTTTGCGGCCCTGTCGCGACCCCAGGTTGCTCCTACGCTGGCTGGGTCGGCAGGCGCTATTCCATCGGCTGATACACCAGCCCGGCa                                                            >  7:125452/1‑148 (MQ=255)
                                                                                        cTCGCGAGCGGATTTCTGCCAGCAACGCGCTGAGACAGAGTATTGCCAGGGCCGCTTTGCGGCCCTGTCGCGACCCCAGGTTGCTCCTACGCTGGCTGGGTCGGCAGGCGCTTTTCCATCGGCTGATACACCAGCCCGGCAAACTCAcc                                                    >  6:51473/1‑149 (MQ=255)
                                                                                         tCGCGAGCGGATTTCTGCCAGCAACGCGCTGAGACAGAGTATTGCCAGGGCCGCTTTGCGGCCCTGTCGCGACCCCAGGTTGCTCCTACGCTGGCTGGGTCGGCAg                                                                                               <  5:31447/106‑1 (MQ=255)
                                                                                         tCGCGAGCGGATTTCTGCCAGCAACGCGCTGAGACAGAGTATTGCCAGGGCCGCTTTGCGGCCCTGTCGCGACCCCAGGTTGCTCCTACGCTGGCTGGGTCGGCAg                                                                                               >  6:31447/1‑106 (MQ=255)
                                                                                                    tttCTGCCAGCAACGCGCTGAGACAGAGTATTGCCAGGGCCGCTTTGCGGCCCTGTCGCGACCCCAGGTTGCTCCTACGCTGGCTGGGTCGGCAGGCGCTTTTCCATCGGCTGATACACCAGCCCGGCAAACTCACCCACGTCGCCTGc                                        <  4:239361/149‑1 (MQ=255)
                                                                                                         gCCAGCAACGCGCTGAGACAGAGTATTGCCAGGGCCGCTTTGCGGCCCTGTCGCGACCCCAGGTTGCTCCTACGCTGGCTGGGTCGGCAGGCGCTTTTCCATCGGCTGATACACCAGCCCGGCAAACTCACCCACGTCGCCTGCCAggg                                   >  8:155860/1‑149 (MQ=255)
                                                                                                          ccAGCAACGCGCTGAGACAGAGTATTGCCAGGGCCGCTTTGCGGCCCTGTCGCGACCCCAGGTTGCTCCTACGCTGGCTGGGTCGGCAGGCGCTTTTCCATCGGCTGATACACCAGCCCGGCAAACTCACCCACGTCGCCTGCCAg                                     <  2:32912/146‑1 (MQ=255)
                                                                                                             gCAACGCGCTGAGACAGAGTATTGCCAGGGCCGCTTTGCGGCCCTGTCGCGACCCCAGGTTGCTCCTACGCTGGCTGGGTCGGCAGGCGCTTTTCCATCGGCTGATACACCAGCCCCGGCAAACTCACCCACGTCGCCTGCCAGGGTGa                                <  1:168656/149‑1 (MQ=255)
                                                                                                                          aCAGAGTATTGCCAGGGCCGCTTTGCGGCCCTGTCGCGACCCCAGGTTGCTCCTACGCTGGCTGGGTCGGCAGGCGCTTTTCCATCGGCTGATACACCAGCCCGGCAAACTCACCCACGTCGCCTGCCAGGGTGAAGCCGTAGCGTTgg                  >  2:35193/1‑149 (MQ=255)
                                                                                                                                  ttGCCAGGGCCGCTTTGCGGCCCTGTCGCGACCCCAGGTTGCTCCTACGCTGGCTGGGTCGGCAGGCGCTTTTCCATCGGCTGATACACCAGCCCGGCAAACTCACCCACGTCGCCTGCCAGGGTGAAGCCGTAGCGTTGGTAGGCgg           <  5:51473/148‑1 (MQ=255)
                                                                                                                                  ttGCCAGGGCCGCTTTGCGGCCCTGTCGCGACCCCAGGTTGCTCCTACGCTGGCTGGGTCGGCAGGCGCTTTTCCATCGGCTGATACACCAGCCCGGCAAACTCACCCACGTCGCCTGCCAGGGTGAAGCCGTAGCGTTGGTAGGCgg           <  7:155860/148‑1 (MQ=255)
                                                                                                                                    tccAGGGCCGCTTTGCGGCCCTGTCGCGTCCCGAGGTTGCTCCTTCGCTGGCTGGGTCGGGAGGCGCTTTTCCATCGGCTGATACACCAGCCCGGCAATCTCACCCACGTCGGCTGCCAGGGTGAAGCCGTAGCGTTGGTAGGCGGGt         <  8:125452/147‑1 (MQ=255)
                                                                                                                                    gacAGGGCCGCAAAGCGGCCCTGTCGCGACCCCAGGTTGCTCCTACGCTGGCTGGGTCGGCAGGCGCTTTTCCATCGGCTGATACACCAGCCCGGCAAACTCACCCACGTCGCCTGCCAGGGTGAAGCCGTAGCGTTGGTAGGCGGGTa        >  2:34671/3‑149 (MQ=255)
                                                                                                                                    gacAGGGCCGCAAAGCGGCCCTGTCGCGACCCCAGGTTGCTCCTACGCTGGCTGGGTCGGCAGGCGCTTTTCCATCGGCTGATACACCAGCCCGGCAAACTCACCCACGTCGCCTGCCAGGGTGAAGCCGTAGCGTTGGTAGGCGGGTa        >  7:19226/3‑149 (MQ=255)
                                                                                                                                     ccAGGGCCGCTTTGCGGCCCTGTCGCGACCCCAGGTTGCTCCTACGCTGGCTGGGTCGGCAGGCGCTTTTCCATCGGCTGATACACCAGCCCGGCAAACTCACCCACGTCGCCTGCCAgg                                    <  7:84655/120‑1 (MQ=255)
                                                                                                                                     ccAGGGCCGCTTTGCGGCCCTGTCGCGACCCCAGGTTGCTCCTACGCTGGCTGGGTCGGCAGGCGCTTTTCCATCGGCTGATACACCAGCCCGGCAAACTCACCCACGTCGCCTGCCAgg                                    >  8:84655/1‑120 (MQ=255)
                                                                                                                                     ccAGGGCCGCTTTGCGGCCCTGTCGCGACCCCAGGTTGCTCCTACGCTGGCTGGGTCGGCAGGCGCTTTTCCATCGGCTGATACACCAGCCCGGCAAACTCACCCACGTCGCCTGCCAGGGTGAAGCCGTAGCGTTGGTAGGCGGGTAc       <  8:182869/149‑1 (MQ=255)
                                                                                                                                          gCCGCTTTGCGGCCCTGTCGCGACCCCAGGTTGCTCCTACGCTGGCTGGGTCGGCAGGCGCTTTTCCATCGGCTGATACACCAGCCCGGCAAACTCACCCACGTCGCCTGCCAGGGTGAAGCCGTAGCGTTGGTAGGCGGGTACagaag  >  8:218924/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                
CTGTGTTCACTCGGCGCCTGAAGCGGCAACGGTTGATGGCGGGTTGACAAGCCGCCACACCCTTGGTTGGTTCGGGGTCGGCGGATGCCTCGCGAGCGGATTTCTGCCAGCAACGCGCTGAGACAGAGTATTGCCAGGGCCGCTTTGCGGCCCTGTCGCGACCCCAGGTTGCTCCTACGCTGGCTGGGTCGGCAGGCGCTTTTCCATCGGCTGATACACCAGCCCGGCAAACTCACCCACGTCGCCTGCCAGGGTGAAGCCGTAGCGTTGGTAGGCGGGTACAGAAG  >  NC_002947/3024165‑3024451

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