Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,536,851 A→T 12.4% D149V (GAC→GTC)  PP_3125 → Cro/CI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,536,8510AT12.4% 110.9 / 5.5 40D149V (GAC→GTC) PP_3125Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (18/17);  new base T (3/2);  total (21/19)
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.93e-01

CGCAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCCACACGTTCTCAGGCGTGAAAC  >  NC_002947/3536703‑3536996
                                                                                                                                                    |                                                                                                                                                 
cgcAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGAc                                                                                                                                                  <  1:1005415/150‑1 (MQ=255)
 gcAGCCTGCTGCCGGTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTCTTCACACAGGCCAGGAGGAAGTGGATTTCATGACTGAGAGGGTGTTGCTGGAGCAAGGCGTCg                                                                                                                                                 >  2:142871/1‑150 (MQ=255)
    gCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGccc                                                                                                                                              <  2:1055707/150‑1 (MQ=255)
      ctgctgCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTg                                                                                                                                            >  2:354912/1‑150 (MQ=255)
           gcCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCTTAGACTCTCTTTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGTCGCCCTGctcca                                                                                                                                       <  1:142871/150‑5 (MQ=255)
                   tGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGccc                                                                                                                               >  1:986910/1‑150 (MQ=255)
                   tGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGccc                                                                                                                               >  1:657776/1‑150 (MQ=255)
                          gccgccGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGAcacc                                                                                                                        <  2:729345/150‑1 (MQ=255)
                                 tCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGTCGCCCTGc                                                                                                                                           <  1:132060/124‑1 (MQ=255)
                                 tCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGTCGCCCTGc                                                                                                                                           >  2:132060/1‑124 (MQ=255)
                                 tCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCgg                                                                                                                 <  2:1119644/150‑1 (MQ=255)
                                   cGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCt                                                                                                               >  1:1028422/1‑150 (MQ=255)
                                                ccTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGgt                                                                                                  <  1:354912/150‑1 (MQ=255)
                                                              acacCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCt                                                                                                                                             <  1:1017047/93‑1 (MQ=255)
                                                              acacCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCt                                                                                                                                             >  2:1017047/1‑93 (MQ=255)
                                                                  cAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTa                                                                                >  1:103591/1‑150 (MQ=255)
                                                                       cGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGgt                                                                           >  2:54574/1‑150 (MQ=255)
                                                                         aaGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGgttg                                                                         <  2:103591/150‑1 (MQ=255)
                                                                           gagTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGgttgtt                                                                       <  2:657776/150‑1 (MQ=255)
                                                                           gagTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGgttgtt                                                                       <  1:54574/150‑1 (MQ=255)
                                                                                 cTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGgttgttgt                                                                     <  1:241684/146‑1 (MQ=255)
                                                                                 cTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGgttgttgt                                                                     >  2:241684/1‑146 (MQ=255)
                                                                                   gtttgttCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTc                                                                   <  2:358766/146‑1 (MQ=255)
                                                                                   gtttgttCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTc                                                                   >  1:358766/1‑146 (MQ=255)
                                                                                             acaGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGaa                                                     >  1:847282/1‑150 (MQ=255)
                                                                                                    aGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAgc                                              <  2:986910/150‑1 (MQ=255)
                                                                                                     ggTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTCAGTTTCACCACGCAGACGAATGAAAgcg                                             >  1:907645/1‑150 (MQ=255)
                                                                                                         gAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAg                                         <  2:1028422/150‑1 (MQ=255)
                                                                                                           aGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAggg                                       >  2:722719/1‑150 (MQ=255)
                                                                                                             tGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAggggg                                     <  1:722719/150‑1 (MQ=255)
                                                                                                                aTTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGgttgttg                                                                      >  1:493529/1‑114 (MQ=255)
                                                                                                                aTTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGgttgttg                                                                      <  2:493529/114‑1 (MQ=255)
                                                                                                                    cATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTAt                              >  1:242751/1‑150 (MQ=255)
                                                                                                                    cATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTAt                              >  2:413330/1‑150 (MQ=255)
                                                                                                                    cATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTAt                              <  2:847282/150‑1 (MQ=255)
                                                                                                                      tgagtgagCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGc                            >  1:636213/1‑150 (MQ=255)
                                                                                                                                   gtTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGgtt                                                                          >  1:900696/1‑91 (MQ=255)
                                                                                                                                   gtTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGgtt                                                                          <  2:900696/91‑1 (MQ=255)
                                                                                                                                         aatgCAGGGCGTCGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCCACACGTTCTCAGGc         >  2:488559/4‑150 (MQ=255)
                                                                                                                                                ggCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCCACACGTTCTCAGGCGTGAAAc  >  1:287253/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CGCAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCCACACGTTCTCAGGCGTGAAAC  >  NC_002947/3536703‑3536996

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