Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,127,218 C→G 11.8% V241V (GTG→GTC) ‡ PP_5375 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,127,2180CG11.8% 90.5 / 5.5 34V241V (GTG→GTC) ‡PP_5375LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (15/15);  new base G (2/2);  total (17/17)
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.82e-01

CGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTGGGAGTCTT  >  NC_002947/6127069‑6127354
                                                                                                                                                     |                                                                                                                                        
cGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGcc                                                                                                                                          >  2:16837/1‑150 (MQ=255)
         cGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTgg                                                                                                                                 >  1:518410/1‑150 (MQ=255)
              gccttCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCGTGGCCCTGGcctgg                                                                                                                            <  2:786835/147‑4 (MQ=255)
                tttCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTTCCCCAGGCCAGGGCCACGCCCTGGCCAATGa                                                                                                                          >  2:344823/1‑150 (MQ=255)
                                       gacgacTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCGGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCg                                                                                                   <  1:16837/150‑1 (MQ=255)
                                       gacgacTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCGTGGCCCTGGc                                                                                                                                >  2:152955/1‑121 (MQ=255)
                                         cgacTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGAAGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTa                                                                                                 <  1:481207/150‑1 (MQ=255)
                                                agcagGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCGTGGCCCTGGcc                                                                                                                               <  2:767719/113‑1 (MQ=255)
                                                agcagGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCGTGGCCCTGGcc                                                                                                                               >  1:767719/1‑113 (MQ=255)
                                                agcagGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTg                                                                                          >  1:647178/1‑150 (MQ=255)
                                                    ggcagAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGatat                                                                                      <  1:167590/150‑1 (MQ=255)
                                                        gAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTCGCTGTTGATATTGAc                                                                                  >  2:595092/1‑150 (MQ=255)
                                                          aCTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGc                                                                                                                      >  1:557744/1‑112 (MQ=255)
                                                          aCTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGc                                                                                                                      <  2:557744/112‑1 (MQ=255)
                                                              ggCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATg                                                                            <  1:286259/150‑1 (MQ=255)
                                                                       cTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGt                                                                   <  1:709196/150‑1 (MQ=255)
                                                                         ggAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTggg                                                                      <  1:482317/145‑1 (MQ=255)
                                                                         ggAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTggg                                                                      >  2:482317/1‑145 (MQ=255)
                                                                                cacaTCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTc                                                          >  1:290587/1‑150 (MQ=255)
                                                                                  caTCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTgg                                                                       <  2:561100/135‑1 (MQ=255)
                                                                                  caTCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTgg                                                                       >  1:561100/1‑135 (MQ=255)
                                                                                      gTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTc                                                                  <  2:87948/136‑1 (MQ=255)
                                                                                      gTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTc                                                                  >  1:87948/1‑136 (MQ=255)
                                                                                         gggCTTGACCAGCGCCCCGGTACTCAACGCATCGTCCTCCAGGTTGCCCCAGGCCAGGGCCACCCCTTGGACACTGCTCGCCTACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAg                                                 <  2:517754/150‑1 (MQ=255)
                                                                                         gggCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGAcgc                                                                                >  2:345760/1‑119 (MQ=255)
                                                                                         gggCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGAcgc                                                                                <  1:345760/119‑1 (MQ=255)
                                                                                           gCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGgc                                               <  1:595092/150‑1 (MQ=255)
                                                                                                            gTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGc                              >  2:671910/1‑150 (MQ=255)
                                                                                                                  aaCGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCt                        >  1:566228/1‑150 (MQ=255)
                                                                                                                               ccaGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGCTGCTTGGTTGTCATGGGCTc                                                          <  2:635264/103‑1 (MQ=255)
                                                                                                                               ccaGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTc                                                          >  1:635264/1‑103 (MQ=255)
                                                                                                                                      gCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTGGGAGTc    <  2:290587/150‑1 (MQ=255)
                                                                                                                                        cccAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGACCCGTATCACCGACGACTGACGCATTGGTACCAACCGACCCAGTCGGGCTGGGAGTCtt  >  2:8213/1‑150 (MQ=255)
                                                                                                                                           aggccaggGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGAcgc                                                                                <  2:319162/69‑1 (MQ=255)
                                                                                                                                           aggccaggGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGAcgc                                                                                >  1:319162/1‑69 (MQ=255)
                                                                                                                                                     |                                                                                                                                        
CGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTGGGAGTCTT  >  NC_002947/6127069‑6127354

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: