Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 322,140 T→C 7.3% intergenic (+427/+31) 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_002947322,1400TC 7.3% 154.6 / 2.4 55intergenic (+427/+31)aguA/PP_0267agmatine deiminase/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base T (21/30);  new base C (2/2);  total (23/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

AAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAATGCCGCCCG  >  NC_002947/321995‑322287
                                                                                                                                                 |                                                                                                                                                     
aaGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCcgc                                                                                                                                                    >  6:173181/1‑149 (MQ=255)
     tCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGaa                                                                                                                                                 >  6:75318/1‑147 (MQ=255)
     tCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAa‑ga                                                                                                                                              <  5:173181/149‑1 (MQ=255)
     tCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAa‑ga                                                                                                                                              <  6:250257/149‑1 (MQ=255)
       aGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATAGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAgg                                                                                                                                            >  2:260862/1‑149 (MQ=255)
        ggggCAGGACGACCGACCCACGACGGCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGc                                                                                                                                           <  2:234273/149‑1 (MQ=255)
          gTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑g                                                                                                                                        <  7:220360/149‑1 (MQ=255)
              ggACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTAc                                                                                                                                    <  5:75318/149‑1 (MQ=255)
                 cgaccgacCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAgg                                                                                                                                 <  2:68065/149‑1 (MQ=255)
                  gaccgaccCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAa‑gag                                                                                                                                             >  2:336711/1‑137 (MQ=255)
                  gaccgaccCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAa‑gag                                                                                                                                             <  1:336711/137‑1 (MQ=255)
                         ccACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTaa                                                                                                                         <  1:260862/149‑1 (MQ=255)
                         ccACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTa                                                                                                                          >  4:84820/1‑148 (MQ=255)
                          cACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTaaa                                                                                                                        <  3:90860/149‑1 (MQ=255)
                          cACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATACTCGCGAA‑GAGGCA‑GGTACAGGCAACTTaa                                                                                                                         >  5:315700/1‑148 (MQ=255)
                          cACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTCCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTaaa                                                                                                                        <  7:147459/149‑1 (MQ=255)
                           acgacgTCTGCTTGCCGGGCGATACCCCCAAGGGGTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGAACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTaaaa                                                                                                                       <  5:213192/149‑1 (MQ=255)
                                 tCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGa                                                                                                                  >  2:25575/1‑148 (MQ=255)
                                     cTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGaaa                                                                                                                <  7:54050/146‑1 (MQ=255)
                                        gCCGGGCGATACGCGCATGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGTTGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTGTTCGCGGGCATGCTCGCGAA‑GAGGCATGGTACAGGCAACTTAAAATCAGAAACTGt                                                                                                            <  1:25575/148‑1 (MQ=255)
                                        gCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTa                                                                                                           <  3:93779/148‑1 (MQ=255)
                                            ggCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTAc                                                                                                      <  3:84820/149‑1 (MQ=255)
                                                 tACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGt                                                                                                  >  5:136461/1‑148 (MQ=255)
                                                    ccccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTcac                                                                                               >  5:337131/1‑148 (MQ=255)
                                                    ccccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTcac                                                                                               >  4:289292/1‑148 (MQ=255)
                                                      ccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTcacca                                                                                             >  8:270334/1‑148 (MQ=255)
                                                        aaGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATACTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACaa                                                                                          <  6:315700/149‑1 (MQ=255)
                                                               gATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTAcg                                                                                   >  4:137586/1‑149 (MQ=255)
                                                                           agagATGCAGGCGCTTATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCAc                                                                       <  2:277880/149‑1 (MQ=255)
                                                                           agagATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCAc                                                                       <  3:289292/149‑1 (MQ=255)
                                                                           agagATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCAc                                                                       <  2:281708/149‑1 (MQ=255)
                                                                             agaTGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCt                                                                     <  7:270334/149‑1 (MQ=255)
                                                                              gaTGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTg                                                                    >  1:149725/1‑149 (MQ=255)
                                                                                 gCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACaa                                                                                          <  7:290869/124‑1 (MQ=255)
                                                                                 gCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACaa                                                                                          >  8:290869/1‑124 (MQ=255)
                                                                                         tGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTaaaa                                                                                                                       <  1:154053/87‑1 (MQ=255)
                                                                                         tGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTaaaa                                                                                                                       >  2:154053/1‑87 (MQ=255)
                                                                                         tGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGcgccg                                                         >  6:345839/1‑149 (MQ=255)
                                                                                              gCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTg                                                    <  2:149725/149‑1 (MQ=255)
                                                                                                       cACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCAc                                                      <  4:296121/138‑1 (MQ=255)
                                                                                                       cACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCAc                                                      >  3:296121/1‑138 (MQ=255)
                                                                                                        aCTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCCCCTGGATCTGCGCCGCACTGGTCGCCGAGc                                          >  3:249185/1‑149 (MQ=255)
                                                                                                          tGATCGATTGCACCTGTACCTGCCTCTTCGCGTGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGcgcc                                                          <  6:256197/131‑1 (MQ=255)
                                                                                                          tGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGcgcc                                                          >  5:256197/1‑131 (MQ=255)
                                                                                                            aTCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTg                                        <  5:345839/147‑1 (MQ=255)
                                                                                                                    cACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTgtagta                                   <  8:68675/144‑1 (MQ=255)
                                                                                                                    cACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGtagtagagat                              >  7:68675/1‑145 (MQ=255)
                                                                                                                    cACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTc                               <  6:136461/148‑1 (MQ=255)
                                                                                                                          tACCTGCCTCTTCGCGGGCATGCTCGCGAAGGAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGtgttgt                           >  1:30927/1‑147 (MQ=255)
                                                                                                                          tACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAg                        <  3:137586/149‑1 (MQ=255)
                                                                                                                                  tctTCGCGAGCATGCCCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCGATGTTg                >  7:58124/1‑149 (MQ=255)
                                                                                                                                    ttCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTg                <  2:30927/147‑1 (MQ=255)
                                                                                                                                        cgGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAAt          <  3:65842/149‑1 (MQ=255)
                                                                                                                                                cTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAATGccgcccg  <  7:6284/149‑1 (MQ=255)
                                                                                                                                                cTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAATGccgcccg  >  8:6284/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
AAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAA‑GAGGCA‑GGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAATGCCGCCCG  >  NC_002947/321995‑322287

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