Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,623,823 C→T 9.6% intergenic (+2416/+2148) 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,8230CT 9.6% 110.9 / 2.5 42intergenic (+2416/+2148)PP_5654/PP_4091phosphotriesterase/transposase
Reads supporting (aligned to +/- strand):  ref base C (20/17);  new base T (2/2);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.68e-01

TAGCTTGTATGTGTATGCCACTCCGGACAGCTATGAGCCGCTAGCACGGATCGATGGGTCGATGGGTAGGGAGGTGGTTCAGTATTTTCATACCAACCTGGCCGGGCTTCCCGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGAGGGCTGGGGGAAGACGCGGGATGAGTGGCACAGTCCACACCAGGCGAGAGAACAGAACCTGCGATATCAAGGTCAGTATGTTGACAGAGAAATAAGCCTCCATTACAACACCTTCAGGTTCTACGATTTTACCA  >  NC_002947/4623675‑4623970
                                                                                                                                                    |                                                                                                                                                   
tAGCTTGTATGTGTATGCCACTCCGGACAGCTATGAGCCGCTAGCACGGATCGATGGGTCGATGGGTAGGGAGGTGGTTCAGTATTTTCATACCAACCTGGCCGGGCTTCCCGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGc                                                                                                                                                     >  2:24329/1‑149 (MQ=255)
               tgCCACTCCGGACAGCTATGAGCCGCTAGCACGGATCGATGGGTCGATGGGTAGGGAGGTGGTTCAGTATTTTCATACCAACCTGGCCGGGCTTCCCGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGAgg                                                                                                                                      <  8:196890/149‑1 (MQ=255)
                 ccACTCCGGACAGCTATGAGCCGCTAGCACGGATCGATGGGTCGATGGGTAGGGAGGTGGTTCAGTATTTTCATACCAACCTGGCCGGGCTTCCCGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGAGGGc                                                                                                                                    >  4:112938/1‑149 (MQ=255)
                  cACTCCGGACAGCTATGAGCCGCTAGCACGGATCGATGGGTCGATGGGTAGGGAGGTGGTTCAGTATTTTCATACCAACCTGGCCGGGCTTCCCGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGAGGGCt                                                                                                                                   >  3:66398/1‑149 (MQ=255)
                       cGGACAGCTATGAGCCGCTAGCACGGATCGATGGGTCGATGGGTAGGGAGGTGGTTCAGTATTTTCATACCAACCTGGCCGGGCTTCCCGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGAGGGCTgggg                                                                                                                               <  2:120226/148‑1 (MQ=255)
                        ggACAGCTATGAGCCGCTAGCACGGATCGATGAGTCGATGGGTAGGGAGGTGGTTCAGTATTTTCATACCAACCTGGCCGGGCTTCCCGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGAGGGCTGGGGGa                                                                                                                             <  7:250499/149‑1 (MQ=255)
                                        cTAGCACGGATCGATGGGTCGATGGGTAGGGAGGTGGTTCAGTATTTTCATACCAACCTGGCCGGGCTTCCCGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGAGGGCTGGGGGAAGACGCGGGATGAGTg                                                                                                             <  1:24329/149‑1 (MQ=255)
                                           gCACGGATCGATGGGTCGATGGGTAGGGAGGTGGTTCAGTATTTTCATACCAACCTGGCCGGGCTTCCCGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGAGGGCTGGGGGAAGACGCGGGATGAGTgg                                                                                                            <  6:220584/147‑1 (MQ=255)
                                              cGGATCGATGGGTCGATGGGTAGGGAGGTGGTTCAGTTTTTTCATACCAACCTGGCCGGGCTTCCCGAGCAACTGACTGACGGGAATGGCAATGCGATTTGGCGTAACGATTTTGAGGGCTGGGGGAAGACGCGGAAGGAGTGGcccgg                                                                                                       >  6:100892/1‑145 (MQ=255)
                                                  tcgatgggtcgatgggTAGGGAGGTGGTTCAGTATTTTCATACCAACCTGGCCGGGCTTCCCGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGAGGGCTGGGGGAAGACGCGGGATGAGTGGCACAGTCca                                                                                                   <  4:66398/149‑1 (MQ=255)
                                                     atgggtcgatgggtAGGGAGGTGGTTCAGTTTTTTCATACCAACCTGGCCGGGCTTCCCGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGAGGGCTGGGGGAAGACGCGGGAt                                                                                                                  >  1:28168/1‑131 (MQ=255)
                                                     atgggtcgatgggtAGGGAGGTGGTTCAGTATTTTCATACCAACCTGGCCGGGCTTCCCGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGAGGGCTGGGGGAAGACGCGGGAt                                                                                                                  <  2:28168/131‑1 (MQ=255)
                                                          tcgatgggtAGGGAGGTGGTTCAGTATTTTCATACCAACCTGGCCGGGCTTCCCGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGTGTAACGATTATGAGGGCTGGGGGAAGACGCggg                                                                                                                    <  8:73359/124‑1 (MQ=255)
                                                          tcgatgggtAGGGAGGTGGTTCAGTATTTTCATACCAACCTGGCCGGGCTTCCCGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGTGTAACGATTATGAGGGCTGGGGGAAGACGCggg                                                                                                                    <  8:69030/124‑1 (MQ=255)
                                                          tcgatgggtAGGGAGGTGGTTCAGTATTTTCATACCAACCTGGCCGGGCTTCCCGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGTGTAACGATTATGAGGGCTGGGGGAAGACGCggg                                                                                                                    >  7:73359/1‑124 (MQ=255)
                                                          tcgatgggtAGGGAGGTGGTTCAGTATTTTCATACCAACCTGGCCGGGCTTCCCGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGTGTAACGATTATGAGGGCTGGGGGAAGACGCggg                                                                                                                    >  7:69030/1‑124 (MQ=255)
                                                                  tAGGGAGGTGGTTCAGTATTTTCATACCAACCTGGCCGGGCTTCCCGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGAGGGCTGGGGGAAGACGCGGGATGAGTGGCACAGTCCACACCAGGCGAGagaac                                                                                   <  6:93134/149‑1 (MQ=255)
                                                                     ggAGGTGGTTCAGTATTTTCATACCAACCTGGCCGGGCTTCCCGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGAGGGCTGGGGGAAGAcgcg                                                                                                                      >  2:288656/1‑111 (MQ=255)
                                                                     ggAGGTGGTTCAGTATTTTCATACCAACCTGGCCGGGCTTCCCGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGAGGGCTGGGGGAAGAcgcg                                                                                                                      <  1:288656/111‑1 (MQ=255)
                                                                      gAGGTGGTTCAGTATTTTCATACCAACCTGGCCGGGCTTCCCGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGAGGGCTGGGGGAAGACGCGGGATGAGTGGCACAGTCCACACCAGGCGAGagaacag                                                                                 >  1:184393/1‑147 (MQ=255)
                                                                      gAGGTGGTTCAGTATTTTCATACCAACCTGGCCGGGCTTCCCGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGAGGGCTGGGGGAAGACGCGGGATGAGTGGCACAGTCCACACCAGGCGAGagaacag                                                                                 <  2:184393/147‑1 (MQ=255)
                                                                        ggtggtTCAGTATTTTCATACCAACCTGGCCGGGCTTCCCGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGAGGGCTGGGGGAAGACGCGGGATGAGTGGCACAGTCCACACCAGGCGAGAGAACAGAAcc                                                                             >  3:31586/1‑149 (MQ=255)
                                                                           ggtTCAGTATTTTCATACCAACCTGGCCGGGCTTCCCGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGAGGGCTGGGGGAAGACGCGGGATGAGTGGCACAGTCCACACCAGGCGAGAGAACAGAACCTGc                                                                          >  6:239484/1‑149 (MQ=255)
                                                                                           aCCAACCTGGCCGGGCTTCCCGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGAGGGCTGGGGGAAGACGCGGGATGAGTGGCACAGTCCACACCAGGCGAGAGAACAGAACCTGCGATATCAAGGTCAGTa                                                          >  2:167055/1‑149 (MQ=255)
                                                                                             cAACCTGGCCGGGCTTCCCGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGAGGGCTGGGGGAAGACGCGGGATGAGTGGCACAGTCCACACCAGGCGAGAGAACAGAACCTGCGATATCAAGGTCAGTATg                                                        >  8:217403/1‑149 (MQ=255)
                                                                                                   ggCCGGGCTTCCCGATCAACTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGAGGGCTGGGGGAAGACGCGGGATGAGTGGCACAGACCACACCAGGCGAGAGAACAGAACCTGCGATCTCAAGGTGAGTATGTTGa                                                    >  4:163926/1‑147 (MQ=255)
                                                                                                      cGGGCTTCCCGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGAGGGCTGGGGGAAGACGCGGGATGAGTGGCACAGTCCACACCAGGCGAGAGAACAGAACCTGCGATATCAAGGTCAGTATGTTGACa                                                  >  6:293006/1‑146 (MQ=255)
                                                                                                           ttCCCGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGAGGGCTGGGGGAAGACGCGGGATGAGTGGCACAGTCCACACCAGGCGAGAGAACAGAACCTGCGATATCAAGGTCAGTATGTTGACAGAGAAAt                                           >  3:158428/1‑148 (MQ=255)
                                                                                                             cccGATCAACTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGACGGCTGGGGGAAGACGCGGGATGAGTGGCACAGTCCACACCAGGCGAGAGAACAGAACCTGCGATATCAAGGTCAGTATGTTGACAGAGAAATaa                                         <  5:100892/148‑1 (MQ=255)
                                                                                                             cccGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGAGGGCTGGGGGAAGACGCGGGATGAGTGGCACAGTCCACACCAGGCGAGAGAACAGAACCTGCGATATCAAGGTCAGTATGTTGACAGAGAAATAAg                                        >  8:43325/1‑149 (MQ=255)
                                                                                                             cccGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGAGGGCTGGGGGAAGACGCGGGATGAGTGGCACAGTCCACACCAGGCGAGAGAACAGAACCTGCGATATCAAGGTAAGTATGTTGACAGAGAAATAAg                                        >  3:110405/1‑149 (MQ=255)
                                                                                                              ccGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGAGGGCTGGGGGAAGACGCGGGATGAGTGGCACAGTCCACACCAGGCGAGAGAACAGAACCTGCGATATCAAGGTCAGTATGTTGACAGAGAAATAAg                                        <  3:112938/148‑1 (MQ=255)
                                                                                                                gAGCAACTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGAGGGCTGGGGGAAGACGCGGGATGAGTGGCACAGTCCACACCAGGCGAGAGAACAGAACCTGCGATATCAAGGTCAGTa                                                          <  6:297788/128‑1 (MQ=255)
                                                                                                                gAGCAACTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGAGGGCTGGGGGAAGACGCGGGATGAGTGGCACAGTCCACACCAGGCGAGAGAACAGAACCTGCGATATCAAGGTCAGTa                                                          >  5:297788/1‑128 (MQ=255)
                                                                                                                    aaCTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGAGGGCTGGGGGAAGACGCGGGATGAGTGGCACAGTCCACACCAGGCGAGAGAACAGAACCTGCGATATCAAGGTCAGTATGTTGACAGAGAAATAAGCCTCCAt                                 <  7:217403/149‑1 (MQ=255)
                                                                                                                           tgacGCGAATGGCAATGCGATTTGGCGTAACGATTATGAGGGCTGGGGGAAGACGCGGGATGAGTGGCACAGTCCACACCAGGCGAGAGAACAGAACCTGCCATATCAAGGTCAGTATGTTGACAGAGAAATAAGCCTCCATTacaaca                          >  6:257232/1‑149 (MQ=255)
                                                                                                                                   aTGGCAATGCGATTTGGAGTAACGATTATGATGGCTGGGGGAAGACTCGGGATGAGTGTCACAGTCCACACCATGCGAGAGAACAGAACCTGCGATATCAAGGTCAGTATGTTGACAGAGAAATAAGCCTCCATTACAACACCTTCAgg                  <  3:163926/149‑1 (MQ=255)
                                                                                                                                     ggCAATGCGATTTGGCGTAACGATTATGAGGGCTGGGGGAAGACGCGGGATGAGTGGCACAGTCCACACCAGGCGAGAGAACAGAACCTGCGATATCAAGGTCAGTATGTTGACAGAGAAATAAGCCTCCATTACAACACCTTCAGGt                 >  7:231746/1‑148 (MQ=255)
                                                                                                                                        aaTGCGATTTGGCGTAACGATTATGAGGGCTGGGGGAAGACGCGGGATGAGTg                                                                                                             >  7:174354/1‑53 (MQ=255)
                                                                                                                                        aaTGCGATTTGGCGTAACGATTATGAGGGCTGGGGGAAGACGCGGGATGAGTg                                                                                                             <  8:174354/53‑1 (MQ=255)
                                                                                                                                         aTGCGATTTGGCGTAACGATTATGAGGGCTGGGGGAAGACGCGGGATGAGTGGCACAGTCCACACCAGGCGAGAGAACAGAACCTGCGATATCAAGGTCAGTATGTTGACAGAGAAATAAGCCTCCATTACAACACCTTCAGGTTCTAc            >  6:127082/1‑149 (MQ=255)
                                                                                                                                               tttGGCGTAACGATTATGAGGGCTGGGGGAAGACGCGGGATGAGTGGCACAGTCCACACCAGGCGAGAGAACAGAACCTGCGATATCAAGGTCAGTATGTTGACAGAGAAATAAGCCTCCATTACAACACCTTCAGGTTCTACGAtttt      <  5:127082/149‑1 (MQ=255)
                                                                                                                                                   gCGTAACGATTATGAGGGCTGGGGGAAGACGCGGGATGAGTGGCACAGTCCACACCAGGCGAGAGAACAGAACCTGCGATATCAAGGTCAGTATGTTGACAGAGAAATAAGCCTCCATTACAACACCTTCAGGTTCTACGATTTTACCa  <  2:271444/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TAGCTTGTATGTGTATGCCACTCCGGACAGCTATGAGCCGCTAGCACGGATCGATGGGTCGATGGGTAGGGAGGTGGTTCAGTATTTTCATACCAACCTGGCCGGGCTTCCCGAGCAACTGACTGACGCGAATGGCAATGCGATTTGGCGTAACGATTATGAGGGCTGGGGGAAGACGCGGGATGAGTGGCACAGTCCACACCAGGCGAGAGAACAGAACCTGCGATATCAAGGTCAGTATGTTGACAGAGAAATAAGCCTCCATTACAACACCTTCAGGTTCTACGATTTTACCA  >  NC_002947/4623675‑4623970

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