Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 322,135 G→A 25.0% intergenic (+420/+38) 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,1350GA25.0% 59.4 / 21.3 40intergenic (+420/+38)aguA/PP_0267agmatine deiminase/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base G (15/15);  new base A (5/5);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.49e-01

ACCATGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTC  >  NC_002947/321990‑322278
                                                                                                                                                 |                                                                                                                                               
acCATGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCAt                                                                                                                                             >  1:215129/1‑150 (MQ=255)
  cATGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGc                                                                                                                                           >  1:247569/1‑150 (MQ=255)
           aTCCGGGTCTTGTCGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                  <  2:410239/150‑1 (MQ=255)
           aTCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                  >  2:282037/1‑150 (MQ=255)
           aTCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAaga                                                                                                                                  <  2:247569/150‑1 (MQ=255)
                    aGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTAc                                                                                                                         <  1:151282/150‑1 (MQ=255)
                            cgaccCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAagag                                                                                                                                 >  1:467357/1‑134 (MQ=255)
                            cgaccCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAagag                                                                                                                                 <  2:467357/134‑1 (MQ=255)
                                ccACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTaaa                                                                                                             <  1:269092/150‑1 (MQ=255)
                                       gTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGaa                                                                                                      >  2:24401/1‑150 (MQ=255)
                                          tGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTTATCGATTTGCCCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACt                                                                                                   >  1:318270/1‑150 (MQ=255)
                                          tGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAg                                                                                                                       <  2:38938/130‑1 (MQ=255)
                                          tGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAg                                                                                                                       >  1:38938/1‑130 (MQ=255)
                                                            cccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGaa                                                                                                                                    <  1:156453/99‑1 (MQ=255)
                                                            cccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGaa                                                                                                                                    >  2:156453/1‑99 (MQ=255)
                                                             ccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCCTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACcaca                                                                                >  2:159465/1‑150 (MQ=255)
                                                                      gATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTAcgc                                                                       <  2:335005/150‑1 (MQ=255)
                                                                      gATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCCTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTAcgc                                                                       <  1:159465/150‑1 (MQ=255)
                                                                          aGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAaga                                                                                                                                  <  2:523848/87‑1 (MQ=255)
                                                                          aGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAaga                                                                                                                                  >  1:523848/1‑87 (MQ=255)
                                                                                    agaTGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGGTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGGTAACCGTCACCACAAGGCTACGCGGATCGCCCACCTg                                                         >  1:92602/1‑150 (MQ=255)
                                                                                          aGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAg                                                                                                                             <  2:539487/76‑1 (MQ=255)
                                                                                          aGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAg                                                                                                                             >  1:539487/1‑76 (MQ=255)
                                                                                          aGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTg                                                   <  1:27285/150‑1 (MQ=255)
                                                                                                        gtgtGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTgg                                        >  1:452198/1‑147 (MQ=255)
                                                                                                        gtgtGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTgg                                        <  2:452198/147‑1 (MQ=255)
                                                                                                            ggCACGGATCGATGGCACCTGTACCTGCCTCGTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCGACCCGTCACCACAAGGCTACGCGGATCGCACACCTGGATCTGCGCCGCACTGGTCGCCGa                                 <  2:217762/150‑1 (MQ=255)
                                                                                                               aCTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCt                              <  1:282037/150‑1 (MQ=255)
                                                                                                                       aTTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAg                                                                                                        <  2:193536/68‑1 (MQ=255)
                                                                                                                       aTTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAg                                                                                                        >  1:193536/1‑68 (MQ=255)
                                                                                                                        ttGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGt                     <  2:318270/150‑1 (MQ=255)
                                                                                                                           cACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCtt                  <  1:360219/150‑1 (MQ=255)
                                                                                                                                gTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCACGTACAGGCAACTTAAAATCAGAAACTGTCGGTACCCGTCACCACAAGGCTACGCGAAACGACCACATAGATCAACCCCACACTGGTCGCAGCGCCGTAGTAGATCTTCTCAg             >  2:117532/1‑150 (MQ=255)
                                                                                                                                        ctctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCtt                  <  2:216109/137‑1 (MQ=255)
                                                                                                                                        ctctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCtt                  >  1:216109/1‑137 (MQ=255)
                                                                                                                                         tctTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGc    >  1:482364/1‑150 (MQ=255)
                                                                                                                                         tctTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGc    >  1:122843/1‑150 (MQ=255)
                                                                                                                                         tctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGTAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTgta                           <  1:46370/127‑1 (MQ=255)
                                                                                                                                         tctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGTAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTgta                           >  2:46370/1‑127 (MQ=255)
                                                                                                                                           ttCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTc  >  2:225481/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
ACCATGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTC  >  NC_002947/321990‑322278

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: