Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,623,195 T→C 8.8% intergenic (+1788/+2776) PP_5654 → / ← PP_4091 phosphotriesterase/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,623,1950TC 8.8% 118.9 / 2.3 46intergenic (+1788/+2776)PP_5654/PP_4091phosphotriesterase/transposase
Reads supporting (aligned to +/- strand):  ref base T (25/16);  new base C (2/2);  total (28/18)
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.97e-01

AAGTGCTGCCACTGCTGGACGAGGCCTATCGCTACGATGCCAGCGACAACCTGATTGCCGAGGTACTGACCCAGACCCAGCGCCGCGGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCATGACCTACCCCACACGGGCAAGACTTATTCGGGACACAACCGTTACGGCTACGACCTCAACGAGCAACTGCAAACCATCCAGCAGTCACGCCCCAACTGGCAAGCCACGCAGGTCGAGGATTTCAAGTACGACAAGGC  >  NC_002947/4623047‑4623332
                                                                                                                                                    |                                                                                                                                         
aaGTGCTGCCACTGCTGGACGAGGCCTATCGCTACGATGCCAGCGACAACCTGATTGCCGAGGTACTGACCCAGACCCAGCGCCGCGGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCAt                                                                                                                                           <  8:58079/149‑1 (MQ=255)
    gctgcCACTGCTGGACGAGGCCTATCGCTACGATGCCAGCGACAACCTGATTGCCGAGGTACTGACCCAGACCCAGCGCCGCGGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCATGAc                                                                                                                                        >  2:111617/1‑148 (MQ=255)
     ctgcCACTGCTGGACGAGGCCTATCGCTACGATGCCAGCGACAACCTGATTGCCGAGGTACTGACCCAGACCCAGCGCCGCGGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCATGAcc                                                                                                                                       >  2:134106/1‑148 (MQ=255)
     ctgcCACTGCTGGACGAGGCCTATCGCTACGATGCCAGCGACAACCTGATTGCCGAGGTACTGACCCAGACCCAGCGCCGCGGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCATGAcc                                                                                                                                       >  6:272051/1‑148 (MQ=255)
         cACTGCTGGACGAGGCCTATCGCTACGATGCCAGCGACAACCTGATTGCCGAGGTACTGACCCAGACCCAGCGCCGCGGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCATGACCTAccc                                                                                                                                  <  1:111617/149‑1 (MQ=255)
            tgctgGACGAGGCCCATCGCTACGATGCCAGCGACAACCTGATTGCCGAGGTACTGACCCAGACCCAGCGCAGCGGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCATGCCCTACCCcac                                                                                                                               >  7:276265/1‑149 (MQ=255)
               tgGACGAGGCCTATCGCTACGATGCCAGCGACAACCTGATTGCCGAGGTACTGACCCAGACCCAGCGCCGCGGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCACGACCTACCCCACACg                                                                                                                            >  6:136786/1‑149 (MQ=255)
                       gCCTATCGCTACGATGCCAGCGACAACCTGATTGCCGAGGTACTGACCCAGACCCAGCGCCGCGGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCATGACCTACCCCACACGGGCAAGAc                                                                                                                    >  3:96778/1‑149 (MQ=255)
                            tCGCTACGATGCCAGCGACAACCTGATTGCCGAGGTACTGACCCAGACCCAGCGCCGCGGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCATGACCTACCCCACACGGGCAAGACTTAtt                                                                                                               <  5:272051/149‑1 (MQ=255)
                            acGCTACGAAGCCAGCGACAACCTGATTGCCGAGGTACTGACCCAGACCCAGCGCCGCGGGATGACCAACGCGGCCAATGACGAAAACGCGCACCTCGAACAGATCATCGGCCGCTTTCATGACCTACCCCACACGGGCAAGACTTa                                                                                                                 <  8:225696/146‑1 (MQ=255)
                              gCTACGATGCCAGCGACAACCTGATTGCCGAGGTACTGACCCAGACCCAGCGCCGCGGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCATGACCTACCCCACACGGGCAAGACTTATTCg                                                                                                             <  4:96778/149‑1 (MQ=255)
                               cTACGATGCCAGCGACAACCTGATTGCCGAGGTACTGACCCAGACCCAGCGCCGCGGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCATGACCTACCCCACACGGGCAAGACTTATTCgg                                                                                                            >  5:172323/1‑149 (MQ=255)
                                tACGATGCCAGCGACAACCTGATTGCCGAGGTACTGACCCAGACCCAGCGCCGCGGGATGACCAACGCGGACAATGACGAAAACGAGCATCTCGACCAGATCATCGGCCGATTTCAGGACCTACCCCACACGGGCAAGACTTATTAggg                                                                                                           >  7:287813/1‑149 (MQ=17)
                                tACGATCCCAGCGACAACCTGATTGCCGAGGTACTGACCCAGACCCAGCGCCTCGGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTAGAACAGATCATCGGCCGCTTTCATGACCTACCCCACACGGGCAAGACTTATTCgg                                                                                                            >  2:32709/1‑148 (MQ=255)
                                  cGATGCCAGCGACAACCTGATTGCCGAGGTACTGACCCAGACCCAGCGCCGCGGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCATGACCTAc                                                                                                                                    <  2:208764/122‑1 (MQ=255)
                                  cGATGCCAGCGACAACCTGATTGCCGAGGTACTGACCCAGACCCAGCGCCGCGGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCATGACCTAc                                                                                                                                    >  1:208764/1‑122 (MQ=255)
                                             aCAACCTGATTGCCGAGGTACTGACCCAGACCCAGCGCCGCGGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCATGACCTACCCCACACGGGCAAGACTTATTCGGGACACAACCGTTAc                                                                                              <  1:134106/149‑1 (MQ=255)
                                                 ccTGATTGCCGAGGTACTGACCCAGACCCAGCGCCGCGGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCATGACCTACCCCACACGGGCAAGACTTATTCGGGACACAACCGTTACGGc                                                                                           >  7:172528/1‑148 (MQ=255)
                                                 ccTGATTGCCGAGGTACTGACCCAGACCCAGCGCCGCGGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCACGACCTACCCCACACGGGCAAGACTTATTCGGGACACAACCGTTACGGc                                                                                           >  4:65933/1‑148 (MQ=255)
                                                   tGATTGCCGAGGTACTGACCCAGACCCAGCGCCGCGGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCATGACCTGCCCCACACGGGCAAGACTTATTCGGGACACAACCGTTACGGCTAc                                                                                        >  4:177119/1‑149 (MQ=255)
                                                   tGATTGCCGAGGTAATGACCCAGACCCAGCGCCGCGGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGACGCTTTCATGACCTACCCCACACGGGCAAGACTTATTCGGGACACAACCGTTACGGCTAc                                                                                        >  8:166814/1‑149 (MQ=255)
                                                        gCCGAGGTACTGACCCAGACCCAGCGCCGCGGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCATGACCTACCCCACACGGGCAAGACTTATTCGGGACACAACCGTTACGGCTACGAcc                                                                                    <  6:172323/148‑1 (MQ=255)
                                                               tACTGACCCAGACCCAGCGCCGCGGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCACGACCTACCCCACACGGGCAAGACTTATTCGGGACACAACCGTTACGGCTACGACCTCAACGAg                                                                            <  5:136786/149‑1 (MQ=255)
                                                                aCTGACCCAGACCCAGCGCCGCGGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCATGACCTACCCCACACGGGCAAGACTTATTCGGGACACAACCGTTACGGCTACGACCTCAACGAg                                                                            >  3:17127/1‑148 (MQ=255)
                                                                      ccagacccagCGCCGCGGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCATGACCTACCCCACACGGGCAAGACTTATTCGGGACACAACCGTTACGGCTACGACCTCAAc                                                                               >  5:149560/1‑139 (MQ=255)
                                                                      ccagacccagCGCCGCGGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCATGACCTACCCCACACGGGCAAGACTTATTCGGGACACAACCGTTACGGCTACGACCTCAAc                                                                               <  6:149560/139‑1 (MQ=255)
                                                                        agacccagCGCCGCGGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCATGACCTACCCCACACGGGCAAGACTTATTCGGGGCACAACCGTTACGGCTACGACCTCAACGAGCAACTGCa                                                                    >  5:78646/1‑148 (MQ=255)
                                                                         gacccagCGCCGCGGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCACGACCTACCCCACACGGGCAAGACTTATTCGGGACACAACCGTTACGGCTACGACCTCAACGAGCAACTGCaaa                                                                  <  3:65933/149‑1 (MQ=255)
                                                                                 gccgcGGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCATGACCTACCCCACACGGGCAAGACTTATTCGGGACACAACCGTTACGGCTACGACCTCAACGAGCAACTGCAAACCATCcag                                                          >  5:171586/1‑149 (MQ=255)
                                                                                     ctggATGACCAACGAGGACAATGACGAAAACGCGGATCTCGAACACATCATCGTCCGCTTTCATGACCTACCCCACACGGGCAAGACTTATTCGGGACACAACCGTTACGGCTACGACCTAAACTAGCAACTGCAAACCATCcagca                                                        <  8:287813/145‑1 (MQ=18)
                                                                                     cgGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCATGACCTACCCCACACGGGCAAGACTTATTCGGGACACAACCGTTACGGCTACGACCTCAACGAGCAACTGCAAACCATCCAGCAGt                                                      <  6:171586/149‑1 (MQ=255)
                                                                                     cgGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTAGAACAGATCATCGGCCGCTTTCATGACCTACCCCACACGGGCAAGACTTATTCGGGACACAACCGTTACGGCTACGACCTCAACGAGCAACTGCAAACCATCCAGCAGt                                                      <  1:32709/149‑1 (MQ=255)
                                                                                        gATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCATGACCTACCCCACACGGGCAAGACTTATTCGGGACACAACCGTTACGGCTACGACCTCAACGAGCAACTGCAAACCATc                                                             >  3:6228/1‑139 (MQ=255)
                                                                                        gATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCATGACCTACCCCACACGGGCAAGACTTATTCGGGACACAACCGTTACGGCTACGACCTCAACGAGCAACTGCAAACCATc                                                             <  4:6228/139‑1 (MQ=255)
                                                                                           gACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCATGACCTACCCCACACGGGCAAGACTTATTCGGGACACAACCGTTACGGCTACGACCTCAACGAGCAACTGCAAACCATc                                                             >  8:210601/1‑136 (MQ=255)
                                                                                              cAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCATGACCTACCCCACACGGGCAAGACTTATTCGGGACACAACCGTTACGGCTACGACCTCAACGAGCAACTGCAAACCATCCAGCAGTCACGcccc                                              >  7:36788/1‑148 (MQ=255)
                                                                                                          tGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCATGACCTACCCCACACGGGCAAGACTTATTCGGGACACAACCGTTACGGCTACGACCTCAACGAGCAACTGCAAACCATCCAGCAGTCACGCCCCAACTGGCAAGCCa                                 >  1:50617/1‑149 (MQ=255)
                                                                                                             cGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCATGACCTACCCCACACGGGCAAGACTTATTCGGGACACAACCGTTACGGCTACGACCTCAACGAGCAACTGCAAACCATCCAGCAGTCACGCCCCAACTGGCAAGCCACGc                              <  7:166814/149‑1 (MQ=255)
                                                                                                                 atcgctCATCTCGAAGAGATCATCGGCCTCTTTCATGACCTACCCCACTCGGGCAAGACTTATTCGGGACACAACCGTAACGGCTACGACCTCAACGAGCATCTGCAACCCATCCAGCAGTCACGCCCCATCTGGCAAGCCACGCAGGt                          <  6:179340/147‑1 (MQ=17)
                                                                                                                      gcATCTTGAACAGATCATCGGCCGCTTTCATGACCTACCCCACACGGGCAAGACTTATTCGGGACACAACCGTTACGGCTACGACCTCAACGAGCAACTGCAAACCATCCAGCAGTCACGCCCCAACTGGCAAGCCACGCAGGTCGAg                      >  8:42413/1‑148 (MQ=255)
                                                                                                                         tctcGAACAGATCATCGGCCGCTTTCATGACCTACCCCACACGGGCAAGACTTATTCGGGACACAACCGTTACGGCTACGACCTCAACGAGCAACTGCAAACCATCCAGCAGTCACGCCCCAACTGGCAAGCCACGCAGGTCGAGGAtt                  >  5:132481/1‑149 (MQ=255)
                                                                                                                                   atcatcGGCCGCTTTCATGACCTACCCCACACGGGCAAGACTTATTCGGGACACAACCGTTACGGCTACGACCTCAACGAGCAACTGCAAACCATCCAGCAGTCACGCCCCAACTGGCAAGCCACGCAGGTCGAGGATTTCAAGTACGa        >  4:63161/1‑149 (MQ=255)
                                                                                                                                       tcGGCCGCTTTCATGACCTACCCCACACGGGCAAGACTTATTCGGGACACAACCGTTACGGCTACGACCTCAACGAGCAACTGCAAACCATCCAGCAGTCACGCCCCAACTGGCAAGCCACGCAGGTCGAGGATTTCAAGTACGACAAg    >  4:58114/1‑149 (MQ=255)
                                                                                                                                         ggCCGCTTTCATGACCTACCCCACACGGGCAAGACTTATTCGGGACACAACCGTTACGGCTACGACCTCAACGAGCAACTGCAAACCATCCAGCAGTCACGCCCCAACTGGCAAGCCACGCAGGTCGAGGATTTCAAGTACGACAAGGc  >  2:201003/1‑149 (MQ=255)
                                                                                                                                         ggCCGCTTTCATGACCTACCCCACACGGGCAAGACTTATTCGGGACACAACCGTTACGGCTACGACCTCAACGAGCAACTGCAAACCATCCAGCAGTCACGCCCCAACTGGCAAGCCACGCAGGTCGAGGATTTCAAGTACGACAAGGc  <  8:276265/149‑1 (MQ=255)
                                                                                                                                                   aTGACCTACCCCACACGGGCAAGACTTATTCGGGACACAACCGTTACGGCTACGACCTCAACGAGCAACTGCAAACCATCCAGCAGTCACGCCCCAACTGGCAAGCCACGCAGGTCGAGGATTTCAAGTACGACaa     <  6:312814/136‑1 (MQ=255)
                                                                                                                                                   aTGACCTACCCCACACGGGCAAGACTTATTCGGGACACAACCGTTACGGCTACGACCTCAACGAGCAACTGCAAACCATCCAGCAGTCACGCCCCAACTGGCAAGCCACGCAGGTCGAGGATTTCAAGTACGACaa     >  5:312814/1‑136 (MQ=255)
                                                                                                                                                    |                                                                                                                                         
AAGTGCTGCCACTGCTGGACGAGGCCTATCGCTACGATGCCAGCGACAACCTGATTGCCGAGGTACTGACCCAGACCCAGCGCCGCGGGATGACCAACGCGGCCAATGACGAAAACGCGCATCTCGAACAGATCATCGGCCGCTTTCATGACCTACCCCACACGGGCAAGACTTATTCGGGACACAACCGTTACGGCTACGACCTCAACGAGCAACTGCAAACCATCCAGCAGTCACGCCCCAACTGGCAAGCCACGCAGGTCGAGGATTTCAAGTACGACAAGGC  >  NC_002947/4623047‑4623332

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