Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,610,332 G→A 9.2% A78A (GCC→GCT)  PP_4929 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,610,3320GA 9.2% 114.5 / 2.9 44A78A (GCC→GCT) PP_4929LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (23/16);  new base A (2/2);  total (26/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

TAGCGGCTCGCCAGTTCGTCGAGGCTGTCGATCACCGGCTGATACGGCATGGCTTCGTCCTGCGCGACCCGCAGCAAGGCCTCCTGCCCGCGCATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTTCGGTTGCCGGCCGCTGCTGCGCTCGAACAGAGTGACCCCCAGGTCGGTTTCGAGCAGCGCGATGGCACTGCTGATCGCCGACTGGGCTTTGTGCTGCTCACGTGCCA  >  NC_002947/5610186‑5610477
                                                                                                                                                  |                                                                                                                                                 
tAGCGGCTCGCCAGTTCGTCGAGGCTGTCGATCACCGGCTGATACGGCATGGCTTCGTCCTGCGCGACCCGCAGCAAGGCCTCCTGCCCGCGCATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCg                                                                                                                                                >  1:779047/1‑150 (MQ=255)
     gCTCGCCAGTTCGTCGAGGCTGTCGATCACCGGCTGATACGGCATGGCTTCGTCCTGCGCGACCCGCAGCAAGGCCTCCTGCCCGCGCATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGtcctc                                                                                                                                           >  1:853052/1‑150 (MQ=255)
         gCCAGTTCGTCGAGGCTGTCGATCACCGGCTGATACGGCATGGCTTCGTCCTGCGCGACCCGCAGCAAGGCCTCCTGCCCGCGCATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGc                                                                                                                                       >  1:280222/1‑150 (MQ=255)
         gCCAGTTCGTCGAGGCTGTCGATCACCGGCTGATACGGCATGGCTTCGTCCTGCGCGACCCGCAGCAAGGCCTCCTGCCCGCGCATCAGGGCCAGGGCGCGGCCATACCGCGGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGc                                                                                                                                       >  1:851325/1‑150 (MQ=255)
           cAGTTCGTCGAGGCTGTCGATCACCGGCTGATACGGCATGGCTTCGTCCTGCGCGACCCGCAGCAAGGCCTCCTGCCCGCGCATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCaa                                                                                                                                     <  1:1121204/150‑1 (MQ=255)
            aGTTCGTCGAGGCTGTCGATCACCGGCTGATACGGCATGGCTTCGTCCTGCGCGACCCGCAGCAAGGCCTCCTGCCCGCGCATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAg                                                                                                                                    <  2:926751/150‑1 (MQ=255)
             gTTCGTCGAGGCTGTCGATCACCGGCTGATACGGCATGGCTTCGTCCTGCGCGACCCGCAGCAAGGCCTCCTGCCCGCGCATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAgg                                                                                                                                   >  1:1100053/1‑150 (MQ=255)
                 gtcgAGGCTGTCGATCACCGGCTGATACGGCATGGCTTCGTCCTGCGCGACCCGCAGCAAGGCCTCCTGCCCGCGCATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCt                                                                                                                               <  2:280222/150‑1 (MQ=255)
                         tGTCGATCACCGGCTGATACGGCATGGCTTCGTCCTGCGCGACCCGCAGCAAGGCCTCCTGCCCGCGCATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGAGCGTCCTCCAGCAAGGCGCTGCCGtttt                                                                                                                       >  1:106940/1‑150 (MQ=255)
                         tGTCGATCACCGGCTGATACGGCATGGCTTCGTCCTGCGCGACCCGCAGCAAGGCCTCCTGCCCGCGCATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGAGCGTCCTCCAGCAAGGCGCTGCCGtttt                                                                                                                       >  1:106949/1‑150 (MQ=255)
                             gATCACCGGCTGATACGGCATGGCTTCGTCCTGCGCGACCCGCAGCAAGGCCTCCTGCCCGCGCATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGt                                                                                                                   >  2:859141/1‑150 (MQ=255)
                              aTCACCGGCTGATACGGCATGGCTTCGTCCTGCGCGACCCGCAGCAAGGCCTCCTGCCCGCGCATCAGGGCCAGGGAGTTGCCATCCAGGCGTTCGCACTGGCATCACACTTGAAGTTCGTCTTGCAGCTAGGCGCTGCCGTTTTCAGTc                                                                                                                  >  1:804717/1‑150 (MQ=255)
                                  ccGGCTGATACGGCATGGCTTCGTCCTGCGCGACCCGCAGCAAGGCCTCCTGCCCGCGCATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAAtt                                                                                                              >  2:177668/1‑150 (MQ=255)
                                     gCTGATACGGCATGGCTTCGTCCTGCGCGACCCGCAGCAAGGCCTCCTGCCCGCGCATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTTCg                                                                                                           <  2:622969/150‑1 (MQ=255)
                                             ggCATGGCTTCGTCCTGCGCGACCCGCAGCAAGGCCTCCTGCCCGCGCATCAGGGCCAGGGCGCGGCCATCCAGGCGTCCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTTCGGTTgccgg                                                                                                   <  1:390346/150‑1 (MQ=255)
                                                         tCCTGCGCGACCCGCAGCAAGGCCTCCTGCCCGCGCATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTTCGGTTGCCGGCCGCTGCTgcgc                                                                                       <  2:1100053/150‑1 (MQ=255)
                                                         tCCTGCGCGACCCGCAGCAAGGCCTCCTGCCCGCGCATCAGGGCAAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCa                                                                                                                     <  1:1039/120‑1 (MQ=255)
                                                         tCCTGCGCGACCCGCAGCAAGGCCTCCTGCCCGCGCATCAGGGCAAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCa                                                                                                                     >  2:1039/1‑120 (MQ=255)
                                                              cgcgACCCGCAGCAAGGCCTCCTGCCCGCGCATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTTCGGTTGCCGGCCGCTGCTGCGCTCGaa                                                                                  >  1:89572/1‑150 (MQ=255)
                                                                    ccGCAGCAAGGCCTCCTGCCCGCGCATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTTCGGTTGCCGGCCGCTGCTGCGCTCGAACAGAGt                                                                            >  2:1101882/1‑150 (MQ=255)
                                                                      gcagcaAGGCCTCCTGCCCGCGCATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTTCGGTTGCCGGCCGCTGCTGCGCTCGAACaga                                                                              >  2:612876/1‑146 (MQ=255)
                                                                      gcagcaAGGCCTCCTGCCCGCGCATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTTCGGTTGCCGGCCGCTGCTGCGCTCGAACaga                                                                              <  1:612876/146‑1 (MQ=255)
                                                                        agcaAGGCCTCCTGCCCGCGCATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTTCGGTTGCCGGCCGCTGCTGCGCTCGAACAGAGTGAcc                                                                        >  1:587257/1‑150 (MQ=255)
                                                                         gcaAGGCCTCCTGCCCGCGCATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTTCGGTTGCCGGCCGCTGCTGCGCTCGAACAGAGTGAccc                                                                       <  2:587257/150‑1 (MQ=255)
                                                                           aaGNCCTCCTGCCCGCGCATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCAGTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTACGGATGCCGGCCGCTGCTGCGCTCGAACAGAGTGAccccc                                                                     <  1:819951/150‑1 (MQ=255)
                                                                            aGGCCTCCTGCCCGCGCATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTTCGGTTGCCGGCCGCTGCTGCGCTCGAACAGAGTGACCCCCa                                                                    >  1:479001/1‑150 (MQ=255)
                                                                              gCCTCCTGCCCGCGCATCAGGGCCAGGGCTCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTTCGGTTGCCGGCCGCTGCTGCGCTCGAACAGAGTGACCCCCAgg                                                                  >  1:426691/1‑150 (MQ=255)
                                                                              gCCTCCTGCCCGCGCATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTTCGGTTGCCGGCCGCTGCTGCGCTCGAACAGAGTGACCCCCAgg                                                                  >  2:745818/1‑150 (MQ=255)
                                                                                 tcctGCCCGCGCATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTTCGGTTGCCGGCCGCTGCTGCGCTCGAACAGAGTGACCCCCAGGTCg                                                               >  1:337786/1‑150 (MQ=255)
                                                                                    tGCCCGCGCATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTTCGGTTGCCGGCCGCTGCTGCGCTCGAACAGAGTGACCCCCAGGTCGGtt                                                            >  2:218773/1‑150 (MQ=255)
                                                                                    ggCCCGCGCATCAGGGGCAGGGCGCGGCCAAACAGACCTTCTCACTTGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTTCGGTTGCCGGCCGCTGCTGCGCTCGAACAGAGTGACCCCCAGGTCGGtt                                                            <  2:851325/149‑1 (MQ=255)
                                                                                     gCCCGCGCATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGAGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTTCGGTTGCCGGCCGCTGCTGCGCTCGAACAGAGTGACCCCCAGGTCGGttt                                                           <  2:106940/150‑1 (MQ=255)
                                                                                     gCCCGCGCATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGAGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTTCGGTTGCCGGCCGCTGCTGCGCTCGAACAGAGTGACCCCCAGGTCGGttt                                                           <  2:106949/150‑1 (MQ=255)
                                                                                      cccGCGCATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTTCGGTTGCCGGCCGCTGCGGCGCTCGAACAGAGTGACCCCCAGGTCGGTTTc                                                          <  2:53311/150‑1 (MQ=255)
                                                                                           gcATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTTCGGTTGCCGGCCGCTGCTGCGCTCGAACAGAGTGACCCCCAGGTCGGTTTCGagca                                                     <  1:129574/150‑1 (MQ=255)
                                                                                            cATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTTCGGTTGCCGGCCGCTGCTGCGCTCGAACAGAGTGACCCCCAGGTCGGTTTCGagcag                                                    >  1:567041/1‑150 (MQ=255)
                                                                                            cATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTTCGGTTGCCGGCCGCTGCTGCGCTCGAACAGAGTGACCCCCAGGTCGGTTTCGagcag                                                    <  2:779047/150‑1 (MQ=255)
                                                                                             aTCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTTCGGTTGCCGGCCGCTGCTGCGCTCGAACAGAGTGACCCCCAGGTCGGTTCCGagcagc                                                   >  1:385973/1‑150 (MQ=255)
                                                                                              tCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTTCGGTTGCCGGCCGCTGCTGCGCTCGAACAGAGTGACCCCCAGGTCGGTTTCGAGCAgcg                                                  >  2:1044399/1‑150 (MQ=255)
                                                                                                    ccagggcTCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTTCGGTTGCCGGCCGCTGCTGCGCTCGAACAGAGTGACCCCCAGGTCGGTGTCGAGCAGCGCGGTga                                            <  2:568964/150‑2 (MQ=255)
                                                                                                      agggcGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTTCGGTTGCCGGCCGCTGCTGCGCTCGAACAGAGTGACCCCCAGGTCGGTTTCGAGCAGCGCGATGGCa                                          <  1:177668/150‑1 (MQ=255)
                                                                                                                                   gCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTTCGGTTGCCGGCCGCTGCTACGCTCGAACAGAGTGACCCCCAGGTCGGTTTCGAGCAGCGCGATGGCACTGCTGATCGCCGACTGGGCTTTGtgctg             >  2:587189/1‑150 (MQ=255)
                                                                                                                                          aaTTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTTCGGTTGCCGGCCGCTGCTGCGCTCGAACAGAGTGACCCCCAGGTCGGTTTCGAGCAGTGCGATGGCACTGCTGATCGCCGACTGGGCTTTGTGCTGCTCACGt      >  1:962754/1‑150 (MQ=255)
                                                                                                                                              cgcgGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTTCGGTTGCCGGCCGCTGCTGCGCTCGAACAGAGTGACCCCCAGGTCGGTTTCGAGCAGCGCGATGGCACTGCTGATCGCCGACTGGGCTTTGTGCTGCTCACGTGCCa  >  2:273313/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TAGCGGCTCGCCAGTTCGTCGAGGCTGTCGATCACCGGCTGATACGGCATGGCTTCGTCCTGCGCGACCCGCAGCAAGGCCTCCTGCCCGCGCATCAGGGCCAGGGCGCGGCCATCCAGGCGTTCGCACTGGCGTAACAATTCGCGGGCGTCCTCCAGCAAGGCGCTGCCGTTTTCAGTCAATTTCGGTTGCCGGCCGCTGCTGCGCTCGAACAGAGTGACCCCCAGGTCGGTTTCGAGCAGCGCGATGGCACTGCTGATCGCCGACTGGGCTTTGTGCTGCTCACGTGCCA  >  NC_002947/5610186‑5610477

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