Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,037,488 T→C 8.3% A33A (GCA→GCG)  PP_2649 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,037,4880TC 8.3% 129.1 / 2.7 48A33A (GCA→GCG) PP_2649LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (22/22);  new base C (2/2);  total (24/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

CAGGCGTGCGTAGCGGATGGCGCAGTGCCCGGCGGCCGTGGGCTCCAGGCCTCGGTTGGTGCGGTTGAACAGCACGGTGCCGAAGGTGGTTTCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGTGGCGGGCGTGCAGGCGGGACATGATCGAATTGAGCGATGGGGTCATCTGGTATTCGTCCTGGTTATACCGATATGGATTGTTCT  >  NC_002947/3037344‑3037624
                                                                                                                                                 |                                                                                                                                        
cAGGCGTGCGTAGCGGATGGCGCAGTGCCCGGCGGCCGTGGGCTCCAGGCCTCGGTTGGTGCGGTTGAACAGCACGGTGCCGAAGGTGGTTTCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTtgctg                                                                                                                                      >  7:121486/1‑149 (MQ=255)
   gcgtgcgtAGCGGATGGCGCAGTGCCCGGCGGCCGTGGGCTCCAGGCCTCGGTTGGTGCGGTTGAACAGCACGGTGCCGAAGGTGGTTTCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCAc                                                                                                                                   >  5:224605/1‑149 (MQ=255)
       gcgtAGCGGATGGCGCAGTGCCCGGCGGCCGTGGGCTCCAGGCCTCGGTTGGTGCGGTTGAACAGCACGGTGCCGAAGGTGGTTTCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGCGTAAGCCCGACG‑CTGCTTGCTGCACCGag                                                                                                                               >  8:39953/1‑149 (MQ=255)
              ggATGGCGCAGTGCCCGGCGGCCGTGGGCTCCAGGCCTCGGTTGGTGCGGTTGAACAGCACGGTGCCGAAGGTGGTTTCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCagcgag                                                                                                                        <  6:224605/149‑1 (MQ=255)
                aTGGCGCAGTGCCCGGCGGCCGTGGGCTCCAGGCCTCGGTTGGTGCGGTTGAACAGCACGGTGCCGAAGGTGGTTTCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCagcgagcg                                                                                                                      >  3:204153/1‑149 (MQ=255)
                   gcgcAGTGCCCGGCGGCCGTGGGCTCCAGGCCTCGGTTGGTGCGGTTGAACAGCACGGTGCCGAAGGTGGTTTCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTg                                                                                                                   <  4:165912/149‑1 (MQ=255)
                    cgcAGTGCCCGGCGGCCGTGGGCTCCAGGCCTCGGTTGGTGCGGTTGAACAGCACGGTGCCGAAGGTGGTTTCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTCGCTGCTCGGTGCagcaagca                                                                                                                      >  6:105713/1‑144 (MQ=255)
                    cgcAGTGCCCGGCGGCCGTGGGCTCCAGGCCTCGGTTGGTGCGGTTGAACAGCACGGTGCCGAAGGTGGTTTCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTCGCTGCTCGGTGCagcaagca                                                                                                                      <  5:105713/145‑2 (MQ=255)
                     gcAGTGCCCGGCGGCCGTGGGCTCCAGGCCTCGGTTGGTGCGGTTGAACAGCACGGTGCCGAAGGTGGTTTCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCt                                                                                                                 >  6:269483/1‑149 (MQ=255)
                      cAGTGCCCGGCGGCCGTGGGCTCCAGGCCTCGGTTGGTGCGGTTGAACAGCACGGTGCCGAAGGTGGTTTCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTc                                                                                                                >  2:75138/1‑149 (MQ=255)
                            ccGGCGGCCGTGGGCTCCAGGCCTCGGTTGGTGCGGTTGAACATCACGGTGCCGAAGGTGGTTTCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCaa                                                                                                          <  8:189979/149‑1 (MQ=255)
                             cggcggCCGTGGGCTCCAGGCCTCGGTTGGTGCGGTTGAACAGCACGGTGCCGAAGGTGGTTTCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACGCCTGCTCGCTGCTCGGTGCagcaagcag                                                                                                                     <  1:203789/138‑3 (MQ=255)
                             cggcggCCGTGGGCTCCAGGCCTCGGTTGGTGCGGTTGAACAGCACGGTGCCGAAGGTGGTTTCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTCGCTGCTCGGTGCagcaagcag                                                                                                                     >  2:203789/1‑135 (MQ=255)
                                 ggcCGTGGGCTCCAGGCCTCGGTTGGTGCGGTTGAACAGCACGGTGCCGAAGGTGGTTTCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGa                                                                                                     <  3:31720/149‑1 (MQ=255)
                                   ccGTGGGCTCCAGGCCTCGGTTGGTGCGGTTGAACAGCACGGTGCCGAAGGTGGTTTCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATc                                                                                                   <  8:121486/149‑1 (MQ=255)
                                        ggCTCCAGGCCTCGGTTGGTGCGGTTGAACAGCACGGTGCCGAAGGTGGTTTCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATcaac                                                                                                >  2:44411/1‑147 (MQ=255)
                                         gCTCCAGGCCTCGGTTGGTGCGGTTGAACAGCACGGTGCCGAAGGTGGTTTCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAg                                                                                             >  8:35498/1‑149 (MQ=255)
                                               ggCCTCGGTTGGTGCGGTTGAACAGCACGGTGCCGAAGGTGGTTTCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAgg                                                                                       >  2:365548/1‑149 (MQ=255)
                                                      gTTGGTGCGGTTGAACAGCACGGTGCCGAAGGTGGTTTCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGTggcggg                                                                                >  4:175234/1‑149 (MQ=255)
                                                         ggTGCGGTTGAACAGCACGGTGCCGAAGGTGGTTTCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGTGGCGGGCGt                                                                             >  4:139216/1‑149 (MQ=255)
                                                              ggTTGAACAGCACGGTGCCGAAGGTGGTTTCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGTGGCGGGCGTGCAgg                                                                        <  3:279611/149‑1 (MQ=255)
                                                               gTTGAACAGCACGGTGCCGAAGGTGGTTTCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGTGGCGGGCGTGCAGGc                                                                       >  7:307297/1‑149 (MQ=255)
                                                                  gAACAGCACGGTGCCGAAGGTGGTTTCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGTGGCGGGCGTGCAGGCgg                                                                     <  3:139216/148‑1 (MQ=255)
                                                                               ccGAAGGTGGTTTCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGTGGCGGGCGTGCAGGCGGGACATGATCGAAtt                                                       <  7:35498/149‑1 (MQ=255)
                                                                                  aaGGTGGTTTCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGTGGCGGGCGTGCAGGCGGGACATGATCGAATTGAg                                                    >  5:269472/1‑149 (MQ=255)
                                                                                        gtTTCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGTGGCGGGCGTGCAGGCGGGACATGATCGAATTGa                                                     <  7:88181/142‑1 (MQ=255)
                                                                                        gtTTCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGTGGCGGGCGTGCAGGCGGGACATGATCGAATTGa                                                     >  8:88181/1‑142 (MQ=255)
                                                                                        gcttCGATTTGCTGCAAGGCCTTGCCGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGTGGCGGGCGTGCAGGCGGGACATGATCGAATTGAGCGATgg                                              <  4:189534/147‑1 (MQ=255)
                                                                                         tttCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGTGGCGGGCGTGCAGGCGGGACATGATCGAATTGAGCGATggg                                             <  7:39953/149‑1 (MQ=255)
                                                                                         tttCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGGGGCGGGCGTGCAGGCGGGACATGATCGAATTGAGCGATggg                                             <  4:204153/149‑1 (MQ=255)
                                                                                             gATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGTGGCGGGCGTGCAGGCGGGACATGATCGAATTGAGCGATGGGGTCa                                         >  8:39068/1‑149 (MQ=255)
                                                                                                    tGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGTGGCGGGCGTGCAGGCGGGACATGATCGAATTGAGCGATg                                               <  3:207038/136‑1 (MQ=255)
                                                                                                    tGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGTGGCGGGCGTGCAGGCGGGACATGATCGAATTGAGCGATg                                               >  4:207038/1‑136 (MQ=255)
                                                                                                    tGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGTGGCGGGCGTGCAGGCGGGACATGATCGAATTGAGCGATGGGGTCATCTGGTa                                  >  1:236794/1‑149 (MQ=255)
                                                                                                     gCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGTGGCGGGCGTGCAGGCGGGACATGATCGAATTGAGCGATGGGGTCATCTGGTAt                                 >  8:145886/1‑149 (MQ=255)
                                                                                                               gCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGTGGCGGGCGTGCAGGCGGGACATGATCGAATTGAGCGATGGGGTCATCTGGTATTCGTCCTg                         <  8:307297/147‑1 (MQ=255)
                                                                                                                    cACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGTGGCGGGCGTGCAGGCGGGACATGATCGAATTGAGCGATGGGGTCATCTGGTATTCGTCCTGGTTATAc                  <  7:145886/149‑1 (MQ=255)
                                                                                                                    cACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGTGGCGGGCGTGCAGGCGGGACATGATCGAATTGAGCGATGGGGTCATCTGGTATTCGTCCTGGTTATAc                  <  1:365548/149‑1 (MQ=255)
                                                                                                                       aaGGCTGTGTAAGCCCGACG‑CTGCTTGCTGAACAGAGCAGCGAGCGGTGCTAGTCCAACTAGATCAACAAGCGAAGGTGGCGGGCGTGCAGGCGGGACATGATCGAATTGAGAGATGGGGTCATATTGTATTCGTCCTGGTtata                   >  3:114952/2‑145 (MQ=255)
                                                                                                                          gCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGTGGCGGGCGTGCAGGCGGGACATGATCGAATTGAGCGATgggg                                            <  8:208338/117‑1 (MQ=255)
                                                                                                                          gCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGTGGCGGGCGTGCAGGCGGGACATGATCGAATTGAGCGATgggg                                            >  7:208338/1‑117 (MQ=255)
                                                                                                                          gCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGTGGCGGGCGTGCAGGCGGGACATGATCGAATTGAGCGATGGGGTCATCTGGTATTCGTCCTGGTTATACCGatat            <  5:269483/149‑1 (MQ=255)
                                                                                                                                  agaccGACG‑CTGCTTGCTGCACCGAGCAGCGAGCTTTTCTCGTACAACTCGATCAACAAGCGCAGTTGGCGGGCTTGCATTCTGGACATTATAGAATTGAGCGATTGTGTCATCTGTTATTCGTCCTTGTTATACCGATATGGATtttt    <  4:114952/146‑1 (MQ=255)
                                                                                                                                    cccGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGTGGCGGGCGTGCAGGCGGGACATGATCGAATTGAGCGATGGGGTCATCTGGTATTCGTCCTGGTTATACCGATATGGATTGTtct  <  1:75138/149‑1 (MQ=255)
                                                                                                                                     ccGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGTGGCGGGCGTGCAGGCGGGACATGATCGAATTGAGCGATGGGGTCATCTGGTATTCGTCCTGGTTATACCGATATGGATTGtt    <  7:311127/146‑1 (MQ=255)
                                                                                                                                     ccGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGTGGCGGGCGTGCAGGCGGGACATGATCGAATTGAGCGATGGGGTCATCTGGTATTCGTCCTGGTTATACCGATATGGATTGtt    >  8:311127/1‑146 (MQ=255)
                                                                                                                                     ccGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGTGGCGGGCGTGCAGGCGGGACATGATCGAATTGAGCGATGGGGTCATCTGGTATTCGTCCTGGTTATACCGATATGGATTGTtct  >  1:168383/1‑148 (MQ=255)
                                                                                                                                          g‑ctgctTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGTggcgggc                                                                               >  8:29477/1‑66 (MQ=255)
                                                                                                                                          g‑ctgctTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGTggcgggc                                                                               <  7:29477/66‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                        
CAGGCGTGCGTAGCGGATGGCGCAGTGCCCGGCGGCCGTGGGCTCCAGGCCTCGGTTGGTGCGGTTGAACAGCACGGTGCCGAAGGTGGTTTCGATTTCCTGCAAGGCCTTGCTGGCACCAGGCTGGGTAAGCCCGACG‑CTGCTTGCTGCACCGAGCAGCGAGCGGTGCTCGTCCAACTCGATCAACAAGCGCAGGTGGCGGGCGTGCAGGCGGGACATGATCGAATTGAGCGATGGGGTCATCTGGTATTCGTCCTGGTTATACCGATATGGATTGTTCT  >  NC_002947/3037344‑3037624

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