Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,621,712 T→C 10.5% intergenic (+305/+4259) PP_5654 → / ← PP_4091 phosphotriesterase/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,621,7120TC10.5% 99.3 / 3.2 38intergenic (+305/+4259)PP_5654/PP_4091phosphotriesterase/transposase
Reads supporting (aligned to +/- strand):  ref base T (14/20);  new base C (2/2);  total (16/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

CTAGCGCAATGTTAAATTACGAGCGCACCACACAAACGCACCTAGCATCCCAAGAAACATAGAAAAAAACGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAACACATCAAATA  >  NC_002947/4621564‑4621859
                                                                                                                                                    |                                                                                                                                                     
cTAGCGCAATGTTAAATTACGAGCGCACCACACAAACGCACCTAGCATCCCAAGAAACATAGAAAAAAACGTTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAAt                                                                                                                                                       <  6:132650/149‑1 (MQ=255)
          gTTAAATTACGAGCGCACCACACAAACGCACCTAGCATCCCAAGAAACATAGAAAAAAACGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAAtt                                                                                                                                             >  3:299198/1‑149 (MQ=255)
              aaTTACGAGCGCACCACACAAACGCACCTAGCTAACCAAGAAACATAGAAAAAAACGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTa                                                                                                                                       <  8:5151/149‑1 (MQ=255)
                   cGAGCGCACCACACAAACGCACCTAGCATCCCAAGAAACATAGAAAAAAACGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACAt                                                                                                                                  <  7:55343/149‑1 (MQ=255)
                    gAGCGCACCACACAAACGCACCTAGCATCCCAAGAAACATAGAAAAAAACGGTTCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACTACTGGCCTCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATc                                                                                                                                 <  3:237598/149‑1 (MQ=255)
                       cgcACCACACAAACGCACCTAGCATCCCAAGAAACATAGAAAAAAACGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAg                                                                                                                              <  4:299198/149‑1 (MQ=255)
                          accacACAAACGCACCTAGCATCCCAAGAAACATAGAAAAAAACGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTc                                                                                                                                                    <  3:31030/126‑1 (MQ=255)
                          accacACAAACGCACCTAGCATCCCAAGAAACATAGAAAAAAACGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTc                                                                                                                                                    >  4:31030/1‑126 (MQ=255)
                          accacACAAACGCACCTAGCATCCCAAGAAACATAGAAAAAAACGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAgtg                                                                                                                            >  2:94458/1‑148 (MQ=255)
                           ccacACAAACGCACCTAGCATCCCAAGAAACATAGAAAAAAACGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTa                                                                                                                          <  4:53250/149‑1 (MQ=255)
                               acaAACGCACCTAGCATCCCAAGAAACATAGAAAAAAACGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAgt                                                                                                                             <  2:184719/142‑1 (MQ=255)
                               acaAACGCACCTAGCATCCCAAGAAACATAGAAAAAAACGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAgt                                                                                                                             >  1:184719/1‑142 (MQ=255)
                                                  cAAGAAACATAGAAAAAAACGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAg                                                                                                               <  3:115601/137‑1 (MQ=255)
                                                   aaGAAACATAGAAAAAAACGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGc                                                                                                                      >  2:247997/1‑129 (MQ=255)
                                                   aaGAAACATAGAAAAAAACGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGc                                                                                                                      <  1:247997/129‑1 (MQ=255)
                                                              aaaaaaaCGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGc                                                                                         >  1:64311/1‑147 (MQ=255)
                                                                            cAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATTTCGCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGgc                                                                           <  1:94458/149‑1 (MQ=255)
                                                                             aaCAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAACGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAgcgc                                                                                                   >  2:275556/1‑122 (MQ=255)
                                                                             aaCAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAACGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAgcgc                                                                                                   >  2:280063/1‑122 (MQ=255)
                                                                             aaCAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAACGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAgcgc                                                                                                   <  1:280063/122‑1 (MQ=255)
                                                                             aaCAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAACGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAgcgc                                                                                                   <  1:275556/122‑1 (MQ=255)
                                                                                  aTAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGc                                                                   >  3:260020/1‑149 (MQ=255)
                                                                                          tGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTa                                                                                                                                       <  4:263372/73‑1 (MQ=255)
                                                                                          tGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTa                                                                                                                                       >  3:263372/1‑73 (MQ=255)
                                                                                                 aCTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCgg                                                     <  1:125240/148‑1 (MQ=255)
                                                                                                         ccTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGAt                                            >  6:238612/1‑149 (MQ=255)
                                                                                                          cTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATc                                           >  5:174368/1‑149 (MQ=255)
                                                                                                          cTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATc                                           <  6:174368/149‑1 (MQ=255)
                                                                                                           tCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATgg                                                                            <  2:109115/115‑1 (MQ=255)
                                                                                                                   acaacTGGCCGCGAACTGTGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCcgc                                  <  4:309998/149‑1 (MQ=255)
                                                                                                                     aacTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATAGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGa                                >  1:139037/1‑149 (MQ=255)
                                                                                                                           ccGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCt                          <  2:64311/149‑1 (MQ=255)
                                                                                                                            cgcgAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTc                         <  4:260020/149‑1 (MQ=255)
                                                                                                                            cgcgAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTc                         <  1:9885/149‑1 (MQ=255)
                                                                                                                                     ggAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGc                >  8:93482/1‑149 (MQ=255)
                                                                                                                                        gCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGaa             <  5:238612/149‑1 (MQ=255)
                                                                                                                                            atatGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAacac          <  2:139037/148‑1 (MQ=255)
                                                                                                                                            atatGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCGGCCGCCTCATCGAACCCCGCGATCCCCTCGGCCGCAGCGAaccc          >  2:83218/1‑146 (MQ=255)
                                                                                                                                                   aTGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAACACATCAAATa  >  4:224415/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTAGCGCAATGTTAAATTACGAGCGCACCACACAAACGCACCTAGCATCCCAAGAAACATAGAAAAAAACGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAACACATCAAATA  >  NC_002947/4621564‑4621859

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