Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,685,134 A→C 8.2% intergenic (‑12/+220) PP_1480 ← / ← patD hypothetical protein/medium chain aldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,685,1340AC 8.2% 144.9 / 3.4 49intergenic (‑12/+220)PP_1480/patDhypothetical protein/medium chain aldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (21/24);  new base C (2/2);  total (23/26)
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.13e-01

GCACCCGCGATCCACGGCCTTCGATCGCCTTGGCGCCTATCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGCGAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATGGCGGTAGTCGTGTTGTGTGGCGAAAATCGAGCGCCGCCAGCGCG  >  NC_002947/1684987‑1685278
                                                                                                                                                   |                                                                                                                                                
gcACCCGCGATCCACGGCCTTCGATCGCCTTGGCGCCTATCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAatat                                                                                                                                                 <  4:236280/149‑1 (MQ=255)
gcACCCGCGATCCACGGCCTTCGATCGCCTTGGCGCCTATCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAatat                                                                                                                                                 <  2:70411/149‑1 (MQ=255)
  aCCCGCGATCCACGGCCTTCGATCGCCTTGGCGCCTATCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGa                                                                                                                                               >  8:280425/1‑149 (MQ=255)
    ccGCGATCCACGGCCTTCGATCGCCTTGGCGCCTATCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATa                                                                                                                                             >  2:44983/1‑149 (MQ=255)
       cgATCCACGGCCTTCGATCGCCTTGGCGCCTATCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCATATTCTGATACCa                                                                                                                                          <  1:187673/149‑1 (MQ=255)
            cACGGCCTTCGATCGCCTTGGCGCCTATCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTg                                                                                                                                      <  4:214957/148‑1 (MQ=255)
             aCGGCCTTCGATCGCCTTGGCGCCTATCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGAt                                                                                                                                    >  7:317793/1‑149 (MQ=255)
                 ccTTCGATCGCCTTGGCGCCTATCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGAttgtt                                                                                                                                <  2:71937/149‑1 (MQ=255)
                     cGATCGCCTTGGCGCCTATCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGt                                                                                                                             <  2:317653/148‑1 (MQ=255)
                      gATCGCCTTGGCGCCTATCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTAt                                                                                                                           <  3:150798/149‑1 (MQ=255)
                          gCCTTGGCGCCTATCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAg                                                                                                                       >  5:368044/1‑149 (MQ=255)
                                    cTATCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGa                                                                                                             <  8:267774/149‑1 (MQ=255)
                                     tATCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAg                                                                                                            >  2:164125/1‑149 (MQ=255)
                                           tatCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTAt                                                                                                      <  4:333857/149‑1 (MQ=255)
                                                gAAATGACTCCAAACGCGGCCACTCCAGCGACCTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTg                                                                                                   >  5:312482/1‑147 (MQ=255)
                                                 aaaTGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAg                                                                                                <  2:305570/149‑1 (MQ=255)
                                                 aaaTGACTCCAAACGCCGCCACTCCAGCGACTTTGGAACCGCCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGAttgttg                                                                                                                               <  8:341726/118‑1 (MQ=255)
                                                       cTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTgg                                                                                          >  1:163823/1‑149 (MQ=255)
                                                       cTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTgg                                                                                          <  1:44983/149‑1 (MQ=255)
                                                        tCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTTTTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGt                                                                                         >  8:308848/1‑149 (MQ=255)
                                                        tCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGt                                                                                            <  8:4919/146‑1 (MQ=255)
                                                        tCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGt                                                                                            >  7:4919/1‑146 (MQ=255)
                                                             aCGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGc                                                                                    <  2:163823/149‑1 (MQ=255)
                                                             aCGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGc                                                                                    <  7:280425/149‑1 (MQ=255)
                                                               gcgGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGAttgtt                                                                                                                                >  3:90113/1‑103 (MQ=255)
                                                               gcgGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGAttgtt                                                                                                                                <  4:90113/103‑1 (MQ=255)
                                                                   ccACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGCGAGCTa                                                                              <  5:143017/149‑1 (MQ=255)
                                                                      cTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGCGAGCTAggg                                                                           <  8:48909/149‑1 (MQ=255)
                                                                      cTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGCGAGCTAggg                                                                           <  8:317793/149‑1 (MQ=255)
                                                                       tCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGCGAGCTAgggg                                                                          <  7:77255/149‑1 (MQ=255)
                                                                         cAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGCGAGCTAgggg                                                                          <  1:350140/147‑1 (MQ=255)
                                                                          aGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGCGAGCTAGGGGTGc                                                                       >  4:172853/1‑149 (MQ=255)
                                                                                    gAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGCGAGCTAGGGGTGCTGGGCGTCGa                                                             <  6:368044/149‑1 (MQ=255)
                                                                                        cGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGCGAGCTAGGGGTGCTGGGCGTCGACGtt                                                         <  1:164125/149‑1 (MQ=255)
                                                                                                   aaaTGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACAAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGCGAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATg                                              >  1:101328/1‑149 (MQ=255)
                                                                                                          ttCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTg                                                                                                                                      <  3:217913/54‑1 (MQ=255)
                                                                                                          ttCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGCGAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATGGCGGTAg                                       >  5:341374/1‑149 (MQ=255)
                                                                                                          ttCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGCGAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATGGCGGTAg                                       >  3:117036/1‑149 (MQ=255)
                                                                                                          ttCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGCGAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATGGCGGTAg                                       >  8:253354/1‑149 (MQ=255)
                                                                                                            cAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGCGAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATGGCGGTAg                                       <  5:252168/147‑1 (MQ=255)
                                                                                                             aaaGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGCGAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATGGCGGTAGTc                                     <  3:47393/148‑1 (MQ=255)
                                                                                                                        attGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGAATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGCGAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATGGCGGTAGTCGTGTTGTGTGGc                         >  6:76686/2‑149 (MQ=255)
                                                                                                                              gatTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGCGAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATGGCGGTAGTCGTGTTGTGTGGCGAAAAt                   >  8:266219/1‑149 (MQ=255)
                                                                                                                                   ttCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGCGAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATGGCGGTAGTCGTGTTGTGTGGCGAAAATCGAg               <  7:308848/148‑1 (MQ=255)
                                                                                                                                   cagctGGTATCATATTCTGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGCGAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATGGCGGTAGTCGTGTTGTGTGGCGAAAATCGAgc              >  8:222264/5‑149 (MQ=255)
                                                                                                                                     gctGGTATCATATTCTGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGCGAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATGGCGGTAGTCGTGTTGTGTGGCGaaaa                    <  1:283452/139‑1 (MQ=255)
                                                                                                                                     gctGGTATCATATTCTGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGCGAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATGGCGGTAGTCGTGTTGTGTGGCGaaaa                    >  2:283452/3‑141 (MQ=255)
                                                                                                                                           aTCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGCGAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATGGCGGTAGTCGTGTTGTGTGGCGAAAATCGAGcgccgcgct      >  4:105949/1‑146 (MQ=255)
                                                                                                                                           aTCAGAATATGATACCAGCTCATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGCGAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATGGCGGTAGTCGTGTTGTGTGGCGAAAATCGACCGCCGCCac      >  6:333508/1‑148 (MQ=255)
                                                                                                                                              aGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGCGAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATGGCGGTAGTCGTGTTGTGTGGCGAAAATCGAGcgccgcgctggc   >  3:338135/1‑143 (MQ=255)
                                                                                                                                               gAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGCGAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATGGCGGTAGTCGTGTTGTGTGGCGAAAATCGAGCGCCGCCAGcgcg  >  3:277290/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GCACCCGCGATCCACGGCCTTCGATCGCCTTGGCGCCTATCGATATCAGAAATGACTCCAAACGCGGCCACTCCAGCGACTTTGGAACCGGCTTTTTGAAAATGGCTTCAAAGGTCGATCTTTGATGATTATTCATGGTATCAGAATATGATACCAGCTGATTGTTGGTATGTAGGATGGGCGGAGCCCTATCTGTAGGAGTGGTCTTGCGAGCTAGGGGTGCTGGGCGTCGACGTTAGCCGTGCATGGCGGTAGTCGTGTTGTGTGGCGAAAATCGAGCGCCGCCAGCGCG  >  NC_002947/1684987‑1685278

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑