Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,070,833 A→T 15.7% intergenic (‑133/+30) PP_5324 ← / ← phoU two‑component system response regulator/phosphate transport‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,070,8330AT15.7% 131.2 / 15.6 51intergenic (‑133/+30)PP_5324/phoUtwo‑component system response regulator/phosphate transport‑associated protein
Reads supporting (aligned to +/- strand):  ref base A (26/17);  new base T (4/4);  total (30/21)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

CTTTACTCATCGATACTCCTTTTGAATCCCGACTAGCCTACTGCCAATGATGGCCACTTGCTACAACCTGGTCACGCCGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCACCCGGAACATTAGCCTTTTCGCCGTTGGCGACAGCAGTGCCCTGCACTTCTGCCGTCATGCGCTTGAGGCCCATGTGGCGCACGTCGGTGCCGCG  >  NC_002947/6070691‑6070976
                                                                                                                                              |                                                                                                                                               
cTTTACTCATCGATACTCCTTTTGAATCCCGACTAGCCTACTGCCAATGATGGCCACTTGCTACAACCTGGTCACGCCGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCgg                                                                                                                                          <  1:1040990/150‑1 (MQ=255)
 tttACTCATCGATACTCCTTTTGAATCCCGACTAGCCTACTGCCAATGATGGCCACTTGCTACAACCTGGTCACGCCGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCggg                                                                                                                                         <  1:1401428/150‑1 (MQ=255)
  ttACTCATCGATACTCCTTTTGAATCCCGACTAGCCTACTGCCAATGATGGCCACTTGCTACAACCTGGTCACGCCGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCTTTGCCGGGc                                                                                                                                        <  2:490585/150‑1 (MQ=255)
     cTCATCGATACTCCTTTTGAATCCCGACTAGCCTACTGCCAATGATGGCCACTTGCTACAACCTGGTCACGCCGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTc                                                                                                                                     <  1:348277/150‑1 (MQ=255)
            aTACTCCTTTTGAATCCCGACTAGCCTACTGCCAATGATGGCAACTACGATGGACCTGGTCACGCCGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGACTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGAGCGTCGATGCTc                                                                                                                              >  1:1106256/1‑150 (MQ=255)
              aCTCCTTTTGAATCCCGACTAGCCTACTGCCAATGATGGCCACTTGCTACAACCTGGTCACGCCGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCgg                                                                                                                            <  1:656086/150‑1 (MQ=255)
                tCCTTTTGAATCCCGACTAGCCTACTGCCAATGATGGCCACTTGCTACAGGCTTTGCACGCCGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGa                                                                                                                          >  1:1169990/1‑150 (MQ=255)
                tCCTTTTGAATCCCGACTAGCCTACTGCCAATGATGGCCACTTGCTACAACCTGGTCACGCCGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCTTGGCCGGGCGtt                                                                                                                                     <  2:840171/139‑2 (MQ=255)
                tCCTTTTGAATCCCGACTAGCCTACTGCCAATGATGGCCACTTGCTACAACCTGGTCACGCCGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCTTGGCCGGGCGtt                                                                                                                                     >  1:840171/1‑138 (MQ=255)
                tCCTTTTGAATCCCGACTAGCCTACTGCCAATGATGGCCACTTGCTACAACCTGGTCACGCCGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGa                                                                                                                          <  2:1056836/150‑1 (MQ=255)
                   ttttGAATCCCGACTAGCCTACTGCCAATGATGGCCACTTGCTACAACCTGGTCACGCCGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCa                                                                                                                       >  1:1020476/1‑150 (MQ=255)
                   ttttGAATCCCGACTAGCCTACTGCCAATGATGGCCACTTGCTACAACCTGGTCACGCCGGGACGTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCa                                                                                                                       >  1:177473/1‑150 (MQ=255)
                                            cAATGATGGCCGCTTGCTACAACCTGGTCACGCCGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCAccc                                                                                              >  1:585891/1‑150 (MQ=255)
                                            cAATGATGGCCACTTGCTACAACCTGGTCACGCCGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCTTGGCCGGGCGt                                                                                                                                      >  2:36389/1‑110 (MQ=255)
                                            cAATGATGGCCACTTGCTACAACCTGGTCACGCCGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCTTGGCCGGGCGt                                                                                                                                      <  1:36389/110‑1 (MQ=255)
                                            cAATGATGGCCACTTGCTACAACCTGGTCACGCCGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCAccc                                                                                              >  1:1099543/1‑150 (MQ=255)
                                            cAATGATGGCCACTTGCTACAACCTGGTCACGCCGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCAccc                                                                                              >  1:170450/1‑150 (MQ=255)
                                                     ccACTTGCTACAACCTGGTCACGCCGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATc                                                                                                           >  1:1296371/1‑128 (MQ=255)
                                                     ccACTTGCTACAACCTGGTCACGCCGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATc                                                                                                           <  2:1296371/128‑1 (MQ=255)
                                                     ccACTTGCTACAACCTGGTCACGCCGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCACCCGGAACAtt                                                                                      >  2:1276731/1‑149 (MQ=255)
                                                     ccACTTGCTACAACCTGGTCACGCCGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCACCCGGAACAtt                                                                                      <  1:1276731/149‑1 (MQ=255)
                                                     ccACTTGCTACAACCTGGTCACGCCGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCACCCGGAACATTa                                                                                     >  2:136807/1‑150 (MQ=255)
                                                      cACTTGCTACAACCTGGTCACGCCGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCACCCGGAACATTAg                                                                                    >  2:811891/1‑150 (MQ=255)
                                                        cTTGCTACAACCTGGTCACGCCGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCACCCGGAACATTAGcc                                                                                  >  2:452196/1‑150 (MQ=255)
                                                            cTACAACCTGGTCACGCCGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCTTGGCCGGGCGttttc                                                                                                                                  <  1:622436/98‑5 (MQ=39)
                                                            cTACAACCTGGTCACGCCGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCTTGGCCGGGCGttttc                                                                                                                                  >  2:622436/1‑94 (MQ=39)
                                                                aaCCTGGTCACGCCGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCACCCGGAACATTAGCCTTTTCGcc                                                                          <  1:327091/150‑1 (MQ=255)
                                                                          cGCCGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCACCCGGAACATTAGCCTTTTCGCCGTTGGCGAca                                                                <  1:843297/150‑1 (MQ=255)
                                                                           gCCGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCACCCGGAAAAATAGCCTTTTCGCCGTTGGCGAcag                                                               >  2:1039420/1‑150 (MQ=255)
                                                                             cGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCACCCGGAACATTAGCCTTTTCGCCGTTGGCGAcagca                                                             <  2:177473/150‑1 (MQ=255)
                                                                                 aCTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATg                                                                                                                                 >  2:452343/1‑78 (MQ=255)
                                                                                 aCTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATg                                                                                                                                 <  1:452343/78‑1 (MQ=255)
                                                                                  cttcttCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCACCCGGAACATTAGCCTTTTCGCCGTTGGCGACAGCAGTGcc                                                        <  2:1099543/150‑1 (MQ=255)
                                                                                   ttcttcACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCACCCGGAACATTAGCCTTTTCGCCGTTGGCGACAGCAGTGccc                                                       >  2:1282602/1‑150 (MQ=255)
                                                                                   ttcttcACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCACCCGGAACATTAGCCTTTTCGCCGTTGGCGACAGCAGTGccc                                                       >  2:283326/1‑150 (MQ=255)
                                                                                   ttcttcACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCACCCGGAACATTAGCCTTTTCGCCGTTGGCGACAGCAGTGccc                                                       <  2:170450/150‑1 (MQ=255)
                                                                                     cttcACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCACCCGGAACATTAGCCTTTTCGCCGTTGGCGAcagca                                                             >  2:1151767/1‑142 (MQ=255)
                                                                                     cttcACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCACCCGGAACATTAGCCTTTTCGCCGTTGGCGAcagca                                                             <  1:1151767/142‑1 (MQ=255)
                                                                                      ttcACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGGTGCTCGGGCGCAACCTTACTTATCGTCAGATTCACCCGGAACATTAGCCTTTTCGCCGTTGGCGACAGCAGTGCCCTgg                                                    >  1:757206/1‑149 (MQ=255)
                                                                                           aTCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCACCCGGAACATCAGCCTTTTCGCCGTTGGCGCCAGCAGTGCCCTGCACTTc                                               >  1:1371931/1‑150 (MQ=255)
                                                                                            tCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCACCCGGAACATTAGCCTTTTCGCCGTTGGCGACAGCAGTGCACTGCACTTCt                                              >  1:113343/1‑150 (MQ=255)
                                                                                            ccGTTCGCGGGTGCCTGCGGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCAGGCGTCGATGCTCTGTAGCACACCTACTTATCGTCAGCTTCACCAGCACCATTCGCATATGCTCCTTTGGCGACAGCAGTGCCCTGCACTTCt                                              <  1:811891/149‑1 (MQ=255)
                                                                                                        gcctgccGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCACCCGGAACATTAGCCTTTTCGCCGTTGGCGACAGCAGTGCCCTGCACTTCTGCCGTCATgcgc                                  >  1:105083/1‑150 (MQ=255)
                                                                                                                  tctcGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCACCCGGAACATTAGCCTTTTCGCCGTTGGCGACAGCAGTGCCCTGCACTTCTGCCGTCATGCGCTTGAGGCCCa                        <  2:585891/150‑1 (MQ=255)
                                                                                                                    tcGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCACCCGGAACATTAGCCTTTTCGCCGTTGGCGACAGCAGTGCCCTGCACTTCTGCCGTCATGCGCTTGAGGCCCAtg                      <  2:1169990/150‑1 (MQ=255)
                                                                                                                     cGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCACCCGGAACATTAGCCTTTTCGCCGTTGGCGACAGCAGTGCCCTGCACTTCTGCCGTCATGCGCTTGAGGCCCAtgt                     >  2:1349155/1‑150 (MQ=255)
                                                                                                                      gACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCACCCGGAACATTAGCCTTTTCGCCGTTGGCGACAGCAGTGCCCTGCACTTCTGCCGTCATGCGCTTGAGGCCCAtgtg                    >  2:648614/1‑150 (MQ=255)
                                                                                                                           cAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCACCCGGAACATTAGCCTTTTCGCCGTTGGCGACAGCAGTGCCCTGCACTTCTGCCGTCATGCGCTTGAGGCCCATGTGGCGCa               >  2:324461/1‑150 (MQ=255)
                                                                                                                                    aCGCCCGGCCTTGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCACCCGGAACATTAGCCTTTTCGCCGTTGGCGACAGCAGTGCCCTGCACTTCTGCCGTCATGCGCTTGAGGCCCATGTGGCGCACGTCGGTGc      >  2:1126487/1‑150 (MQ=255)
                                                                                                                                        ccGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCACCCGGAACATTAGCCTTTTCGCCGTTGGCGACAGCAGTGCCCTGCACTTCTGCCGTCATGCGCTTGAGGCCCATGTGGCGCACGTCGGTGCcgcg  >  1:215112/1‑150 (MQ=255)
                                                                                                                                        ccGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCACCCGGAACATTAGCCTTTTCGCCGTTGGCGACAGCAGTGCCCTGCACTTCTGCCGTCATGCGCTTGAGGCCCATGTGGCGCACGTCGGTGCcgcg  >  1:1398663/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                               
CTTTACTCATCGATACTCCTTTTGAATCCCGACTAGCCTACTGCCAATGATGGCCACTTGCTACAACCTGGTCACGCCGGGACTTCTTCACATCGTTCGCGGGTGCCTGCCGGCTCTCGACCGCAAACGAAAACGCCCGGCCAAGGCCGGGCGTCGATGCTCGGGAGCAACCTTACTTATCGTCAGATTCACCCGGAACATTAGCCTTTTCGCCGTTGGCGACAGCAGTGCCCTGCACTTCTGCCGTCATGCGCTTGAGGCCCATGTGGCGCACGTCGGTGCCGCG  >  NC_002947/6070691‑6070976

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: