Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,337,301 G→A 56.2% G235R (GGG→AGG)  PP_4696 → nitrogen regulation protein NR(I)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,337,3010GA56.2% 7.6 / 54.1 48G235R (GGG→AGG) PP_4696nitrogen regulation protein NR(I)
Reads supporting (aligned to +/- strand):  ref base G (10/11);  new base A (16/11);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 5.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.45e-01

AAAGAGCTGGTGGCCCGCGCCCTGCACAACCTGTCTCGCCGGGCCAAGGCACCGATGATCTCGGTGAACTGCGCGGCCATCCCCGAAACGCTGATCGAATCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCGGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTCGATGAAATTGGCGAACTGCCGCTCGAAGCCCAGGCTCGCCTGCTGCGTGTGCTGCAGGAAGGCGAGATACGCCGGGTAGGTTCGGTGCAGTCGCAAAAGGT  >  NC_002947/5337154‑5337449
                                                                                                                                                   |                                                                                                                                                    
aaaGAGCTGGTGGCCCGCGCCCTGCACAACCTGTCTCGCCGGGCCAAGGCACCGATGATCTCGGTGAACTGCGCGGCCATCCCCGAAACGCTGATCGAATCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCTGGCGCCggt                                                                                                                                                    <  1:80799/150‑2 (MQ=255)
 aaGAGCTGGTGGACCGAGAACTGGAGAAATTGTTTGGTCGGGACAATGCACCGATGAGCTCAAGGAACGCCGACGCCATCCCCTAAACGCTGAGCGAATCCGAACTGTTCGGCCAGGACAGAGGCGCGTTCGCCGGACCCAGCGCCAGgc                                                                                                                                                   >  2:11548/1‑150 (MQ=255)
        ggtggCCCGCGCCCTGCACAACCTGTCTCGCCGGGCCAAGGCACCGATGATCTCGGTGAACTGCGCGGCCATCCCCGAAACGCTGATCGAATCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCGGGCGCGCCgg                                                                                                                                            >  1:730902/1‑150 (MQ=255)
              ccGCGCCCTGCACAACCTGTCTCGCCGGGCCAAGGCACCGATGATCTCGGTGAACTGCGCGGCCATCCCCGAAACGCTGATCGAATCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCAGGCGCGCCGGCTtgg                                                                                                                                       >  1:1109148/1‑149 (MQ=255)
              ccGCGCCCTGCACAACCTGTCTCGCCGGGCCAAGGCACCGATGATCTCGGTGAACTGCGCGGCCATCCCCGAAACGCTGATCGAATCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCAGGCGCGCCGGCTtgg                                                                                                                                       <  2:1109148/149‑1 (MQ=255)
               cgcgCCCGGCGCAACCTGGCTCGCCGGGCGAAGGGACCGATGACCTAGGTGAACTGCGCGGCCATCCCCGAAACGCTGATCGAATCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCAGGCGCGCCGGCTtggtg                                                                                                                                     <  2:400705/150‑1 (MQ=38)
                             ccTGTCTCGCCGGGCCAAGGCACCGATGATCTCGGTGAACTGCGCGGCCATCCCCGAAACGCTGATCGAATCCGAACTGTGCTGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCGGGCGCGCCGGCTGGGTGGAGGCCGCCGATgg                                                                                                                       >  1:801492/1‑150 (MQ=255)
                              cTGTCTCGCCGGGCCAAGGCACCGATGATCTCGGTGAACTGCGCGGCCATCCCCGAAACGCTGATCGAATCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCAGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGc                                                                                                                      >  2:1070157/1‑150 (MQ=255)
                                 tctcGCCGGGCCAAGGCACCGATGATCTCGGTGAACTGCGCGGCCATCCCCGAAACGCTGATCGAATCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCAGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGt                                                                                                                   >  2:6649/1‑150 (MQ=255)
                                   tcGCCGGGCCAAGGCACCGATGATCTCGGTGAACTGCGCGGCCATCCCCGAAACGCTGATCGAATCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCAGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTAc                                                                                                                 >  2:792557/1‑150 (MQ=255)
                                      ccGGGCCAAGGCACCGATGATCTCGGTGAACTGCGCGGCCATCCCCGAAACGCTGATCGAATCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCGGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCt                                                                                                              <  1:923973/150‑1 (MQ=255)
                                      ccGGGCCAAGGCACCGATGATCTCGGTGAACTGCGCGGCCATCCCCGAAACGCTGATCGAATCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCGGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCt                                                                                                              <  1:509098/150‑1 (MQ=255)
                                       cGGGCCAAGGCACCGATGATCTCGGTGAACTGCGCGGCCATCCCCGAAACGCTGATCGAATCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCAGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTg                                                                                                             >  1:501784/1‑150 (MQ=255)
                                           ccAAGGCACCGATGATCTCGGTGAACTGCGCGGCCATCCCCGAAACGCTGATCGAATCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCGGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTcc                                                                                                         >  2:567216/1‑150 (MQ=255)
                                            cAAGGCACCGATGATCCCGGTGAACTGCGCGGCCATCCCCGAAACGCTGATCGAATCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCAGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCt                                                                                                        >  2:549596/1‑150 (MQ=255)
                                             aaGGCACCGATGATCTCGGTGAACTGCGCGGCCATCCCCGAAACGCTGATCGAATCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCAGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTc                                                                                                       >  1:909139/1‑150 (MQ=255)
                                                gCACCGATGATCTCGGTGAACTGCGCGGCCATCCCCGAAACGCTGATCGAATCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCAGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTCGAt                                                                                                    >  2:961395/1‑150 (MQ=255)
                                                  aCCGATGATCTCGGTGAACTGCGCGGCCATCCCCGAAACGCTGATCGAATCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCGGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTCGATGa                                                                                                  >  1:387890/1‑150 (MQ=255)
                                                       tgatCTCGGTGAACTGCGCGGCCATCCCCGAAACGCTGATCGAATCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCAGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTCGATGAAATTg                                                                                             <  1:311788/150‑1 (MQ=255)
                                                            tcGGTGAACTGCGCGGCCATCCCCGAAACGCTGATCGAATCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCGGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTCGATGAAATTGGCGaa                                                                                        >  1:144972/1‑150 (MQ=255)
                                                                     tGCGCGGCCATCCCCGAAACGCTGATCGAATCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCGGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTCGATGAAATTGGCGAACTGCCGCTc                                                                               <  2:387890/150‑1 (MQ=255)
                                                                        gcgGCCATCCCCGAAACGCTGATCGAATCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCGGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTCGATGAAATTGGCGAACTGCCGCTCGaa                                                                            <  1:939602/150‑1 (MQ=255)
                                                                            ccATCCCCGAAACGCTGATCGAATCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCAGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTCGATGAAATTGGCGAACTGCCGCTCGAAGccc                                                                        <  1:792557/150‑1 (MQ=255)
                                                                                 cccGAAACGCTGATCGAATCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCAGGCGCGCCGGCTTGGTGGAGGCCGCCGAt                                                                                                                         >  1:795270/1‑96 (MQ=255)
                                                                                 cccGAAACGCTGATCGAATCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCAGGCGCGCCGGCTTGGTGGAGGCCGCCGAt                                                                                                                         <  2:795270/96‑1 (MQ=255)
                                                                                  ccGAAACGCTGATCGAATCCGAACTGTTCGGCCATCAAAAAGGCGCGTTCACCGGCGCCAGCGCCAGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTCGATGAAATTGGCGAACTGCCGCTCGAAGCCCAGGCTc                                                                  >  1:640855/1‑150 (MQ=255)
                                                                                  ccGAAACGCTGATCGAATCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCAGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTCGATGAAATTGGCGAACTGCCGCTCGAAGCCCAGGCTc                                                                  <  1:1070157/150‑1 (MQ=255)
                                                                                  ccGAAACGCTGATCGAATCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCAGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTCGATGAAATTGGCGAACTGCCGCTCGAAGCCCAGGCTc                                                                  >  1:198104/1‑150 (MQ=255)
                                                                                       aCGCTGATCGAATCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCAGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTCGATGAAATTGGCGAACTGCCGCTCGAAGCCCAGGCTCGCctg                                                             <  1:961395/150‑1 (MQ=255)
                                                                                            gATCGAATCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCAGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTCGATGAAATTGGCGAACTGCCGCTCGAAGCCCAGGCTCGCCTGCTGCg                                                        >  2:232176/1‑150 (MQ=255)
                                                                                                gAATCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCGGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTCGATGAAATTGGCGAACTGCCGCTCGAAGCCCAGGCTCGCCTGCTGCGtgtg                                                    >  2:1098198/1‑150 (MQ=255)
                                                                                                   tCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCAGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCg                                                                                                                     >  1:1015255/1‑82 (MQ=255)
                                                                                                   tCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCAGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCg                                                                                                                     <  2:1015255/82‑1 (MQ=255)
                                                                                                     cGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCGGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTCGATGAAATTGGCGAACTGCCGCTCGAAGCCCAGGCTCGCCTGCTGCGTGTGCTGCa                                               <  1:1029428/150‑1 (MQ=255)
                                                                                                              cGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCAGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTCGATGAAATTGGCGAACTGCCGCTCGAAGCCCAGGCTCGCCTGCTGCGTGTGCTGCAGGAAGGCga                                      <  2:198104/150‑1 (MQ=255)
                                                                                                              cGGCCACGAAAAAGGCGCGTTCACCGGAACCCGCATGAGGCGCGTCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTCGATGAAATTGGCGAACTGCCGCTCGAAGCCCAGGCTCGCCTGCTGCGTGTGCTGCAGGAAGGCga                                      <  2:640855/150‑1 (MQ=255)
                                                                                                                gCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCGGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTCGATGAAATTGGCGAACTGCCGCTCGAAGCCCAGGCTCGCCTGCTGCGTGTGCTGCAGGAAGGCgaga                                    <  1:567216/150‑1 (MQ=255)
                                                                                                                  cACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCAGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTCGATGAAATTGGCGAACTGCCGCTCGAAGACCAGGCTCGCCTGCTGCGTGTGCTGCAGGAAGGCGAGATa                                  >  2:687219/1‑150 (MQ=255)
                                                                                                                      aaaaaGGCGCGTTCACCGGCGCCAGCGCCAGGCGCGCCGGCTTGGTGNAGGCCGCCGATGGCGGTACGCTGTTCCTCGATGAAATTGGCGAACTGCCGCTCGAAGCCCAGGCTCGCCTGCTGCGTGTGCTGCAGGAAGGCGAGATACGcc                              >  1:54516/1‑150 (MQ=255)
                                                                                                                       aaaaGGCGCGTTCACCGGCGCCAGCGCCGGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTCGATGAAATTGGCGaa                                                                                        <  2:1131440/91‑1 (MQ=255)
                                                                                                                       aaaaGGCGCGTTCACCGGCGCCAGCGCCGGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTCGATGAAATTGGCGaa                                                                                        >  1:1131440/1‑91 (MQ=255)
                                                                                                                            gcgcGTTCACCGGCGCCAGCGCCGGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTCGATGAAATTGGCGAACTGCCGCTCGAAGCCCAGGCTCGCCTGCTGCGTGTGCTGCAGGAAGGCGAGATACGCCGGGTAg                        >  2:1063338/1‑150 (MQ=255)
                                                                                                                             cgcgTTCACCGGCGCCAGCGCCGGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTCGATGAAATTGGCGAACTGCCGCTCGAAGCCCAGGCTCGCCTGCTGCGTGTGCTGCAGGAAGGCGAGATACGCCGGGTatg                       <  2:801492/150‑3 (MQ=255)
                                                                                                                             cgccTTCACCGGCGCCAGCGCCAGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTCGATGAAATTGGCGAACTGCCGCTCGAAGCCCAGGCTCGCCTGCTGCGTGTGCTGCAGGAAGGCGAGATACGCCGGGTAgg                       <  1:549596/150‑1 (MQ=255)
                                                                                                                                gTTCACCGGCGCCAGCGCCGGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTCGATGAAATTGGCGAACTGCCGCTCGAAGCCCAGGCTCGCCTGCTGCGTGTGCTGCAGGAAGGCGAGATACGCCGGGTAGGTTc                    >  1:82859/1‑150 (MQ=255)
                                                                                                                                   cACCGGCGCCAGCGCCGGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTCGATGAAATTGGCGAACTGCCGCTCGAAGCCCAGGCTCGCCTGCTGCGTGTGCTGCAGGAAGGCGAGATACGCCGGGTAGGTTCGGt                 <  2:1097000/150‑1 (MQ=255)
                                                                                                                                      cGGCGCCAGCGCCGGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTCGATGAAATTGGCGAACTGCCGCTCGAAGCCCAGGCTCGCCTGCTGCGTGTGCTGCAGGAAGGCGAGATACGCCGGGTAGGTTCGGTGCa              >  2:1100338/1‑150 (MQ=255)
                                                                                                                                                gcCGGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTCGATGAAATTGGCGAACTGCCGCTCGAAGCCCAGGCTCGCCTGCTGCGTGTGCTGCAGGAAGGCGAGATACGCCGGGTAGGTTCGGTGCAGTCGCAAAAg    <  1:1063338/150‑1 (MQ=255)
                                                                                                                                                  cGGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTCGATGAAATTGGCGAACTGCCGCTCGAAGCCCAGGCTCGCCTGCTGCGTGTGCTGCAGGAAGGCGAGATACGCCGGGTAGGTTCGGTGCAGTCGCAAAAGGt  <  2:730902/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AAAGAGCTGGTGGCCCGCGCCCTGCACAACCTGTCTCGCCGGGCCAAGGCACCGATGATCTCGGTGAACTGCGCGGCCATCCCCGAAACGCTGATCGAATCCGAACTGTTCGGCCACGAAAAAGGCGCGTTCACCGGCGCCAGCGCCGGGCGCGCCGGCTTGGTGGAGGCCGCCGATGGCGGTACGCTGTTCCTCGATGAAATTGGCGAACTGCCGCTCGAAGCCCAGGCTCGCCTGCTGCGTGTGCTGCAGGAAGGCGAGATACGCCGGGTAGGTTCGGTGCAGTCGCAAAAGGT  >  NC_002947/5337154‑5337449

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