Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,671,702 G→T 9.0% intergenic (+121/‑49) PP_4978 → / → PP_4979 hypothetical protein/substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,671,7020GT 9.0% 121.9 / 2.1 43intergenic (+121/‑49)PP_4978/PP_4979hypothetical protein/substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (22/17);  new base T (2/2);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

CTGCGATCCACCCAAGCCTGATGTGCCCGCAAAACTCATTGTCTTGCGCGGTCTTGTGGGAGCGGGCGTGCCCGCGAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGCCCGTTGGCTCTGGGCGAACGCTTGCGCCTGGTAACCGACGACTGGGCAC  >  NC_002947/5671561‑5671850
                                                                                                                                             |                                                                                                                                                    
ctgcGATCCACCCAAGCCTGATGTGCCCGCAAAACTCATTGTCTTGCGCGGTCTTGTGGGAGCGGGCGTGCCCGCGAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATa                                                                                                                                                   <  4:97154/145‑1 (MQ=255)
   cGATCCACCCACGCCTGATGTGCCCGCAAAACTCATTGTCTTGCGCGGTCTTGTGGGAGCGGGCGTGCCCGCTAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAatgcatg                                                                                                                                            <  5:203517/149‑1 (MQ=255)
          cccAAGCCTGATGTGCCCGCAAAACTCATTGTCTTGCGCGGTCTTGTGGGAGCGGGCGTGCCCGCGAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTAc                                                                                                                                     <  6:309706/149‑1 (MQ=255)
                           cGCAAAACTCATTGTCTTGCGCGGTCTTGTGGGAGCGGGCGTGCCCGCGAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAgcg                                                                                                                    >  1:241279/1‑149 (MQ=255)
                           cGCAAAACTCATTGTCTTGCGCGGTCTTGTGGGAGCGGGCGTGCCCGCGAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAgcg                                                                                                                    >  6:161610/1‑149 (MQ=255)
                                   tCATTGTCTTGCGCGGTCTTGTGGGAGCGGGCGTGCCCGCGAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACCCCGGCCGAgcgagatagg                                                                                                             >  7:56644/1‑141 (MQ=255)
                                   tCATTGTCTTGCGCGGTCTTGTGGGAGCGGGCGTGCCCGCGAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAgc                                                                                                                     <  8:56644/140‑1 (MQ=255)
                                   tCATTGTCTTGCGCGGTCTTGTGGGAGCGGGCGTGCCCGCGAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTAGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAAtc                                                                                                             >  7:4956/1‑148 (MQ=255)
                                      ttGTCTTGCGCGGTCTTGTGGGAGCGGGCGTGCCCGCGAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCg                                                                                                         <  5:161610/149‑1 (MQ=255)
                                        gTCTTGCTCGGTCTTGTGGGAGCGGGCGTGCCCGCGAAGAGAGCAACCCGATCCTTGCCCCCGGCTTCGCCGGTGCTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTAGACCGGCCGAGCGCGTAATCTCCgg                                                                                                        <  5:248370/148‑1 (MQ=38)
                                           ttGCGCGGTCTTGTGGGAGCGGGCGTGCCCGCGAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTAGCGGGCACGCCCGCTCCCACATCTATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGcc                                                                                                      >  6:260080/1‑147 (MQ=255)
                                                       GTGGGAGCGGGCGTGCCCGCGAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGc                                                                                        <  4:89144‑M1/149‑1 (MQ=255)
                                                               GGGCGTGCCCGCGAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGt                                                                                <  2:241558‑M1/149‑1 (MQ=255)
                                                                ggcgtgcccgcgaacaccggcgaagccggtgccaacgatcgggttgctctcttcgcgggcacgcccgctcccacaTCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGtt                                                                               >  5:126973‑M2/76‑149 (MQ=255)
                                                                      cccgcgaacaccggcgaagccggtgccaacgatcgggttgctctcttcgcgggcacgcccgctcccacaTCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGAt                                                                         >  6:191719‑M2/70‑149 (MQ=255)
                                                                        CGCGAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATcc                                                                       >  5:167880‑M1/1‑149 (MQ=255)
                                                                         GCGAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCct                                                                      >  2:241279‑M1/1‑149 (MQ=255)
                                                                             acaccggcgaagccggtgccaacgatcgggttgctctcttcgcgggcacgcccgctcccacaTCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCctgctg                                                                  <  6:126973‑M2/87‑1 (MQ=255)
                                                                               agagCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTgctgct                                                                >  6:49062/1‑149 (MQ=255)
                                                                                      cccGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCTATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGc                                                         >  1:60172/1‑149 (MQ=255)
                                                                                      cccGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGc                                                         >  5:241563/1‑149 (MQ=255)
                                                                                                  cACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGCCCGTTg                                             <  6:241563/149‑1 (MQ=255)
                                                                                                  cACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGCCCGTTg                                             >  6:235922/1‑149 (MQ=255)
                                                                                                   aCCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGCCCGTTgg                                            >  4:194815/1‑149 (MQ=255)
                                                                                                         ttCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGCCCGTTg                                             <  1:43269/142‑1 (MQ=255)
                                                                                                         ttCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGCCCGTTg                                             >  2:43269/1‑142 (MQ=255)
                                                                                                         ttCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGCCCGTTGGCTCTgg                                      >  6:293988/1‑149 (MQ=255)
                                                                                                           cGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCTATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGGCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGCCCGTTGGCTCTgg                                      <  5:260080/147‑1 (MQ=255)
                                                                                                           cGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGCCCGTTGGCTCTgg                                      <  2:224122/147‑1 (MQ=255)
                                                                                                            gCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGCCCGTTGGCTCTGGGc                                    >  8:226605/1‑148 (MQ=255)
                                                                                                            gCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCAACATATCGACCGTCTGTTCTGGATCCTGCTCCTCGTTTGCCTGAGCCCGTAGGCTCTGGGc                                    >  6:201620/1‑148 (MQ=255)
                                                                                                               ggTGTTCGCGGGCACGCCCGCTCCCACATCTATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGCCCGTTGGCTCTGGGCGAAc                                <  2:60172/149‑1 (MQ=255)
                                                                                                               ggTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGCCCGTTGGCTCTGGGCGAAc                                <  3:89144/149‑1 (MQ=255)
                                                                                                                  gtTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTGCTGGATCCTGCTGCTCGTTTGCCTGAGCCCGTTGGCTCTGGGCGAACGCt                             >  4:92418/1‑149 (MQ=255)
                                                                                                                     cgcgggcacgcccgctcccacaTCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGCCCGTTGGCTCTGGGCGAACGCtt                            <  5:191719‑M2/125‑1 (MQ=255)
                                                                                                                     cgcgggcacgcccgctcccacaTCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGCCCGTTGGCTCTGGGCGAACGCTTgc                          >  7:60564‑M2/23‑149 (MQ=255)
                                                                                                                      gcgggcacgcccgctcccacaTCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGCCCGTTGGCTCTGGGCGAACGCTTgcg                         >  6:194010‑M2/22‑149 (MQ=255)
                                                                                                                         ggcacgcccgctcccacaTAGATAATGCATGCCTTTACCGTACACCGGCCGATCGCGAAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGCATCCTGCTGCTCGTTTGACTCAGCCCGATGGCTCTGGGCGAACGCATGCCCCt                      >  7:86347‑M2/19‑149 (MQ=255)
                                                                                                                           cacgcccgctcccacaTCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGCCCGTTGGCTCTGGGCGAACGCTTGCGCCTgg                    <  5:49062‑M2/133‑1 (MQ=255)
                                                                                                                            acgcccgctcccacaTCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGCCCGTTGGCTCTGGGCGAACGCTTGCGCCTGGt                   <  5:111291‑M2/134‑1 (MQ=255)
                                                                                                                                      ccacaTCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGCCCGTTGGCTCTGGGCGAACGCTTGCGCCTGGTAACCgacgac         >  5:303650‑M2/6‑149 (MQ=255)
                                                                                                                                          aTCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGCCCGTTGGCTCTGGGCGAACGCTTGCGCCTGGTAACCGACGACTggg     >  7:69998‑M2/2‑149 (MQ=255)
                                                                                                                                          aTCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGCCCGTTGGCTCTGGGCGAACGCTTGCGCCTGGTAACCGACGACTggg     <  6:303650‑M2/148‑1 (MQ=255)
                                                                                                                                             gATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGCCCGTTGGCTCTGGGCGAACGCTTGCGCCTGGTAACCGACGACTGGGCa   >  3:96860/1‑148 (MQ=255)
                                                                                                                                             gATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGCCCGTTGGCTCTGGGCGAACGCTTGCGCCTGGTAACCGACGACTGGGCAc  <  5:293988/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
CTGCGATCCACCCAAGCCTGATGTGCCCGCAAAACTCATTGTCTTGCGCGGTCTTGTGGGAGCGGGCGTGCCCGCGAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGCCCGTTGGCTCTGGGCGAACGCTTGCGCCTGGTAACCGACGACTGGGCAC  >  NC_002947/5671561‑5671850

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: ‑