Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 321,743 G→A 11.2% intergenic (+28/+430) aguA → / ← PP_0267 agmatine deiminase/ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947321,7430GA11.2% 110.6 / 5.4 45intergenic (+28/+430)aguA/PP_0267agmatine deiminase/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base G (21/18);  new base A (3/2);  total (24/21)
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.49e-01

GCAAGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCTCACGGGCCAG  >  NC_002947/321602‑321892
                                                                                                                                             |                                                                                                                                                     
gCAAGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGc                                                                                                                                               <  2:427636/150‑1 (MQ=255)
ccctggATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCAGGCGCCGGTCAAGCGCTGAGCCAGTAATGAATGCGCGGCCCGTCATGTGACGGGc                                                                                                                                               >  1:26054/5‑150 (MQ=255)
  aaGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCct                                                                                                                                             >  1:732509/1‑150 (MQ=255)
             ccGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCttt                                                                                                                                         >  2:424684/1‑143 (MQ=255)
             ccGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCttt                                                                                                                                         <  1:424684/143‑1 (MQ=255)
              cGACCATGAGTTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCCAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGcgcg                                                                                                                                 >  1:199509/1‑150 (MQ=255)
                    tGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGAc                                                                                                                           >  1:110907/1‑150 (MQ=255)
                        gtggtgATGATCCCGGGCCGCAAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCAGTAATGAATGAGAGGCCCGTCATGTGACGGGCATCTTTTTTGCGCGGGCGACGcac                                                                                                                       >  2:176665/1‑148 (MQ=255)
                           gtgATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGACCGTAACATATGTGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCAGGCGGCGCGCGGa                                                                                                                    >  2:551864/1‑150 (MQ=255)
                           gtgATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACACAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGa                                                                                                                    >  2:666105/1‑150 (MQ=255)
                            tgatgaTCCCGGGCCGCGTGTATCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCACCCGGCGCCGGTCAAGCGCTGAGCCCGTAATTAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGaa                                                                                                                   <  1:413971/150‑1 (MQ=255)
                                gatCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCttt                                                                                                                                         >  1:322097/1‑124 (MQ=255)
                                gatCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCttt                                                                                                                                         <  2:322097/124‑1 (MQ=255)
                                gatCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCg                                                                                                               >  2:438975/1‑150 (MQ=255)
                                gatCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCACATGACgg                                                                                                                                                 <  2:408112/116‑1 (MQ=255)
                                gatCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCACATGACgg                                                                                                                                                 >  1:408112/1‑116 (MQ=255)
                                  tCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAg                                                                                                             >  1:89806/1‑150 (MQ=255)
                                   cccGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGc                                                                                                            <  2:89806/150‑1 (MQ=255)
                                     cGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTg                                                                                                          >  1:104034/1‑150 (MQ=255)
                                      gggCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTgg                                                                                                         <  1:454845/150‑1 (MQ=255)
                                             gAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGt                                                                                                  >  2:37688/1‑150 (MQ=255)
                                                  gctCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGt                                                                                             >  1:457984/1‑150 (MQ=255)
                                                   ctctTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGAGACGGTCTAGTGCTGATCCCGTAAAAAAAGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTa                                                                                            <  1:355934/150‑1 (MQ=255)
                                                       tGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGtt                                                                                       <  2:385237/150‑1 (MQ=255)
                                                                 cAACATTCACTTCCTGACCCAACCGCATCCGGCGCCGGTCAATCCCTTAGCCCGTAATGAATTAGAGGCCCGTCATGTTACGGGCCTCTTTTTTACGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCa                                                                              <  1:176665/150‑1 (MQ=255)
                                                                 cAACATTCACTGCCTGACCCAACAGAAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAAAAAACGAGAGACCAGTCATGTGAAGGCACACTTTTTTGCGTGGTCGACGCCTTGTAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCa                                                                              <  2:273966/150‑1 (MQ=255)
                                                                         aCTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCaa                                                                      >  1:326782/1‑150 (MQ=255)
                                                                             ccTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATg                                                                  >  2:675512/1‑150 (MQ=255)
                                                                                   ccAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGAACTGCCTCCGCCGATCTCCTTGAGGCCTGCCGCGGGGGGTCGCGGGCGACGACCGGAAGGAGAGCTGTCATACGTTTTGTATTGTTTTTGAAGCAGATAGCAAGATGGCACgg                                                            >  1:578496/1‑150 (MQ=255)
                                                                                            gcCGGCGCCGGTCGCGAGCTGCGAACGTAATCTATTAGACGCCCTTCTTTTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCcgc                                                   <  2:457984/150‑1 (MQ=255)
                                                                                                     ggTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGAt                                          <  2:266320/150‑1 (MQ=255)
                                                                                                          aGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTa                                     >  2:380838/1‑150 (MQ=255)
                                                                                                           gcgcTGAGCACGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCGGGAGGCGGAGCTGGCATGAGTTTTGTATTGTTTTCCCCGCAGAGAGCAACATAGCACGGTCAGGCCGCGCCTCTGATCTGTaa                                    <  2:26054/150‑1 (MQ=255)
                                                                                                              cTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACaa                                 <  1:534853/150‑1 (MQ=255)
                                                                                                                     cGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTAc                          <  1:337145/150‑1 (MQ=255)
                                                                                                                           gaatgaGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAAt                    <  2:473154/150‑1 (MQ=255)
                                                                                                                                 gagGCCCGTCACATGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGtttt                                                                                     <  1:533820/79‑1 (MQ=255)
                                                                                                                                 gagGCCCGTCACATGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGtttt                                                                                     >  2:533820/1‑79 (MQ=38)
                                                                                                                                  agGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGc             >  2:399010/1‑150 (MQ=255)
                                                                                                                                  agGCCCGTCACATGGCGGGCCTTGTTTTGGGGCGTGCGACTGCCGGAAGCCGAGCTGGCGTACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATCGCACGGCCAGCCCTCGTCTCTGATCTTTAACAATACCTACGTAAATTAGCCGc             >  2:163958/1‑150 (MQ=255)
                                                                                                                                      ccGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTCCGTAAATTAGCCGCTCAc         >  1:118583/1‑150 (MQ=255)
                                                                                                                                      ccGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCTCAc         <  1:37688/150‑1 (MQ=255)
                                                                                                                                       cGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTaa                                    >  2:273139/1‑122 (MQ=255)
                                                                                                                                       cGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTaa                                    <  1:273139/122‑1 (MQ=255)
                                                                                                                                        gTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCTCACgg       <  1:399010/150‑1 (MQ=255)
                                                                                                                                             gtgACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCTCACGGGCCAg  <  2:326782/150‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                     
GCAAGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCTCACGGGCCAG  >  NC_002947/321602‑321892

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