Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,006,777 G→A 8.5% G279G (GGC→GGT PP_5261 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,006,7770GA 8.5% 132.0 / 2.9 47G279G (GGC→GGTPP_5261LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (22/21);  new base A (2/2);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TACAGGGCATACGACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTCAAGCGCTTCTGCGCCACATCTTCGGCCACCATCAGCAGGTCGCCCATGGATATCCCGTAGCCACGCGCCGCGGCAA  >  NC_002947/6006630‑6006925
                                                                                                                                                   |                                                                                                                                                    
tACAGGGCATACGACCTGTGGGAGCGGGCGAGCCGGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCTGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCCCAACCGCTCgccgcc                                                                                                                                                    >  2:139407/1‑150 (MQ=255)
               cTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCgg                                                                                                                                     >  2:717760/1‑150 (MQ=255)
                        cGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAg                                                                                                                            <  1:592797/150‑1 (MQ=255)
                         gggCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAgg                                                                                                                           <  1:717760/150‑1 (MQ=255)
                          ggCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGgt                                                                                                                          >  2:486065/1‑150 (MQ=255)
                                cccGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAAt                                                                                                                    >  1:437988/1‑150 (MQ=255)
                                 ccGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATc                                                                                                                   >  2:618941/1‑150 (MQ=255)
                                 ccGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATc                                                                                                                   <  2:437988/150‑1 (MQ=255)
                                 ccGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATc                                                                                                                   >  1:281550/1‑150 (MQ=255)
                                                    ggTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCAccgcc                                                                                                <  2:71475/150‑1 (MQ=255)
                                                     gTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCAccgccg                                                                                               >  2:686155/1‑150 (MQ=255)
                                                        aTCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCca                                                                                                                                  >  1:313169/1‑112 (MQ=255)
                                                        aTCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCca                                                                                                                                  <  2:313169/112‑1 (MQ=255)
                                                         tCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCgg                                                                                           >  1:320221/1‑150 (MQ=255)
                                                           aaaGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTAGTTCTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGcc                                                                                         >  2:104686/1‑150 (MQ=255)
                                                                    gCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGc                                                                                <  2:281550/150‑1 (MQ=255)
                                                                              cacaTCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTCAAGCGCtt                                                                      >  2:76773/1‑150 (MQ=255)
                                                                                caTCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTCAAGCGCTTCt                                                                    <  1:76773/150‑1 (MQ=255)
                                                                                  tCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTCAAGCGCTTCTgc                                                                  <  1:422431/150‑1 (MQ=255)
                                                                                    cGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCACCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTCAAGCGCTTCTgcgc                                                                >  1:677869/1‑150 (MQ=255)
                                                                                    cGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAAAGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTACTCCATGCCGCTGGGCACCGCCGTCGGCCGTGGCAGGCTCAAGCGCTTCTgcgc                                                                >  2:265005/1‑150 (MQ=255)
                                                                                     ggCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTCAAGCGCTTCTGCGcc                                                               >  1:445599/1‑150 (MQ=255)
                                                                                      gCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTCAAGCGCTTCTGCGCca                                                              <  2:445599/150‑1 (MQ=255)
                                                                                           tCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTgg                                                                                                        <  2:638658/103‑1 (MQ=255)
                                                                                           tCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTgg                                                                                                        >  1:638658/1‑103 (MQ=255)
                                                                                           tCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTc                                                                              <  2:508036/129‑1 (MQ=255)
                                                                                           tCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTc                                                                              >  1:508036/1‑129 (MQ=255)
                                                                                           tCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTCAAGCGCTTCTGCGCCACATCt                                                         >  1:330859/1‑150 (MQ=255)
                                                                                            ccATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTCAAGCGCTTCTGCGCCCCATCtt                                                        <  1:618941/150‑1 (MQ=255)
                                                                                               ttGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTCAAGCGCTTCt                                                                    >  1:229144/1‑135 (MQ=255)
                                                                                               ttGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTCAAGCGCTTCt                                                                    <  2:229144/135‑1 (MQ=255)
                                                                                                   cgccACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTCAAGCGCTTCTGCGCCACATCTTCGGccac                                                 <  1:686155/150‑1 (MQ=255)
                                                                                                     ccACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTCAAGCGCTTCTGCGCCACATCTTCGGccacca                                               >  2:627426/1‑150 (MQ=255)
                                                                                                      cACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTCAAGCGCTTCTGCGCCACATCTTCGGCCACCAt                                              <  2:264998/150‑1 (MQ=255)
                                                                                                              cTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTCAAGCGCTTCTGCGCCACATCTTCg                                                      >  1:551753/1‑134 (MQ=255)
                                                                                                              cTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTCAAGCGCTTCTGCGCCACATCTTCg                                                      <  2:551753/134‑1 (MQ=255)
                                                                                                                  cAGGAACGCACTCAACCGCCGCAACCGCTCGCCACCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTCAAGCGCTTCTgc                                                                  >  2:80645/1‑118 (MQ=255)
                                                                                                                  cAGGAACGCACTCAACCGCCGCAACCGCTCGCCACCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTCAAGCGCTTCTgc                                                                  <  1:80645/118‑1 (MQ=255)
                                                                                                                   aGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTCAAGCGCTTCTGCGCCACATCTTCGGCCACCATCAGCAGGTc                                     <  1:661179/146‑1 (MQ=255)
                                                                                                                   aGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTCAAGCGCTTCTGCGCCACATCTTCGGCCACCATCAGCAGGTc                                     >  2:661179/1‑146 (MQ=255)
                                                                                                                      aaCGCACTCAACCGCCGCAACCGCTCGCCACCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTCAAGCGCTTCTGCGCCACATCTTCGGCCACCATCAGCAGGTCGCCCATg                              <  2:677869/150‑1 (MQ=255)
                                                                                                                        cGCACCCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGATCAAGCGCTTCTCCGCCACATCTTCGGCCACCATCAGCAGGTCGCCCATGGa                            <  1:139407/150‑1 (MQ=255)
                                                                                                                          cACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTCAAGCGCTTCTGCGCCACATCTTCGGCCACCATCAGCAGGTCGCCCATGGata                          >  2:398614/1‑150 (MQ=255)
                                                                                                                           aCTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTCAAGCGCTTCTGCGCCACATCTTCGGCCACCATCAGCAGGTCGCCCATGGatat                         >  2:571152/1‑150 (MQ=255)
                                                                                                                           aCTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTCAAGCGCTTCTGCGCCACATCTTCGGCCACCATCAGCAGGTCGCCCATGGatat                         >  1:354184/1‑150 (MQ=255)
                                                                                                                              cAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTCAAGCGCTTCTGCGCCACATCTTCGGCCACCATCAGCAGGTCGCCCATGGATATccc                      <  2:299746/150‑1 (MQ=255)
                                                                                                                                aCCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTCAAGCGCTTCTGCGCCACATCTTCGGCCACCATCAGCAGGTCGCCCATGGATATCCCGt                    <  2:330859/150‑1 (MQ=255)
                                                                                                                                  cgccgcAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTCAAGCGCTTCTGCGCCACATCTTCGGCCACCATCAGCAGGGCGCCCATGGATATCCCGTAg                  <  1:509144/150‑1 (MQ=255)
                                                                                                                                                  cgccTCGCCGGGTGCGCGGCCACAACAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTCAAGCGCTTCTGCGCCACATCTTCGGCCACCATCAGCAGGTCGCCCATGGATATCCCGTAGCCACGCGCCGCGGCaa  <  2:553514/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TACAGGGCATACGACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGCGCGGCCACACCAGGTAGTAATCCATGCCGCTGGGCACCGCCGTCGGCCATGGCAGGCTCAAGCGCTTCTGCGCCACATCTTCGGCCACCATCAGCAGGTCGCCCATGGATATCCCGTAGCCACGCGCCGCGGCAA  >  NC_002947/6006630‑6006925

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