Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,787,722 G→A 8.2% A156T (GCC→ACC)  pyrH → uridylate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,787,7220GA 8.2% 132.2 / 2.8 49A156T (GCC→ACC) pyrHuridylate kinase
Reads supporting (aligned to +/- strand):  ref base G (21/24);  new base A (2/2);  total (23/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.37e-01

GGTCATGTCGGCCATTTCCATGGTGGGTGTCACCGATCATTACGATCGCCGCAAAGCTATTCGTCACCTCAACTCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGATGGCGTGTACACTGCCGATCCATTCAAGGACCCGCATGCCGAGAAATTCGATCACCTGACCTACGATGAGGTTCTGGA  >  NC_002947/1787574‑1787846
                                                                                                                                                    |                                                                                                                            
ggTCATGTCGGCCATTTCCATGATAGCTGTCACCGATCATTACGATCGCCGCAAAGCTATTCGTCACCTCAACTCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCAGCCTGCCTgcgcgc                                                                                                                             >  1:150731/1‑150 (MQ=255)
           ccATTTCCATGGTGGGTGTCACCGATCATTACGATCGCCGCAAAGCTATTCGTCACCTCAACTCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCg                                                                                                                  >  1:16626/1‑150 (MQ=255)
             aTTTCCATGGTGGGTGTCACCGATCATTACGATCGCCGCAAAGCTATTCGTCACCTCAACTCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGAt                                                                                                                <  1:773295/150‑1 (MQ=255)
             aTTTCCATGGTGGGTGTCACCGATCATTACGATCGCCGCAAAGCTATTCGTCACCTCAACTCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGAt                                                                                                                <  1:242620/150‑1 (MQ=255)
             aTTTCCATGGTGGGTGTCACCGACCATTACGACCGCCAAAAAGCAATAGGAAACACTAAATGAGGGGATGATGTTTTTTTTGCCGACGGTACCGGCAACTCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCAAAATCGAt                                                                                                                >  2:126594/1‑150 (MQ=255)
                      gtgggtgTCACCGATCATTACGATCGCCGCAAAGCTATTCGTCACCTCAACTCCGGGGATGTGGTAATTTTCTCCCGCGTTACCGGCAACCCGTTCTTCACCACAGACTCCGCGGCCGGCCGGCGCGCCATCGAAATCGATGCTGACgtg                                                                                                       >  1:518782/1‑150 (MQ=255)
                        gggtgTCACCGATCATTACGATCGCCGCAAAGCTATTCGTCACCGCAACTCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACgtggt                                                                                                     >  2:224380/1‑150 (MQ=255)
                             tCACCGATCATTACGAGCGCCGCAAAGCTATTCGTCACCTCAACTCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGa                                                                                                >  1:214077/1‑150 (MQ=255)
                              cACCGATCATTACGATCGCCGCAAAGCTATTCGTCACCTCAACTCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGaa                                                                                               <  1:273428/150‑1 (MQ=255)
                               aCCGATCATTACGATCGCCGCAAAGCTATTCGTCACCTCAACTCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAg                                                                                              >  1:517247/1‑150 (MQ=255)
                               aCCGATCATTACGATCGCCGCAAAGCTATTCGTCACCTCAACTCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAg                                                                                              >  1:1087050/1‑150 (MQ=255)
                                ccGATCATTACGATCGCCGCAAAGCTATTCGTCACCTCAACTCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCCCCGACCCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAgg                                                                                             <  1:302102/150‑1 (MQ=255)
                                ccGATCATTACGATCGCCGCAAAGCTATTCGTCACCTCAACTCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAgg                                                                                             >  1:471080/1‑150 (MQ=255)
                                   aTCATTACGATCGCCGCAAAGCTATTCGTCACCTCAACTCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGa                                                                                          <  2:517247/150‑1 (MQ=255)
                                          cGATCGCCGCAAAGCTATTCGTCACCTCAACTCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGt                                                                                   <  2:16626/150‑1 (MQ=255)
                                             tCGCCGCAAAGCTATTCGTCACCTCAACTCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGa                                                                                <  1:449639/150‑1 (MQ=255)
                                              cgccgcAAAGCTATTCGTCACCTCAACTCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGAt                                                                               <  2:646569/150‑1 (MQ=255)
                                                ccgcAAAGCTATTCGTCACCTCAACTCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGATgg                                                                             >  2:503387/1‑150 (MQ=255)
                                                    aaaGCTATTCGTCACCTCACCTCCGGGGATGTGGTAATTCTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGATGGCgtg                                                                         <  1:1029990/150‑1 (MQ=255)
                                                          aTTCGTCACCTCAACTCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGATGGCGTGTACACt                                                                   <  1:503387/150‑1 (MQ=255)
                                                           ttCGTCACCTCAACTCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCACCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGATGGCGTGTACACTg                                                                  >  2:865709/1‑150 (MQ=255)
                                                           ttCGTCACCTCAACTCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCACCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGATGGCGTGTACACTg                                                                  >  2:865829/1‑150 (MQ=255)
                                                            tCGTCACCTCAACTCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGATGGCGTGTACACTGc                                                                 <  1:791300/150‑1 (MQ=255)
                                                             cGTCACCTCAACTCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCACCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAATGTCGATGGCGTGTACACTGcc                                                                <  1:865709/150‑1 (MQ=255)
                                                             cGTCACCTCAACTCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCACCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGATGGCGTGTACACTGcc                                                                <  1:865829/150‑1 (MQ=255)
                                                                 cttccgaCTCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGATGGCGTGTACACTGCCGATc                                                            >  1:38698/7‑150 (MQ=255)
                                                                   cTCAACTCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGATGGCGTGTACACTGCCGa                                                              >  1:164563/1‑146 (MQ=255)
                                                                   cTCAACTCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGATGGCGTGTACACTGCCGa                                                              <  2:164563/146‑1 (MQ=255)
                                                                   cTCAACTCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGATGGCGTGTACACTGCCGATCCa                                                          >  1:627066/1‑150 (MQ=255)
                                                                        cTCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGATGGCGTGTACACTGCCGATCCATTCaa                                                     <  1:573751/150‑1 (MQ=255)
                                                                        cTCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGATGGCGTGTACACTGCCGATCCATTCaa                                                     <  2:87434/150‑1 (MQ=255)
                                                                         tCCNGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGGCGTGGTGCTGAAGGCGACCAAGGTCGATGGCGTGTACACTGCCGATCCATTCAAg                                                    >  2:702036/1‑150 (MQ=255)
                                                                         tCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGATGGCGTGTACACTGCCGATCCATTCAAg                                                    >  2:421848/1‑150 (MQ=255)
                                                                         tCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGATGGCGTGTACACTGCCGATCCATTCAAg                                                    >  1:350728/1‑150 (MQ=255)
                                                                           cGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGACATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGATGGCGTGTACACTGCCGATCCATTCAAGGa                                                  >  2:1113600/1‑150 (MQ=255)
                                                                            ggggATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGATGGCGTGTACACTGCCGATCCATTCAAGGAc                                                 <  1:421848/150‑1 (MQ=255)
                                                                              ggATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGATGGCGTGTACACTGCCGATCCATTCAAGGAccc                                               <  2:471080/150‑1 (MQ=255)
                                                                                              tcCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGATGGCGTGTACACTGCCGATCCATTCAAGGACCCGCATGCCGAGAAATTc                               <  2:350728/150‑1 (MQ=255)
                                                                                                         cGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGATGGCGTGTACACTGCCGATCCATTCAAGGACCCGCATGCCGAGAAATTCGATCACCTGAc                    <  1:1113600/150‑1 (MQ=255)
                                                                                                          ggCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGATGGCGTGTACACTGCCGATCCATTCAAGGACCCGCATGCCGAGAAATTCGATCACCTGAcc                   <  2:518782/150‑1 (MQ=255)
                                                                                                           gCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGATGGCGTGTACACTGCCGATCCATTCAAGGACCCGCATGCCGAGAAATTCGATCACCTGACCt                  >  1:929398/1‑150 (MQ=255)
                                                                                                               cccGTTCTTCACCACCGCCTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGATGGCGTGTACACTGCCGATCCATTCTAGGACCCGCATGCCGAGAAATTCGATCACCTGACCTACGa              <  2:293010/150‑1 (MQ=255)
                                                                                                                ccGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGATGGCGTGTACACTGCCGATCCATTCAAGGACCCGCATGCCGAGaa                                   >  2:947027/1‑128 (MQ=255)
                                                                                                                ccGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGATGGCGTGTACACTGCCGATCCATTCAAGGACCCGCATGCCGAGaa                                   <  1:947027/128‑1 (MQ=255)
                                                                                                                 cGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGCACTTCGTGCTTCCGTGGTGTTGAAGGCGACCAAGGTCGATGGCGTGTACACTGCCGATCCATTCAAGGACCCGCATGCCGAGAAATTCGATCACCTGACCTACGATg            <  1:1098548/150‑1 (MQ=255)
                                                                                                                  gTTCTTCCCCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGATGGCGTGTACACTGCCGATCCATTCAAGGACCCGCATGCCGAGAAATTCGATCACCTGACCTACGATGa           <  2:65082/150‑1 (MQ=255)
                                                                                                                        caccacCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGATGGCGTGTACACTGCCGATCCATTCAAGGACCCGCATGCCGAGAAATTCGATCACCTGACCTACGATGAGGTTCt     >  1:947035/1‑150 (MQ=255)
                                                                                                                           caccGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGATGGCGTGTACACTGCCGATCCATTCAAGGACCCGCATGCCGAGAAATTCGATCACCTGACCTACGATGAGGTTCTGGa  >  1:687175/1‑150 (MQ=255)
                                                                                                                                     gcgGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGGAGGc                                                                                            >  2:802751/1‑50 (MQ=255)
                                                                                                                                     gcgGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGGAGGc                                                                                            <  1:802751/50‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                            
GGTCATGTCGGCCATTTCCATGGTGGGTGTCACCGATCATTACGATCGCCGCAAAGCTATTCGTCACCTCAACTCCGGGGATGTGGTAATTTTCTCCGCCGGTACCGGCAACCCGTTCTTCACCACCGACTCCGCGGCCTGCCTGCGCGCCATCGAAATCGATGCTGACGTGGTGCTGAAGGCGACCAAGGTCGATGGCGTGTACACTGCCGATCCATTCAAGGACCCGCATGCCGAGAAATTCGATCACCTGACCTACGATGAGGTTCTGGA  >  NC_002947/1787574‑1787846

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