Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,528,353 T→A 13.4% intergenic (+28/‑48) purB → / → PP_4017 adenylosuccinate lyase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,528,3530TA13.4% 120.7 / 11.1 52intergenic (+28/‑48)purB/PP_4017adenylosuccinate lyase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (19/26);  new base A (4/3);  total (23/29)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

GCTGCTGACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGCTGCTGGTACGC  >  NC_002947/4528211‑4528499
                                                                                                                                              |                                                                                                                                                  
gctgctGACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCgg                                                                                                                                              >  5:165835/1‑149 (MQ=255)
  tgctgACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGc                                                                                                                                            >  5:61941/1‑149 (MQ=255)
  tgctgACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGc                                                                                                                                            >  3:6831/1‑149 (MQ=255)
    ctgACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCg                                                                                                                                           <  7:119531/148‑1 (MQ=255)
         cTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGTGCGTTTTTa                                                                                                                                     <  5:184671/149‑1 (MQ=255)
         cTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTa                                                                                                                                     >  3:167956/1‑149 (MQ=255)
         cTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGCCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTa                                                                                                                                     >  1:260676/1‑149 (MQ=255)
            cATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTAtt                                                                                                                                   <  3:85041/148‑1 (MQ=255)
              tCGACGGCCTGGAAATGCCAGCCGAAGCCAAGTCCTATCTAAAGCTGCTTACCCCCTCTACCTAAATCGGTAAATCGGCATACCAGGCCAAACGCATCTAATCTTAGCGTCTCGTAAAACGCCCGGCCTGGCCGGGCGTTTTTATTccc                                                                                                                                <  4:196498/149‑1 (MQ=255)
                           cATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTAc                                                                                                                   <  8:32649/149‑1 (MQ=255)
                                  gCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGtttttt                                                                                                            <  1:290758/149‑1 (MQ=255)
                                   ccGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTc                                                                                                           <  3:92753/149‑1 (MQ=255)
                                       aGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATa                                                                                                       <  6:61941/149‑1 (MQ=255)
                                        tccAAGGGCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGGCCATGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCGGTCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTTCGTTTTTTCTATAg                                                                                                      <  2:93702/148‑1 (MQ=255)
                                           aaGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGtt                                                                                                   >  7:161861/1‑149 (MQ=255)
                                           aaGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACCCCCGGCCAGCCCGGACg                                                                                                                                           >  8:241487/1‑109 (MQ=255)
                                                       aaGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATcc                                                                                       >  3:210710/1‑149 (MQ=255)
                                                                 cccccGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCa                                                                             <  4:6831/149‑1 (MQ=255)
                                                                   cccGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACt                                                                           >  4:328829/1‑149 (MQ=255)
                                                                   ccagCTACCTACATCGGTAACGCGGCAGCCCATTCCAAACGCATATAAGCTGAGCGTAGCGTACAACGCCCTGCCTGGCCGGGCTTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTTAACATTAATCCTGATACTCCACt                                                                           <  3:116454/146‑1 (MQ=255)
                                                                     cGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGc                                                                         <  8:161861/149‑1 (MQ=255)
                                                                       ctacctacATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTTCAg                                                                       <  6:172078/149‑1 (MQ=255)
                                                                        tacctacATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAgc                                                                      >  6:342714/1‑149 (MQ=255)
                                                                           ctacATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTAACTTTTTTCAATAGTTTGACAATGAATCCTGATACTCCACAGCAgctt                                                                    >  6:29180/1‑147 (MQ=255)
                                                                                tCGGTAACGCGGGAGCCCATGCCAAACTCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTcggc                                                              <  6:162204/149‑1 (MQ=255)
                                                                                tCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATTATCCACTGCAGCTGCTcggc                                                              >  3:204907/1‑149 (MQ=255)
                                                                                tCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTcggc                                                              <  6:165835/149‑1 (MQ=255)
                                                                                tCGGTAACGCGGCAGCCCAGGACAAACGCATCTAAGCTGAGCGTCTCGTACAACGCCCGGCCTGGGCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATATGTTGAACATTAATCCTGTTACTCCACTGCAGCTGGTcggc                                                              <  2:314254/149‑1 (MQ=255)
                                                                                      aCGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCAtc                                                        <  7:75011/149‑1 (MQ=255)
                                                                                         cgtcAGCCCAGTCCAAACGCATCTAAGCTTAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCg                                                     <  6:91408/146‑1 (MQ=255)
                                                                                         cgGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCg                                                     <  5:342714/149‑1 (MQ=255)
                                                                                          ggCAGCCCAGGCCAAACACATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGGGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACAGCCGATGCTCGGCGGCATCTAgg                                                    >  5:43982/1‑149 (MQ=255)
                                                                                                  aGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGa                                                                                    <  4:190368/109‑1 (MQ=255)
                                                                                                  aGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGa                                                                                    >  3:190368/1‑109 (MQ=255)
                                                                                                      cAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCCGCTGCTCGGCGGCATCTCGGCCCGCGGATTcc                                        >  4:47706/1‑148 (MQ=255)
                                                                                                             aTCTAAGGTGAGCGTCGCGTAGAACGCCCGGCCTTGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGa                                 <  4:210710/149‑1 (MQ=255)
                                                                                                                  aGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGactac                             <  3:328829/148‑1 (MQ=255)
                                                                                                                     tGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGc                         <  2:260676/149‑1 (MQ=255)
                                                                                                                       aGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTAGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTgg                          >  8:153972/1‑146 (MQ=255)
                                                                                                                             gcgTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGaagaag                 >  1:39542/1‑149 (MQ=255)
                                                                                                                                tACAACGCCCGGCCTGGCCTGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAAGCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCg              <  5:29180/149‑1 (MQ=255)
                                                                                                                                tACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCg              >  6:134003/1‑149 (MQ=255)
                                                                                                                                 aaaaCGCCCGGCCAGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAgaagaa                  >  4:187481/3‑144 (MQ=255)
                                                                                                                                 aaaaCGCCCGGCCAGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAgaagaa                  <  3:187481/142‑1 (MQ=255)
                                                                                                                                  aaaCGCCCGGCCAGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTg                           <  4:304234/133‑1 (MQ=255)
                                                                                                                                  aaaCGCCCGGCCAGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTg                           >  3:304234/2‑134 (MQ=255)
                                                                                                                                   aaCGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTc                                                                               <  2:311750/81‑1 (MQ=255)
                                                                                                                                    aCGCCAGGCCAGGCAGGGCGTTTTTATTCCCGTACGAAAATTACTTTTTTTCAATAGGTTGAACATGAATACTGATACTCCACTGCAGCTGCTCGGCGGCATATCGGCCCGCGAATTCATGAGCGACTACTGGCAGAAGAAGCCcctac          >  6:116816/1‑147 (MQ=255)
                                                                                                                                      gaccGGCCAGGCCGGGCGTTTTTATTCCCGTACGTAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGATCGGCGGAATCTCGGCCAGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCTcttctg        <  5:116816/147‑1 (MQ=255)
                                                                                                                                        ccGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGctgctg        >  1:167149/1‑147 (MQ=255)
                                                                                                                                        ccGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGctgctg        <  2:167149/147‑1 (MQ=255)
                                                                                                                                             cTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGCTGCTGGTACGc  >  3:225626/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
GCTGCTGACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGCTGCTGGTACGC  >  NC_002947/4528211‑4528499

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑