Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,350,531 C→T 8.9% intergenic (+20/+62) PP_4704 → / ← PP_4705 hypothetical protein/glyoxalase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,350,5310CT 8.9% 155.3 / 5.4 56intergenic (+20/+62)PP_4704/PP_4705hypothetical protein/glyoxalase
Reads supporting (aligned to +/- strand):  ref base C (31/20);  new base T (3/2);  total (34/22)
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.46e-01

TCAAGGACTATCTGGTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGCTGCGTGC  >  NC_002947/5350387‑5350679
                                                                                                                                                |                                                                                                                                                      
tCAAGGACTATCTGGTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCcac                                                                                                                                                   >  1:1210988/1‑150 (MQ=255)
       cTATCTGGTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTg                                                                                                                                          <  1:621545/150‑1 (MQ=255)
              gTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTa                                                                                                                                   >  2:1323889/1‑150 (MQ=255)
               ggtgacggcgcacgcagcttccacctgctggacccggccgggaacaaggttcaaccgctgtatcacccggcggtgtcggggcagcggttggtttgaagcccgtgctggcctcttcgcgggcacGCCCGCTCCCACAGGCCTGCGCACTAc                                                                                                                                  >  2:731697‑M1/124‑150 (MQ=255)
                 gACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGc                                                                                                                                               <  1:587981/137‑1 (MQ=255)
                 gACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGc                                                                                                                                               >  2:587981/1‑137 (MQ=255)
                    aaGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCtgtg                                                                                                                             >  2:861421/1‑150 (MQ=255)
                     aGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTgg                                                                                                                            <  1:125165/150‑1 (MQ=255)
                       cgcacgcagcttccacctgctggacccggccgggaacaaggttcaaccgctgtatcacccggcggtgtcggggcagcggttggtttgaagcccgtgctggcctcttcgcgggcacGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGa                                                                                                                          >  2:1249711‑M1/116‑150 (MQ=255)
                         cAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGc                                                                                                                        <  1:659021/150‑1 (MQ=255)
                          aaGCTTGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCg                                                                                                                       <  1:340934/150‑1 (MQ=255)
                               ggAAGTCGACTAGGTGTACTTCATTGATACTAGCCTGGGGCAGCCTCAAAGAGATGACCGGGTAGTGTCACCGATGTGCCTGGATACAGAGTCATCCCCTCTGGCGGGACAGCCCCGCGAGGCCTGCGCGCTATCTGTGGGAGGGGGCGt                                                                                                                  >  2:681911/1‑150 (MQ=255)
                               gcttccacctgctggacccggccgggaacaaggttcaaccgctgtatcacccggcggtgtcggggcagcggttggtttgaagcccgtgctggcctcttcgcgggcacGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGt                                                                                                                  >  2:782840‑M1/108‑150 (MQ=255)
                                             tgtACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGc                                                                                                    <  1:1323889/150‑1 (MQ=255)
                                               tACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCa                                                                                                  >  1:762406/1‑150 (MQ=255)
                                                ccggccgggaacaaggttcaaccgctgtatcacccggcggtgtcggggcagcggttggtttgaagcccgtgctggcctcttcgcgggcacGCCCGCCCCCACAGGCCTGCGCACTACCtgtg                                                                                                                             >  1:363118‑M1/91‑122 (MQ=255)
                                                ccggccgggaacaaggttcaaccgctgtatcacccggcggtgtcggggcagcggttggtttgaagcccgtgctggcctcttcgcgggcacGCCCGCCCCCACAGGCCTGCGCACTACCtgtg                                                                                                                             <  2:363118‑M1/32‑1 (MQ=255)
                                                   gccgggaacaaggttcaaccgctgtatcacccggcggtgtcggggcagcggttggtttgaagcccgtgctggcctcttcgcgggcacGCCCGCCCCCACAGGCCTGCGCACTACCTGTgg                                                                                                                            <  1:451920‑M1/33‑1 (MQ=255)
                                                   gccgggaacaaggttcaaccgctgtatcacccggcggtgtcggggcagcggttggtttgaagcccgtgctggcctcttcgcgggcacGCCCGCCCCCACAGGCCTGCGCACTACCTGTgg                                                                                                                            >  2:451920‑M1/88‑120 (MQ=255)
                                                     aTCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACgg                                                                                            >  2:1396755/1‑150 (MQ=255)
                                                      tCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGc                                                                                                                                               >  1:180035/1‑100 (MQ=255)
                                                      tCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGc                                                                                                                                               <  2:180035/100‑1 (MQ=255)
                                                      tCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACggg                                                                                           <  2:606593/150‑1 (MQ=255)
                                                      tCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACggg                                                                                           >  1:1344005/1‑150 (MQ=255)
                                                                         gCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCCCTGcc                                                                        >  1:141852/1‑150 (MQ=255)
                                                                           cTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGcccc                                                                      >  2:1254468/1‑150 (MQ=255)
                                                                            tCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTTGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCg                                                                     >  1:336708/1‑150 (MQ=255)
                                                                             cAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGa                                                                    <  1:1254468/150‑1 (MQ=255)
                                                                               acccggcggtgtcggggcagcggttgttttgaagcccgtgctggcctcttcgcgggcacGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCa                                                                                                  >  2:903558‑M1/60‑118 (MQ=255)
                                                                                    aTGGCCGGGTTGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACAccgcc                                                             >  1:1265637/1‑150 (MQ=255)
                                                                                        ccGGGTGGTGACACCGATGCACCTGGATGCCGAGTCAGCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGt                                                         <  2:336708/150‑1 (MQ=255)
                                                                                                   cacCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCAAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGt                                              <  2:346893/150‑1 (MQ=255)
                                                                                                    acCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGtt                                             >  1:33114/1‑150 (MQ=255)
                                                                                                        aTGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAAc                                         <  1:8298/150‑1 (MQ=255)
                                                                                                        aTGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAAc                                         <  2:1344005/150‑1 (MQ=255)
                                                                                                              cTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGtt                                   >  2:618929/1‑150 (MQ=255)
                                                                                                              cTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGtt                                   <  2:762406/150‑1 (MQ=255)
                                                                                                              cTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGtt                                   >  1:992257/1‑150 (MQ=255)
                                                                                                                      cGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCccggccg                           >  2:124273/1‑150 (MQ=255)
                                                                                                                      cGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCccggccg                           <  2:1210988/150‑1 (MQ=255)
                                                                                                                        tggcctcttcgcgggcacGCCCGCTCCCACAGGTAGTGCGCAGGCCTGTGGGGGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGt                                                         >  1:903558‑M1/19‑118 (MQ=255)
                                                                                                                             tCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCt                                                                           >  1:1149929/1‑95 (MQ=255)
                                                                                                                             tCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCt                                                                           <  2:1149929/95‑1 (MQ=255)
                                                                                                                               cccTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCcagca                  >  2:1215650/1‑150 (MQ=255)
                                                                                                                               cccTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCcagca                  >  1:217633/1‑150 (MQ=255)
                                                                                                                                ccTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCcagcag                 >  1:125291/1‑150 (MQ=255)
                                                                                                                                 ctctGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAgg                >  1:58632/1‑150 (MQ=255)
                                                                                                                                   cgggcacGCCCGCTCCCACAGGTAGTGCGCAGGCCTGTGGGGGCGGGCCTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTg              <  2:731720‑M1/143‑1 (MQ=255)
                                                                                                                                    gggcacGCCCGCTCCCACAGGTAGTGCGCAGGCCTGTGGGGGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGt                        <  1:705306‑M1/133‑1 (MQ=255)
                                                                                                                                    gggcacGCCCGCTCCCACAGGTAGTGCGCAGGCCTGTGGGGGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGt                        >  2:705306‑M1/7‑139 (MQ=255)
                                                                                                                                       cAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCg                                                <  1:1031191/112‑1 (MQ=255)
                                                                                                                                       cAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCg                                                >  2:1031191/1‑112 (MQ=255)
                                                                                                                                        aGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGc         >  1:940848/1‑150 (MQ=255)
                                                                                                                                        aGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGc         >  2:1211884/1‑150 (MQ=255)
                                                                                                                                          GCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGCtg       <  1:851214‑M1/150‑1 (MQ=255)
                                                                                                                                               CCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGCtgcgtgc  >  1:1396755‑M1/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
TCAAGGACTATCTGGTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGCTGCGTGC  >  NC_002947/5350387‑5350679

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