Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,822,935 G→A 8.2% L136L (CTG→CTA rdgB → dITP/XTP pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,822,9350GA 8.2% 136.6 / 3.3 49L136L (CTG→CTArdgBdITP/XTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base G (23/22);  new base A (2/2);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

TCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGATCCGCTGTTCTGGGTACCGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGC  >  NC_002947/5822786‑5823079
                                                                                                                                                     |                                                                                                                                                
tCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCtg                                                                                                                                                  >  1:1125916/1‑150 (MQ=255)
       gcTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAg                                                                                                                                           <  1:162144/150‑1 (MQ=255)
         tATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGc                                                                                                                                         <  1:13794/150‑1 (MQ=255)
         tATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGc                                                                                                                                         <  1:440792/150‑1 (MQ=255)
             ccGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCtgt                                                                                                                                     >  1:1074610/1‑150 (MQ=255)
                    cAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACgg                                                                                                                              <  1:294677/150‑1 (MQ=255)
                     aaaGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACggg                                                                                                                             >  2:948515/1‑150 (MQ=255)
                           gACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATc                                                                                                                       <  2:1074610/150‑1 (MQ=255)
                                ggCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGTCACGGGCGCATCCTGtt                                                                                                                  <  2:618354/150‑1 (MQ=255)
                                 gCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGtt                                                                                                                  >  2:1099684/1‑149 (MQ=255)
                                 gCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGtt                                                                                                                  <  1:1099684/149‑1 (MQ=255)
                                 gCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTATGCGAAGGCCTGTGGCACGGGCGCATCCTGTTc                                                                                                                 >  2:992836/1‑150 (MQ=255)
                                 gCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTATGCGAAGGCCTGTGGCACGGGCGCATCCTGTTc                                                                                                                 >  2:392293/1‑150 (MQ=255)
                                    aacaacGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTATGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAg                                                                                                              <  1:392293/150‑1 (MQ=255)
                                    aacaacGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTATGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAg                                                                                                              <  1:992836/150‑1 (MQ=255)
                                      caacGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTgg                                                                                                                                   >  1:158108/1‑127 (MQ=255)
                                      caacGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTgg                                                                                                                                   <  2:158108/127‑1 (MQ=255)
                                       aacGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGcc                                                                                                           >  1:160839/1‑150 (MQ=255)
                                          gCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGgc                                                                                                                            >  1:205128/1‑130 (MQ=255)
                                          gCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGgc                                                                                                                            <  2:205128/130‑1 (MQ=255)
                                           ccAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAgcg                                                                                                       <  2:120933/150‑1 (MQ=255)
                                                    tCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGaa                                                                                                                                            >  1:193461/1‑104 (MQ=255)
                                                    tCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGaa                                                                                                                                            <  2:193461/104‑1 (MQ=255)
                                                     cGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAggg                                                                                             <  2:1125916/150‑1 (MQ=255)
                                                      gAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGc                                                                                            <  2:607653/150‑1 (MQ=255)
                                                                         cTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAg                                                                                                                                           <  1:789573/84‑1 (MQ=255)
                                                                         cTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAg                                                                                                                                           >  2:789573/1‑84 (MQ=255)
                                                                               cTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGATCCGCTGTTCTggg                                                                   >  2:996432/1‑150 (MQ=255)
                                                                                    cgcAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGATCCGCTGTTCTGGGTACCg                                                              <  2:932867/150‑1 (MQ=255)
                                                                                                    tCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGAt                                                                                <  2:90137/116‑1 (MQ=255)
                                                                                                    tCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGAt                                                                                >  1:90137/1‑116 (MQ=255)
                                                                                                    tCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGATCCGCTGTTCTGGGTACCGGAGCGCAACTGCTCCa                                              >  1:1103032/1‑150 (MQ=255)
                                                                                                      tGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGATCCGCTGTTCTGGGTACCGGAGCGCAACTGCTCCAgc                                            >  1:555005/1‑150 (MQ=255)
                                                                                                                cgcaccagcgcCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGATCCGCTGTTCTGGGTACCGGAGCGCAACTGCTCCAGCGCCGATCTAg                                  <  2:965663/143‑1 (MQ=255)
                                                                                                                 gcTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGATCCGCTGTTCTGGGTACCGGAGCGCAACTGCTCCAGCGCCGATCTAGc                                 >  1:202098/1‑150 (MQ=255)
                                                                                                                     gTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGATCCGCTGTTCTGGGTACCGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCca                             >  2:703590/1‑150 (MQ=255)
                                                                                                                          ccACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGATCCGCTGTTCTGGGTACCGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGa                        >  1:897271/1‑150 (MQ=255)
                                                                                                                               ccGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGATCCGCTGTTCTGGGTACCGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGa                   >  1:1198332/1‑150 (MQ=255)
                                                                                                                                cgacgaCCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGATCCGCTGTTCTGGGTACCGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGaa                  >  2:1132461/1‑150 (MQ=255)
                                                                                                                                 gacgacCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGATCCGCTGTTCTGGGTACCGGAGCGCAACTGCTCCAGCGCCGATATAGCCCCAGCAGACAAGAAc                 >  2:231629/1‑150 (MQ=255)
                                                                                                                                  ccgacCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGATCCGCTGTTCTGGGTACCGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAAcc                <  2:1198332/149‑1 (MQ=255)
                                                                                                                                       ccGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCTACCAGGGCTTTGGCTACGATCCGCTGTTCTGGGTACCGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTc           >  2:593288/1‑150 (MQ=255)
                                                                                                                                       ccGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGATCCGCTGtt                                                                        <  2:505705/89‑1 (MQ=255)
                                                                                                                                       ccGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGATCCGCTGtt                                                                        >  1:505705/1‑89 (MQ=255)
                                                                                                                                              cGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGATCCGCTGTTCTGGGTACCGGAGCGCAACTGCTCCAGCGCCGAt                                      >  2:154/1‑116 (MQ=255)
                                                                                                                                              cGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGATCCGCTGTTCTGGGTACCGGAGCGCAACTGCTCCAGCGCCGAt                                      <  1:154/116‑1 (MQ=255)
                                                                                                                                              cGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGATCCGCTGTTCTGGGTACCGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCAcc    >  2:287739/1‑150 (MQ=255)
                                                                                                                                                aTCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGATCCGCTGTTCTGGGTACCGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACcgc  <  2:202098/150‑1 (MQ=255)
                                                                                                                                                   cTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGATCCGCTGTTCTGGGTACCGGAGCGCAACTGCTCCAGCGCCGATCTAGccc                               <  2:100912/118‑1 (MQ=255)
                                                                                                                                                   cTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGATCCGCTGTTCTGGGTACCGGAGCGCAACTGCTCCAGCGCCGATCTAGccc                               >  1:100912/1‑118 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
TCGGCGCGCTATGCCGATGGCAAAGGTGACGCGGCAAACAACGCCAAGCTGCTCGAAGCCTTGAAAGACGTGCCTGAAGCTGAGCGCAGTGCACAGTTCGTCTGCGTTCTGGCGCTGGTGCGCCACGCCGACGACCCGTTGCCGATCCTGTGCGAAGGCCTGTGGCACGGGCGCATCCTGTTCGAGGCCAGCGGCGACCAGGGCTTTGGCTACGATCCGCTGTTCTGGGTACCGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGC  >  NC_002947/5822786‑5823079

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