Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,361,447 G→A 9.7% V169M (GTG→ATG)  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,4470GA 9.7% 107.1 / 3.1 41V169M (GTG→ATG) phoPtwo‑component system DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (20/17);  new base A (2/2);  total (22/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 = 8.76e-01

GTGGTTTCACCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATCGACACGGTGCGCGGCCT  >  NC_002947/1361301‑1361589
                                                                                                                                                  |                                                                                                                                              
gTGGTTTCACCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCACCACCTGCTGGtgatgg                                                                                                                                             <  1:35674/150‑1 (MQ=255)
 tGGTTTCACCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGGCCTCAACCGCAAGCAGGCCACCCTCGACTAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGc                                                                                                                                            >  2:1157466/1‑150 (MQ=255)
    tttCACCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCaa                                                                                                                                         >  2:566673/1‑150 (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                                                                                                    <  2:982785/143‑1 (MQ=255)
                                                tCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGaa                                                                                                    >  1:982785/1‑143 (MQ=255)
                                                tCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGcgcg                                                                                             >  1:836204/1‑150 (MQ=255)
                                                tCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGcgcg                                                                                             >  2:1121111/1‑150 (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                                                                                                            >  1:402604/1‑83 (MQ=255)
                                                                                                    cGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAg                                                                                                            <  2:402604/83‑1 (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                 <  2:813649/145‑1 (MQ=255)
                                                                                                                              ttggcCACCACCTGCTGGTGATGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATCGa                 >  1:813649/4‑148 (MQ=255)
                                                                                                                                   cATCACCAGCAGGTGGTGGCCAAGGATCGCCAGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATCGACACGGTg          >  1:527458/1‑150 (MQ=255)
                                                                                                                                        ccAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCa                     <  1:838680/134‑1 (MQ=255)
                                                                                                                                        ccAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCa                     >  2:838680/1‑134 (MQ=255)
                                                                                                                                           gcagGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATCGACACGGTGCGCGGCCt  <  2:836204/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
GTGGTTTCACCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATCGACACGGTGCGCGGCCT  >  NC_002947/1361301‑1361589

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: