Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,914,197 T→C 8.0% intergenic (+33/+120) spuI → / ← argA glutamylpolyamine synthetase/amino acid N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,914,1970TC 8.0% 144.8 / 2.7 50intergenic (+33/+120)spuI/argAglutamylpolyamine synthetase/amino acid N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (22/24);  new base C (2/2);  total (24/26)
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.82e-01

ATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTT  >  NC_002947/5914054‑5914340
                                                                                                                                               |                                                                                                                                               
aTCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTga                                                                                                                                           <  1:798831/150‑1 (MQ=255)
     tCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCg                                                                                                                                      <  1:1609344/150‑1 (MQ=255)
      ccGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCgg                                                                                                                                     <  1:660806/150‑1 (MQ=255)
               aTCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCct                                                                                                                            >  1:1522269/1‑150 (MQ=255)
                           gtcgCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCAc                                                                                                                <  1:441802/150‑1 (MQ=255)
                                 tGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGc                                                                                                          >  2:802130/1‑150 (MQ=255)
                                    aaGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTcc                                                                                                       >  1:8949/1‑150 (MQ=255)
                                          aGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAgg                                                                                                 <  2:701777/150‑1 (MQ=255)
                                          aGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAgg                                                                                                 <  2:1522269/150‑1 (MQ=255)
                                            cGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGt                                                                                               <  1:802130/150‑1 (MQ=255)
                                                        gTTCGAGCACTCGATCTCCGACCTCGATTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTga                                                                                                                                           <  2:666422/94‑1 (MQ=255)
                                                        gTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTga                                                                                                                                           >  1:666422/1‑94 (MQ=255)
                                                            gAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCCGCCGCCGGGCACTTTCGCAggatcga                                                                               >  2:300201/1‑145 (MQ=255)
                                                                  tCGATCTCCCACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCtgt                                                                         <  1:213773/150‑1 (MQ=255)
                                                                    gATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTcgc                                                                                                                      >  2:573043/1‑103 (MQ=255)
                                                                    gATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTcgc                                                                                                                      <  1:573043/103‑1 (MQ=255)
                                                                    gATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGa                                                                       <  1:300201/150‑1 (MQ=255)
                                                                    gATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGa                                                                       >  2:1301884/1‑150 (MQ=255)
                                                                       ctcCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACtt                                                                    >  1:13007/1‑150 (MQ=255)
                                                                        tcCGACGGCGAGTACGACTGGTACCTGGATAGCGGGTAAAGGAAAACGCCGCCCCTCAAGGGCGGCGTTGTTGTTTGAGACCGGTGTGGCCCTCTTGGCGGGCAGGCCCGCTGCCGCAGGGGAATTTGGAGGGCTGAAGCGTGTGACTTa                                                                   >  2:1256453/1‑150 (MQ=255)
                                                                        tcCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTa                                                                   <  1:1301884/150‑1 (MQ=255)
                                                                         ccGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTAc                                                                  <  2:13007/150‑1 (MQ=255)
                                                                             ccTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCtgt                                                                         <  1:1119899/139‑1 (MQ=255)
                                                                             ccTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCACGCAGGGTAATTTCGAAGGCTGAAGCCtgt                                                                         >  2:1119899/1‑139 (MQ=255)
                                                                                    tACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCtgt                                                                         >  2:461498/1‑132 (MQ=255)
                                                                                    tACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCtgt                                                                         <  1:461498/132‑1 (MQ=255)
                                                                                          tGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTg                                                 >  1:1077977/1‑150 (MQ=255)
                                                                                                cTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCcgcg                                           <  2:1077977/150‑1 (MQ=255)
                                                                                                     tACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCggg                                                     <  2:262793/135‑1 (MQ=255)
                                                                                                     tACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCggg                                                     >  1:262793/1‑135 (MQ=255)
                                                                                                         tctcaaaCAAAAACGCCGCCCTTGAGGGGCGGCGTTTTCGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCtgtg                                                             <  1:640915/119‑1 (MQ=37)
                                                                                                         tctcaaaCAAAAACGCCGCCCTTGAGGGGCGGCGTTTTCGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCtgtg                                                             >  2:640915/5‑123 (MQ=37)
                                                                                                          tgtAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCtgt                                                                         <  2:1503967/110‑1 (MQ=255)
                                                                                                          tgtAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCtgt                                                                         >  1:1503967/1‑110 (MQ=255)
                                                                                                          tgtAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCtgt                                                              >  1:367081/1‑121 (MQ=255)
                                                                                                          tgtAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCtgt                                                              <  2:367081/121‑1 (MQ=255)
                                                                                                              aaCGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGt                                                  <  1:1550189/129‑1 (MQ=255)
                                                                                                              aaCGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGt                                                  >  2:1550189/1‑129 (MQ=255)
                                                                                                               acaaaaaCGCCGCCCTTGAGGGGCGGCGTTTTCGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGcc                                                                                                             >  2:1584957/4‑69 (MQ=17)
                                                                                                               acaaaaaCGCCGCCCTTGAGGGGCGGCGTTTTCGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGcc                                                                                                             <  1:1584957/66‑1 (MQ=17)
                                                                                                                cGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAgg                                     <  2:1557815/140‑1 (MQ=255)
                                                                                                                cGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCGCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAgg                                     >  1:1557815/1‑140 (MQ=255)
                                                                                                                    aaCGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAgg                                     >  1:408693/1‑136 (MQ=255)
                                                                                                                    aaCGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAgg                                     <  2:408693/136‑1 (MQ=255)
                                                                                                                     aCGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGaa                                                                                      >  2:87274/1‑86 (MQ=255)
                                                                                                                     aCGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGaa                                                                                      <  1:87274/86‑1 (MQ=255)
                                                                                                                         cgccCCTCAAGGGCGGCGTTTTTGTGAGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAggg                  >  1:631556/1‑150 (MQ=255)
                                                                                                                             ccTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCcgc                                            >  2:630396/1‑120 (MQ=255)
                                                                                                                             ccTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCcgc                                            <  1:630396/120‑1 (MQ=255)
                                                                                                                                         cgTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGAAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCtt  >  1:205678/1‑150 (MQ=255)
                                                                                                                                           tttttGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAgcgc                                   <  1:1642938/115‑1 (MQ=255)
                                                                                                                                           tttttGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGggc                                   >  2:1642938/1‑115 (MQ=255)
                                                                                                                                               |                                                                                                                                               
ATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTT  >  NC_002947/5914054‑5914340

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