Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,897,725 T→A 13.4% intergenic (+99/‑51) PP_2551 → / → PP_2552 LysR family transcriptional regulator/DOPA decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,897,7250TA13.4% 101.8 / 5.5 43intergenic (+99/‑51)PP_2551/PP_2552LysR family transcriptional regulator/DOPA decarboxylase
Reads supporting (aligned to +/- strand):  ref base T (18/19);  new base A (4/2);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GCGCCAGTCGCAGTTCGGGGTGCAGGGCGGGCATGCGGGTTCGCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGG‑CCACCAACTGATCGACCTGATTGCCGACTACCGCCAGACCGTGGGCGAACGCCCGGTCATGGCC  >  NC_002947/2897578‑2897868
                                                                                                                                                   |                                                                                                                                                
gcgcCAGTCGCAGTTCGGGGTGCAGGGCGGGCATGCGGGTTCGCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTa                                                                                                                                                 >  4:276854/1‑149 (MQ=255)
                  ggTGCAGGGCGGGCATGCGGGTTCGCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGca                                                                                                                               >  8:108731/1‑149 (MQ=255)
                      cAGGGCGGGCATGCGGGTTCGCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAgg                                                                                                                           <  2:97311/149‑1 (MQ=255)
                          gcgggcATGCGGGTTCGCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAAt                                                                                                                       <  2:133605/149‑1 (MQ=255)
                               cATGCGGGTTCGCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGAt                                                                                                                   >  6:51315/1‑148 (MQ=255)
                                       ttCGCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAg                                                                                                          >  4:63002/1‑149 (MQ=255)
                                          gCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGga                                                                                                        <  5:315347/148‑1 (MQ=255)
                                            cATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGagag                                                                                                     >  6:325526/1‑149 (MQ=255)
                                              tAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCt                                                                                                   <  1:260560/149‑1 (MQ=255)
                                              tAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTCGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCt                                                                                                   <  3:195661/149‑1 (MQ=255)
                                                 cagcTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGc                                                                                                <  7:108731/149‑1 (MQ=255)
                                                              aGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAAc                                                                                   <  2:204290/149‑1 (MQ=255)
                                                              aGCGAGGTTATTCGGCTGGGGATAGGGGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGaa                                                                                    <  8:46833/148‑1 (MQ=255)
                                                                    gTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTcc                                                                             <  3:63002/149‑1 (MQ=255)
                                                                    gTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTAAACCcgc                                                                                                                                            <  2:266074/86‑1 (MQ=255)
                                                                    gTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGAATGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTAAACCcgc                                                                                                                                            >  1:266074/1‑86 (MQ=255)
                                                                       aTTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGc                                                                           <  8:112873/148‑1 (MQ=255)
                                                                       aTTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCAGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCGAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGACCAATTCCGcc                                                                          >  8:105504/1‑149 (MQ=255)
                                                                         tCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAg                                                                        >  6:42187/1‑149 (MQ=255)
                                                                            gCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTAc                                                                     >  3:43334/1‑149 (MQ=255)
                                                                                     aGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGG‑CCACCaa                                                           >  6:90171/1‑149 (MQ=255)
                                                                                      ggTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAAATCCGCCAGTACGG‑CCACCAAc                                                          <  4:43334/149‑1 (MQ=255)
                                                                                       gTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGG‑CCACCAACt                                                         >  5:343701/1‑149 (MQ=255)
                                                                                           cTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCACCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGG‑CCACCAACTGATc                                                     <  5:51315/149‑1 (MQ=255)
                                                                                              gACTGGGGCAGTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGG‑CCACCAACTGATCGAc                                                  <  2:287155/149‑1 (MQ=255)
                                                                                                cTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGG‑CCACCAACTGATCGACCt                                                >  4:328785/1‑149 (MQ=255)
                                                                                                   gggCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGG‑CCACCAACTGATCGACCTGAt                                             >  2:7896/1‑149 (MQ=255)
                                                                                                               gATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCTCAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGG‑CCACCAACTGATCGACCTGATTGCCGACTAcc                                  <  4:339209/148‑1 (MQ=255)
                                                                                                               gATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCCACCAACTGATCGACCTGATTGCCGACTa                                    <  1:22162/147‑1 (MQ=255)
                                                                                                               gATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGG‑CCACCAACTGATCGACCTGATTGCCGACTa                                    >  2:22162/1‑146 (MQ=255)
                                                                                                                    tgCGGCTTCTTCGCGGCTGTAGGCGCAGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGG‑CCACCAACTGATCGACCTGATTGCCGACTACCGCCAGa                            >  5:159277/1‑149 (MQ=255)
                                                                                                                          ttcttcGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGG‑CCACCAACTGATCGACCTGATTGCCGACTACCGCCAGACCGTgg                      >  1:147856/1‑149 (MQ=255)
                                                                                                                             ttcGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGG‑CCTCCAACTGATCGACCTCATTTCCGACTACCGCCAGACCGTGGGc                    >  7:184623/1‑148 (MQ=255)
                                                                                                                                       tAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGG‑CCACCAACTGATCGACCTGATTGCCGACTACCGCCAGACCGTGGGCGAACGCCCg           <  5:325526/147‑1 (MQ=255)
                                                                                                                                       tAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGG‑CCACCAACTGATCGACCTGATTGCCGACTACCGCCAGACCGTGGGCGAACGCCCg           <  5:90171/147‑1 (MQ=255)
                                                                                                                                        tttcgcgGGTAAACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGG‑CCACCAACTGATCGACCTGATTGCCGACTACCGCCAg                             >  1:354743/4‑128 (MQ=255)
                                                                                                                                        tttcgcgGGTAAACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGG‑CCACCAACTGATCGACCTGATTGCCGACTACCGCCAg                             <  2:354743/125‑1 (MQ=255)
                                                                                                                                        tttcgcgGGTAAACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGG‑CCACCAACTGATCGACCTGATTGCCGACTACCGCCAGACCGTGGGCGAACGCCCGGt         >  3:78663/4‑148 (MQ=255)
                                                                                                                                        tttcgcgGGTAAACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGG‑CCACCAACTGATCGACCTGATTGCCGACTACCGCCAGACCGTGGGCGAACGCCCGGt         >  2:160208/4‑148 (MQ=255)
                                                                                                                                         ggCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGG‑CCACCAACTGATCGACCTGATTGCCGACTACCGCCAGACCGTTGGCGAACGCCCGGTc        <  6:343701/148‑1 (MQ=255)
                                                                                                                                            gcgGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGG‑CCACCAACTGATCGACCTGATTGCCGACTACCGCCAGACCGTGGGCGAACGCCCGGTCATgg    <  3:328785/149‑1 (MQ=255)
                                                                                                                                            gcgGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCCTGACCCCCGAACAATTCCGCCAGTACGG‑CCACCAACTGATCGACCTGATTGCCGACTACCGCCAGACCGTGGGCGAACGCCCGGTCATgg    >  3:258534/1‑149 (MQ=255)
                                                                                                                                                gTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGG‑CCACCAACTGATCGACCTGATTGCCGCCTCCCGCCAGACCGTGGGCGAACGCCCGGTCATGGcc  >  4:54769/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GCGCCAGTCGCAGTTCGGGGTGCAGGGCGGGCATGCGGGTTCGCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGG‑CCACCAACTGATCGACCTGATTGCCGACTACCGCCAGACCGTGGGCGAACGCCCGGTCATGGCC  >  NC_002947/2897578‑2897868

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