Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,361,434 T→C 11.6% H164H (CAT→CAC phoP → two‑component system DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,361,4340TC11.6% 110.1 / 6.0 43H164H (CAT→CACphoPtwo‑component system DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (21/17);  new base C (2/3);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.77e-01

CTGTTGCGCCGCTCCAGTGGTTTCACCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATCGACACGGTG  >  NC_002947/1361285‑1361581
                                                                                                                                                     |                                                                                                                                                   
cTGTTGCGCCGCTCCAGTGGTTTCACCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCAt                                                                                                                                                     >  1:1174820/1‑150 (MQ=255)
                gTGGTTTCACCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCACCACCTGCTGGtgatgg                                                                                                                                     <  1:35674/150‑1 (MQ=255)
                 tGGTTTCACCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGGCCTCAACCGCAAGCAGGCCACCCTCGACTAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGc                                                                                                                                    >  2:1157466/1‑150 (MQ=255)
                    tttCACCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCaa                                                                                                                                 >  2:566673/1‑150 (MQ=255)
                        aCCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCAc                                                                                                                                                  >  1:282833/1‑129 (MQ=255)
                        aCCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCAc                                                                                                                                                  <  2:282833/129‑1 (MQ=255)
                           cAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCCCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAGTACCTCATTCGCCATCACCAGGCGGGGGTGGCCAgggggcgc                                                                                                                          >  2:975589/1‑145 (MQ=255)
                            aGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGcc                                                                                                                         <  2:1174820/150‑1 (MQ=255)
                                    aTCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAg                                                                                                                 <  2:385529/150‑1 (MQ=255)
                                       gccgcTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAg                                                                                                              <  1:566673/150‑1 (MQ=255)
                                        ccgcTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGt                                                                                                             <  2:1082348/150‑1 (MQ=255)
                                            tGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATGGCCTAATGGAGCAGTTGTa                                                                                                         >  2:578088/1‑150 (MQ=255)
                                                  cTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAgg                                                                                                   >  1:35137/1‑150 (MQ=255)
                                                      gTGCTTGACCTCAACCGCAAGCAGGCCACGTCCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCACGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGAc                                                                                               >  1:31118/1‑150 (MQ=255)
                                                      gTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGAc                                                                                               <  2:35137/150‑1 (MQ=255)
                                                                tCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGaa                                                                                            >  1:982785/1‑143 (MQ=255)
                                                                tCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGaa                                                                                            <  2:982785/143‑1 (MQ=255)
                                                                tCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGcgcg                                                                                     >  1:836204/1‑150 (MQ=255)
                                                                tCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGcgcg                                                                                     >  2:1121111/1‑150 (MQ=255)
                                                                  aaCCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACcagca                                                                                                                                             <  1:1127122/92‑1 (MQ=255)
                                                                  aaCCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACcagca                                                                                                                                             >  2:1127122/1‑92 (MQ=255)
                                                                     cGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGc                                                                                                                          <  2:363978/108‑1 (MQ=255)
                                                                     cGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGc                                                                                                                          >  1:363978/1‑108 (MQ=255)
                                                                              gCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATc                                                                       <  1:769246/150‑1 (MQ=255)
                                                                              gCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGAAGAAGAGCGCGACCACAATGTCATc                                                                       >  1:69031/1‑150 (MQ=255)
                                                                                   ccTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGt                                                                  >  2:219678/1‑150 (MQ=255)
                                                                                       gacgaACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATc                                                                                                                            <  1:13909/88‑1 (MQ=255)
                                                                                       gacgaACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATc                                                                                                                            >  2:13909/1‑88 (MQ=255)
                                                                                       gacgaACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTGCCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTg                                                              >  1:467217/1‑150 (MQ=255)
                                                                                        acgaACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTgg                                                             <  1:219678/150‑1 (MQ=255)
                                                                                                      gCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTgcgc                                               <  1:1157466/150‑1 (MQ=255)
                                                                                                              cGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCt                                       >  1:347129/1‑150 (MQ=255)
                                                                                                                   aCGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAg                                  <  1:1121111/150‑1 (MQ=255)
                                                                                                                    cGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAg                                                                                                    <  2:402604/83‑1 (MQ=255)
                                                                                                                    cGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAg                                                                                                    >  1:402604/1‑83 (MQ=255)
                                                                                                                                  tggAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGGTGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCAGCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCa                   >  2:89368/3‑150 (MQ=255)
                                                                                                                                       tACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTGCCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAAccc              <  2:467217/150‑1 (MQ=255)
                                                                                                                                          cTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGtt                                                                                                            >  2:137011/1‑53 (MQ=255)
                                                                                                                                          cTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGtt                                                                                                            <  1:137011/53‑1 (MQ=255)
                                                                                                                                              ttggcCACCACCTGCTGGTGATGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGc                       <  1:815321/131‑1 (MQ=255)
                                                                                                                                              ttggcCACCACCTGCTGGTGATGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGc                       >  2:815321/4‑134 (MQ=255)
                                                                                                                                              ttggcCACCACCTGCTGGTGATGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATCGa         >  1:813649/4‑148 (MQ=255)
                                                                                                                                              ttggcCACCACCTGCTGGTGATGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATCGa         <  2:813649/145‑1 (MQ=255)
                                                                                                                                                   cATCACCAGCAGGTGGTGGCCAAGGATCGCCAGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATCGACACGGTg  >  1:527458/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CTGTTGCGCCGCTCCAGTGGTTTCACCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATCGACACGGTG  >  NC_002947/1361285‑1361581

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