Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,102,760 A→G 8.2% M1M (ATG→GTG) † phhB → pterin‑4‑alpha‑carbinolamine dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,102,7600AG 8.2% 160.6 / 3.2 57M1M (ATG→GTG) †phhBpterin‑4‑alpha‑carbinolamine dehydratase
Reads supporting (aligned to +/- strand):  ref base A (28/24);  new base G (3/2);  total (31/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 = 1.00e+00

ACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGAAGTACGTGACGGCCACATGGAGCTGGAGCGCGTGTT  >  NC_002947/5102611‑5102905
                                                                                                                                                     |                                                                                                                                                 
aCATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCa                                                                                                                                                   <  1:1078588/150‑1 (MQ=255)
  aTGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCGTg                                                                                                                                                 >  1:91774/1‑150 (MQ=255)
           cTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTg                                                                                                                                        >  1:520821/1‑150 (MQ=255)
             gTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGACAACACCTCATGAATGCCTTGaa                                                                                                                                      <  1:583885/150‑1 (MQ=255)
                  aCCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAg                                                                                                                                 >  1:1042539/1‑150 (MQ=255)
                         ccACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTg                                                                                                                          <  2:250685/150‑1 (MQ=255)
                           aCGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCg                                                                                                                        >  1:742255/1‑150 (MQ=255)
                             gAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGaa                                                                                                                      >  2:1055685/1‑150 (MQ=255)
                                   atcatGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCCCCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGc                                                                                                                >  1:412744/1‑150 (MQ=255)
                                   atcatGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCt                                                                                                                                          >  1:371227/1‑124 (MQ=255)
                                   atcatGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCt                                                                                                                                          <  2:371227/124‑1 (MQ=255)
                                       tGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCcgcg                                                                                                            >  1:611843/1‑150 (MQ=255)
                                           cATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGaa                                                                                                                                               <  1:268900/111‑1 (MQ=255)
                                           cATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGaa                                                                                                                                               >  2:268900/1‑111 (MQ=255)
                                                   aTAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGa                                                                                                >  1:505776/1‑150 (MQ=255)
                                                     aaaGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGaaa                                                                                              <  2:520821/150‑1 (MQ=255)
                                                         ccATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGGCTGCCGCGCCGATGCGCCGAAAGTCa                                                                                          >  2:4725/1‑150 (MQ=255)
                                                           aTGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGaaa                                                                                              >  2:1058724/1‑144 (MQ=255)
                                                           aTGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGaaa                                                                                              <  1:1058724/144‑1 (MQ=255)
                                                                     gTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGc                                                                              <  2:909085/150‑1 (MQ=255)
                                                                        tgcatgcaCCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTgg                                                                           <  1:735047/150‑1 (MQ=255)
                                                                           atgcaCCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCg                                                                        >  2:684710/1‑150 (MQ=255)
                                                                               aCCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCt                                                                    >  1:386800/1‑150 (MQ=255)
                                                                               aCCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCt                                                                    >  1:564526/1‑150 (MQ=255)
                                                                               aCCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCt                                                                    >  1:564519/1‑150 (MQ=255)
                                                                                   aaGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATc                                                                >  2:1099130/1‑150 (MQ=255)
                                                                                           aCCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTTAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAAt                                                        >  2:513360/1‑150 (MQ=255)
                                                                                               aaGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTTAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCg                                                    <  1:513360/150‑1 (MQ=255)
                                                                                               aaGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCg                                                    <  1:36036/150‑1 (MQ=255)
                                                                                                 ggTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGa                                                  <  2:611843/150‑1 (MQ=255)
                                                                                                     gCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGTCGAGCTGATCCGCGAAATCCCGGACTgg                                              <  1:1055685/150‑1 (MQ=255)
                                                                                                     gCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTgg                                              <  2:412744/150‑1 (MQ=255)
                                                                                                      cTGCCTGAGCGACCTTGCTGATAACAACTCGAACCGGAAAACACCTCGTGAATGCCTTGAACCAAGCCCATTGCGAAGCCTTCCGCGCCGATGCGCCGAAAGTCa                                                                                          >  2:533612/1‑105 (MQ=255)
                                                                                                      cTGCCTGAGCGACCTTGCTGATAACAACTCGAACCGGAAAACACCTCGTGAATGCCTTGAACCAAGCCCATTGCGAAGCCTTCCGCGCCGATGCGCCGAAAGTCa                                                                                          >  2:533602/1‑105 (MQ=255)
                                                                                                      cTGCCTGAGCGACCTTGCTGATAACAACTCGAACCGGAAAACACCTCGTGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCa                                                                                          <  1:533602/105‑1 (MQ=255)
                                                                                                      cTGCCTGAGCGACCTTGCTGATAACAACTCGAACCGGAAAACACCTCGTGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCa                                                                                          <  1:533612/105‑1 (MQ=255)
                                                                                                       tGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGaa                                            >  1:982994/1‑150 (MQ=255)
                                                                                                              gCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGaa                                     <  2:1042539/150‑1 (MQ=255)
                                                                                                               cGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGATCGGTCCCGGTTTCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCNAAg                                    <  2:705021/150‑1 (MQ=255)
                                                                                                                  ccTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGAAGTAc                                 >  1:783244/1‑150 (MQ=255)
                                                                                                                   cTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGAAGTACg                                <  1:509506/150‑1 (MQ=255)
                                                                                                                    ttGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGaa                                     <  1:344428/144‑1 (MQ=255)
                                                                                                                    ttGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGaa                                     >  2:344428/1‑144 (MQ=255)
                                                                                                                    ttGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGAAGTACGt                               <  2:742255/150‑1 (MQ=255)
                                                                                                                    ttGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGAAGTACGt                               >  2:437067/1‑150 (MQ=255)
                                                                                                                    ttGCCGATAACAACTCGAACCGGAAAACAACCCATAAAAGCTTATCATGGATATCGTTGCGAAGCATGCCTGTCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGCACATCGAAGTActt                               >  1:193368/1‑148 (MQ=255)
                                                                                                                        cGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGAAGTACGTGACg                           >  2:895042/1‑150 (MQ=255)
                                                                                                                         gATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGAAGTACGTGACgg                          >  1:104655/1‑150 (MQ=255)
                                                                                                                               aacTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTgg                                                                           <  2:752923/95‑1 (MQ=255)
                                                                                                                               aacTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTgg                                                                           >  1:752923/1‑95 (MQ=255)
                                                                                                                                acTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGAAGTAc                                 <  1:640337/136‑1 (MQ=255)
                                                                                                                                acTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGAAGTAc                                 >  2:640337/1‑136 (MQ=255)
                                                                                                                                  tCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGAAGTACGTGACGGCCACAtgga                 <  2:982994/150‑1 (MQ=255)
                                                                                                                                  tCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGAAGTACGTGACGGCCACAtgga                 >  2:293085/1‑150 (MQ=255)
                                                                                                                                            aaaCACCTCAGGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGAAGTACGTGACGGCCACATGGAGCTGGAgcgc       <  1:461763/150‑1 (MQ=255)
                                                                                                                                             aaCACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGAAGTACGTGACGGCCACATGGAGCTGGAGcgcg      >  2:760492/1‑150 (MQ=255)
                                                                                                                                                acCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGAAGTACGTGACGGCCACATGGAGCTGGAGTGCgtgt   <  1:895042/150‑1 (MQ=255)
                                                                                                                                                 ccTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAACCACGGACTGGAACATCGAAGTACGTGACGGCCACATGGAGCTGGAGCGCGTGtt  <  1:202383/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
ACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGAAGTACGTGACGGCCACATGGAGCTGGAGCGCGTGTT  >  NC_002947/5102611‑5102905

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: