Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,812,002 C→T 10.2% R153H (CGC→CAC)  PP_1615 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,812,0020CT10.2% 104.0 / 3.1 39R153H (CGC→CAC) PP_1615LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (15/20);  new base T (2/2);  total (17/22)
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.77e-01

GTAGCGCCCACAAGTCCGAACTCCCCACCAGGCAATTGTGCCGCGCCAGCTCCGACAGGCTGTGTGGCCGGCCATAGCGCTCAAGATAGGCCGGCGAAACACACAGGTACATGCGCCTTGGCGCCAGGCGGGTAGCCACCAGGCGAGAGTCGGCCAGGCGGCCAAGACGGATGGCCAGGTCCATGCCTTCATGCAGCAGGTCGAGGGTGTGGTTGCTCAGCTCCACATCCACGCGCAACTGCGGGTACAGCGCCATGAATCGCGTCACCAGCGGCACGATGAAGCG  >  NC_002947/1811859‑1812144
                                                                                                                                               |                                                                                                                                              
gTAGCGCCCACAAGTCCGAACTCCCCACCAGGCAATTGTGCCGCGCCAGCTCCGACAGGCTGTGTGGCCGGCCATAGCGCTCAAGATAGGCCGGCGAAACACACAGGTACATGCGCCTTGGCGCCAGGCGGGTAGCCACCAGGCGAGAGt                                                                                                                                          >  2:113427/1‑150 (MQ=255)
      cccACAAGTCCGAACTCCCCACCAGGCAATTGTGCCGCGCCAGCTCCGACAGGCTGTGTGGCCGGCCATAGCGCTCAAGATAGGCCGGCGAAACACACAGGTACATGCGCCTTGGCGCCAGGCGGGTAGCCACCAGGCGAGAGTCGGCCa                                                                                                                                    >  1:241937/1‑150 (MQ=255)
           aaGTCCGAACTCCCCACCAGGCAATTGTGCCGCGCCAGCTCCGACAGGCTGTGTGGCCGGCCATAGCGCTCAAGATAGGCCGGCGAAACACACAGGTACATGCGCCTTGGCGCCAGGCGGGTAGCCACCAGGCGAGAGTCGGCCAggcgg                                                                                                                               <  1:237736/150‑1 (MQ=255)
                  aaCTCCCCACCAGGCAATTGTGCCGCGCCAGCTCCGACAGGCTGTGTGGCCGGCCATAGCGCTCAAGATAGGCCGGCGAAACACACAGGTACATGCGCCTTGGCGCCAGGCGGGTAGCCACCAGGCGAGAGTCGGCCAGGCGGCCAAGAc                                                                                                                        <  1:406673/150‑1 (MQ=255)
                      ccccACCAGGCAATTGTGCCGCGCCAGCTCGGACAGGCTGGGTGGCCGGCCATAGCGCTCAAGATAGGCCGGCGAAACACACAGGTACATGCGCCTTGGCGCCAGGCGGGTAGCCACCAGGCGAGAGTCGGCCAGGCGGCCAAGACGGAt                                                                                                                    <  1:222344/150‑1 (MQ=255)
                         caccaGGACGTTGTGCCGCGCCAGCTCCGACAGGCTGTGTGGCCGGCCATAGCGCTCAAGATAGGCCGGCGAAACACACAGGTACATGCGCCTTGGCGCCAGGCGGGTAGCCACCAGGCGAGAGTCGGCCAGGCGGCCAAGACGGATGGc                                                                                                                 <  1:306989/150‑1 (MQ=255)
                            caGGCAATTGTGCCGCGCCAGCTCCGACAGGCTGTGTGGCCGGCCATAGCGCTCAAGATAGGCCGGCGAAACACACAGGTACATGCGCCTTGGCGCCAGGCGGGTAGCCACCAGGCGAGAGTCGGCCAGGCGGCCAAGACGGATGGCCAg                                                                                                              <  1:287490/150‑1 (MQ=255)
                            caGGCAATTGTGCCGCGCCAGCTCCGACAGGCTGTGTGGCCGGCCATAGCGCTCAAGATAGGCCGGCGAAACACACAGGTACATGCGCCTTGGCGCCAGGCGCGTAGCCACCACGCGAGAGTCGGCCAGGCGGCCAAGACGAATGGCCAg                                                                                                              >  2:214318/1‑150 (MQ=255)
                                          gcgcCAGCTCCGACAGGCTGTGTGGCCGGCCATAGCGCTCAAGATAGGCCGGCGAAACACACAGGTACATGCGCCTTGGCGCCAGGCGGGTAGCCACCAGGCGAGAGTCGGCCAGGCGGCCAAGACGGATGGCCAGGTCCATGCCTTCAt                                                                                                <  2:312933/150‑1 (MQ=255)
                                             ccAGCTCCGACAGGCTGTGTGGCCGGCCATAGCGCTCAAGATAGGCCGGCGAAACACACAGGTACATGCGCCTTGGCGCCAGGCGGGTAGCCACCAGGCGAGAGTCGGCCAGGCGGCCa                                                                                                                            <  1:502339/119‑1 (MQ=255)
                                             ccAGCTCCGACAGGCTGTGTGGCCGGCCATAGCGCTCAAGATAGGCCGGCGAAACACACAGGTACATGCGCCTTGGCGCCAGGCGGGTAGCCACCAGGCGAGAGTCGGCCAGGCGGCCa                                                                                                                            >  2:502339/1‑119 (MQ=255)
                                              cAGCTCCGACAGGCTGTGTGGCCGGCCATAGCGCTCAAGATAGGCCGGCGAAACACACAGGTACATGCGCCTTGGCGCCAGGCGGGTAGCCACCAGGCGAGAGTCGGCCAGGCGGCCAAGACGGATGGCCAGGTCCATGc                                                                                                      <  1:333194/140‑1 (MQ=255)
                                              cAGCTCCGACAGGCTGTGTGGCCGGCCATAGCGCTCAAGATAGGCCGGCGAAACACACAGGTACATGCGCCTTGGCGCCAGGCGGGTAGCCACCAGGCGAGAGTCGGCCAGGCGGCCAAGACGGATGGCCAGGTCCATGc                                                                                                      >  2:333194/1‑140 (MQ=255)
                                                gCTCCGACAGGCTGTGTGGCCGGCCATAGCGCTCAAGATAGGCCGGCGAAACACACAGGTACATGCGCCTTGGCGCCAGGCGGGTAGCCACCAGGCGAGAGTCGGCCAGGCGGCCAAGACGGATGGCCAGGTCCATGCCTTCATgcagca                                                                                          >  2:245831/1‑150 (MQ=255)
                                                       cAGGCTGTGTGGCCGGCCATAGCGCTCAAGATAGGCCGGCGAAACACACAGGTACATGCGCCTTGGCGCCAGGCGGGTAGCCACCAGGCgaga                                                                                                                                            >  1:301821/1‑93 (MQ=255)
                                                       cAGGCTGTGTGGCCGGCCATAGCGCTCAAGATAGGCCGGCGAAACACACAGGTACATGCGCCTTGGCGCCAGGCGGGTAGCCACCAGGCgaga                                                                                                                                            <  2:301821/93‑1 (MQ=255)
                                                       cAGGCTGTGTGGCCGGCCATAGCGCTCAAGATAGGCCGGCGAAACACACAGGTACATGCGCCTTGGCGCCAGGCGGGTAGCCACCAGGCGAGAGTCGGCCAggcgg                                                                                                                               <  2:395199/106‑1 (MQ=255)
                                                       cAGGCTGTGTGGCCGGCCATAGCGCTCAAGATAGGCCGGCGAAACACACAGGTACATGCGCCTTGGCGCCAGGCGGGTAGCCACCAGGCGAGAGTCGGCCAggcgg                                                                                                                               >  1:395199/1‑106 (MQ=255)
                                                        aGGCTGTGTGGCCGGCCATAGCGCTCAAGATAGGCCGGCGAAACACACAGGTACATGCGCCTTGGCGCCAGGCGGGTAGCCACCAGGCGAGAGTCGGCCAGGCGGCCAAGACGGATGGCCAGGTCCATGCCTTCATGCAGCAGGTCGAgg                                                                                  <  2:241937/150‑1 (MQ=255)
                                                                tgGCCGGCCATAGCGCTCAAGATAGGCCGGCGAAACACACAGGTACATGCGCCTTGGCGCCAGGCGGGTAGCCACCAGGCGAGAGTCGGCCAGGCGGCCAAGACGGATGGCCAGGTCCATGCCTTCATGCAGCAGGTCGAGGGTGTGGtt                                                                          <  1:245831/150‑1 (MQ=255)
                                                                     ggccATAGCGCTCAAGATAGGCCGGCGAAACACACAGGTACATGGGCCTTGGCGCCAGGCGGGTAGCCACCAGGCGAGAGTCGGCCAGGCGGCCAAGACGGATGGCCAGGTCCATGCCTTCATGCAGCAGGTCGAGGGTGTGGTTGCTCa                                                                     <  2:406634/150‑1 (MQ=255)
                                                                      gccATAGCGCTCAAGATAGGCCGGCGAAACACACAGGTACATGCGCCTTGGCGCCAGGCGGGTAGCCACCAGGTGAGAGTCGGCCAGGCGGCCAAGACGGATGGCCAGGTCCATGCCTTCATGCAGCAGGTCGAGGGTGTGGTTGCTCAg                                                                    >  2:340569/1‑150 (MQ=255)
                                                                      gccATAGCGCTCAAGATAGGCCGGCGAAACACACAGGTACATGCGCCTTGGCGCCAGGCGGGTAGCCACCAGGTGAGAGTCGGCCAGGCGGCCAAGACGGATGGCCAGGTCCATGCCTTCATGCAGCAGGTCGAGGGTGTGGTTGCTCAg                                                                    >  2:340539/1‑150 (MQ=255)
                                                                              gcTCAAGATAGGCCGGCGAAACACACAGGTACATGCGCCTTGGCGCCAGGCGGGTAGCCACCAGGCGAGAGTCGGCCAGGCGGCCAAGACGGATGGCCAGGTCCATGCCTTCATGCAGCAGGTCGAGGGTGTGGTTGCTCAGCTCCACAt                                                            <  1:113427/150‑1 (MQ=255)
                                                                               cTCAAGATAGGCCGGCGAAACACACAGGTACATGCGCCTTGGCGCCAGGCGGGTAGCCACCAGGCGAGAGTCGGCCAGGCGGCCAAGACGGATGGCCAGGTCCATGCCTTCATGCAGCAGGTCGAGGGTGTGGTTGCTCAGCTCCACATc                                                           <  2:433827/150‑1 (MQ=255)
                                                                                    gATAGGCCGGCGAAACACACAGGTACATGCGCCTTGGCGCCAGGCGGGTAGCCACCAGGCGAGAGTCgg                                                                                                                                       <  2:438621/69‑1 (MQ=255)
                                                                                    gATAGGCCGGCGAAACACACAGGTACATGCGCCTTGGCGCCAGGCGGGTAGCCACCAGGCGAGAGTCgg                                                                                                                                       >  1:438621/1‑69 (MQ=255)
                                                                                               gAAACACACAGGTACATGCGCCTTGGCGCCAGGCGGGTAGCCACCAGGCGAGAGTCGGCCAGGCGGCCAAGACGGATGGCCAGGTCCATGCCTTCATGCAGCAGGTCGAGGGTGTGGTTGCTCAGCTCCACATCCACGCGCAACTGCggg                                           >  2:237183/1‑150 (MQ=255)
                                                                                                  acacacAGGTACATGCGCCTTGGCGCTAGGCGGGTAGCCACCAGGCGAGAGTCGGCCAGGCGGCCAAGACGGATGGCCAGGTCCATGCCTTCa                                                                                                 <  1:14149/93‑1 (MQ=255)
                                                                                                  acacacAGGTACATGCGCCTTGGCGCTAGGCGGGTAGCCACCAGGCGAGAGTCGGCCAGGCGGCCAAGACGGATGGCCAGGTCCATGCCTTCa                                                                                                 >  2:14149/1‑93 (MQ=255)
                                                                                                    acacaGGTACATGCGCCTTGGCGCCAGGCGGGTAGCCACCAGGCGAGAGTCGGCCAGGCGGCCAAGACGGATGGCCAGGTCCATGCCTTCATGCAGCAGGTCGAGGGTGTGGTTGCTCAGCTCCACATCCACGCGCAACTGCGGGTACAg                                      >  1:473698/1‑150 (MQ=255)
                                                                                                                gcgcCTTGGCGCCAGGCGGGTAGCCACCAGGCGAGAGTCGGCCAGGCGGCCAAGACGGATGGCCAGGTCCATGCCTTCATGCAGCAGGTCGAGGGTGTGGTTGCTCAGCTCCACATCCACGCGCAACTGCGGGTACAGCGCCATGAATcg                          >  1:330205/1‑150 (MQ=255)
                                                                                                                    cTTGGCGCCAGGCGGGTAGCCACCAGGTGAGAGTCGGCCAGGCGGCCAAGACGGATGGCCAGGTCCATGCCTTCATGCAGCAGGTCGAGGGTGTGGTTGCTCAGCTCCACATCCACGCGCAACTGCGGGTACAGCGCCATGAATCGCGTc                      <  1:340569/150‑1 (MQ=255)
                                                                                                                    cTTGGCGCCAGGCGGGTAGCCACCAGGTGAGAGTCGGCCAGGCGGCCAAGACGGATGGCCAGGTCCATGCCTTCATGCAGCAGGTCGAGGGTGTGGTTGCTCAGCTCCACATCCACGCGCAACTGCGGGTACAGCGCCATGAATCGCGTc                      <  1:340539/150‑1 (MQ=255)
                                                                                                                     ttGGCGCCAGGCGGGTAGCCACCAGGCGAGAGTCGGCCAGGCGGCCAAGACGGATGGCCAGGTCCATGCCTTCATGCAGCAGGTCGAGGGTGTGGTTGCTCAGCTCCACATCCACGCGCAACTGCGGGTACAGCGCCATGAATCGCGTCa                     <  1:237183/150‑1 (MQ=255)
                                                                                                                               gCGGGTAGCCACCAGGCGAGAGTCGGCCAGGCGGCCAAGACGGATGGCCAGGTCCATGCCTTCATGCAGCAGGTCGAGGGTGTGGTTGCTCAGCTCCACATCCACGCGCAACTGCGGGTACAGCGCCATGAATCGCGTCACCAGCGGCAc           <  2:244343/150‑1 (MQ=255)
                                                                                                                                        caccaGGCGAGAGTCGGCCAGGCGGCCAAGACGGATGGCCAGGTCCATGCCTTCATGCAGCAGGTCGAGGGTGTGGTTGCTCAGCTCCACATCCACGCGCAACTGCGGGTACAGCGCCATGAATCGCGTCACCAGCGGCACGATGAAgcg  >  2:492415/1‑150 (MQ=255)
                                                                                                                                         accaGGCGAGAGTCGGCCAGGCGGCCAAGACGGATGGCCAGGTCCATGCCTTCATGCAGCAGGTCGAGGGTGTGGTTGCTCAGCTCCACATCCACGCGCAACTGCGGGTACAGCGCCATGAATc                           <  1:171142/124‑1 (MQ=255)
                                                                                                                                         accaGGCGAGAGTCGGCCAGGCGGCCAAGACGGATGGCCAGGTCCATGCCTTCATGCAGCAGGTCGAGGGTGTGGTTGCTCAGCTCCACATCCACGCGCAACTGCGGGTACAGCGCCATGAATc                           >  2:171142/1‑124 (MQ=255)
                                                                                                                                               |                                                                                                                                              
GTAGCGCCCACAAGTCCGAACTCCCCACCAGGCAATTGTGCCGCGCCAGCTCCGACAGGCTGTGTGGCCGGCCATAGCGCTCAAGATAGGCCGGCGAAACACACAGGTACATGCGCCTTGGCGCCAGGCGGGTAGCCACCAGGCGAGAGTCGGCCAGGCGGCCAAGACGGATGGCCAGGTCCATGCCTTCATGCAGCAGGTCGAGGGTGTGGTTGCTCAGCTCCACATCCACGCGCAACTGCGGGTACAGCGCCATGAATCGCGTCACCAGCGGCACGATGAAGCG  >  NC_002947/1811859‑1812144

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: