Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,958,799 G→C 10.3% intergenic (+65/+77) PP_2589 → / ← PP_2590 aldehyde dehydrogenase family protein/ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,958,7990GC10.3% 117.5 / 5.2 47intergenic (+65/+77)PP_2589/PP_2590aldehyde dehydrogenase family protein/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base G (23/19);  new base C (2/3);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.48e-01

GGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGT  >  NC_002947/2958654‑2958942
                                                                                                                                                 |                                                                                                                                               
gggTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCt                                                                                                                                              <  7:252986/149‑1 (MQ=255)
  gTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGa                                                                                                                                            <  1:21022/149‑1 (MQ=255)
        ccGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTtg                                                                                                                                      >  3:111522/1‑149 (MQ=255)
          gcgATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTtgtg                                                                                                                                    >  7:238662/1‑149 (MQ=255)
          gagATCTGTCACTGGATTCGTTCGACAAGTATTCGCATCTCAAGACCTCCTTGTACGAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTGCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTtgtg                                                                                                                                    <  3:232431/147‑1 (MQ=255)
                tGTCACTGCATTCGTTCGACAAGTATTAGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGAGTCGTCTGCTTCGCGGGCATGCCAGCTCCCACAGGTTCAGAACAAGTCTGAGGGTTGTGCTGGAc                                                                                                                              >  5:203098/1‑149 (MQ=255)
                  tCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCGCGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAAcc                                                                                                                             <  5:54009/148‑1 (MQ=25)
                  tCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTg                                                                                                                                 <  7:101093/144‑1 (MQ=255)
                  tCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGGCGTCTGCTTCGCGGGCGCGCCCGCTCCCACCGGTCCAGCACAACCCTCAGACTTGTGCTGAACCt                                                                                                                            <  6:64405/149‑1 (MQ=25)
                    aCTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCACGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCtgt                                                                                                                          <  8:190625/149‑1 (MQ=25)
                     cTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgtg                                                                                                                         >  8:211887/1‑149 (MQ=255)
                      tGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTgg                                                                                                                        >  1:250658/1‑149 (MQ=255)
                        cATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCGGAGGGTTGTGCTGGACCTGTgggc                                                                                                                      >  3:43067/1‑148 (MQ=255)
                          ttcgttcgACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGc                                                                                                                    <  3:359441/149‑1 (MQ=255)
                             gttcgACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACAGCGTCGTCTGCTTCGCGGTCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGAGCTGGACCTGTGGGAGCggg                                                                                                                 >  1:320589/1‑149 (MQ=255)
                             gttcgACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACAGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGAGCTGGACCTGTGGGAGCggg                                                                                                                 >  1:236475/1‑149 (MQ=255)
                              ttcgACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGgc                                                                                                                <  8:191137/149‑1 (MQ=255)
                                       tATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACGTGTGGGAGCGGGCGCGCCcgc                                                                                                        >  8:311683/1‑148 (MQ=255)
                                         ttCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGa                                                                                                      >  6:2876/1‑148 (MQ=255)
                                            gCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAagag                                                                                                  >  1:81109/1‑149 (MQ=255)
                                               aCTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGcc                                                                                               >  3:277686/1‑149 (MQ=255)
                                                 tCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAACTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCg                                                                                              >  6:112595/1‑148 (MQ=255)
                                                                 cctgTTGCGCGGCTGATGTTTTGATCTTTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACACGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCgcggcg                                                                             <  4:277686/146‑1 (MQ=255)
                                                                  cAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGcggcgg                                                                            >  4:125197/1‑149 (MQ=255)
                                                                          cgGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGt                                                                    <  2:250658/149‑1 (MQ=255)
                                                                            gttgttGTTTTTATCATTGACCGGGTCGTCTGCTTCGCGGTAATGCCGGCTCCCACAGGTTCAGCACATGTGTGTGGTTTGTGCTGGACCTTTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTa                                                                  <  2:81109/144‑1 (MQ=37)
                                                                                aTGTATTGTTCTTTGTCCGCGTCGTCTGCTTCGCGGGCTTGCCGGCTCCCACAGTTTGAGCAGAAGTCTGAGGGTGGGGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGaa                                                               <  2:236475/148‑1 (MQ=37)
                                                                                         tCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTa                                                     <  3:125197/149‑1 (MQ=255)
                                                                                            ttGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGATCTTGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGcc                                                  >  3:149653/1‑149 (MQ=255)
                                                                                                ccGCGTGGGCTGCTCGGCGGGCATCCCCGCTCCCACCGGTTCAGCACACGGCTGCGGGTTGTGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGaaa                                              <  4:149653/149‑1 (MQ=37)
                                                                                                  gcgTCGTCTGGTTCGCGGGCATGCCCGCTCCCACTGGTTCCGCACACGGCTGAGGGTTGTGCTGGACTTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGc                                            <  5:22865/149‑1 (MQ=39)
                                                                                                   cgtcgtCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCGGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCgcggc                                                                              >  7:135222/1‑114 (MQ=255)
                                                                                                       ggcTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCCGCCCACGTCTGAGGGTGGTGCGGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATa                                       <  8:238662/147‑1 (MQ=38)
                                                                                                         cTGCTTCGCGGGCATACCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGAt                                        <  6:212845/146‑1 (MQ=255)
                                                                                                         cTGCTTCGCGGGCATACCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGAt                                        >  5:212845/1‑146 (MQ=255)
                                                                                                                     ctggCCCGCTCCCACCGGTGCAGCACACGTCGGAGGGGGGTGCGGGACCTGGGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGccc                         <  5:112595/146‑1 (MQ=255)
                                                                                                                        gCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCCCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGccccc                       >  4:272651/1‑148 (MQ=255)
                                                                                                                        gCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCGGTGGGAGCGGGCGCGCCCGCGAAGCGTCCGTATCTGGCTCCTCGGCGGGACCCGGTTAGAAGTCGACGCTAGCCGAAAGCTGCTATGTCCGGGGCGCCCCCa                      >  4:267823/1‑149 (MQ=255)
                                                                                                                              tCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAAccc                  >  3:104505/1‑147 (MQ=21)
                                                                                                                                ccACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAgg              >  2:265110/1‑149 (MQ=255)
                                                                                                                                ccACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAgg              >  4:181525/1‑149 (MQ=255)
                                                                                                                                ccACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAgg              >  5:366310/1‑149 (MQ=255)
                                                                                                                                    cggTTCAGCACATGTCTGAGGGTTGGGCGGGACCTGGGGGCGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCt          <  7:311683/148‑1 (MQ=255)
                                                                                                                                           gCACATGTCTGTGGGTTGTGGTGGTGCTGTGGGTGGGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCGGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAg   <  2:320589/149‑1 (MQ=255)
                                                                                                                                           gCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAggg             >  2:82916/1‑139 (MQ=255)
                                                                                                                                           gCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAggg             <  1:82916/139‑1 (MQ=255)
                                                                                                                                            cacaAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGt  >  4:114305/1‑149 (MQ=255)
                                                                                                                                                 gTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTgg        <  6:336779/138‑1 (MQ=255)
                                                                                                                                                 gTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTgg        >  5:336779/1‑138 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGT  >  NC_002947/2958654‑2958942

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