Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,712,951 T→A 12.6% intergenic (+34/‑106) adk → / → tsaB adenylate kinase/endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,712,9510TA12.6% 88.4 / 4.6 38intergenic (+34/‑106)adk/tsaBadenylate kinase/endopeptidase
Reads supporting (aligned to +/- strand):  ref base T (16/17);  new base A (3/2);  total (19/19)
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

CTACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCACCGAAGCCTGTT  >  NC_002947/1712809‑1713099
                                                                                                                                              |                                                                                                                                                    
cTACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCggg                                                                                                                                                >  5:13420/1‑149 (MQ=255)
   ccAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGAGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCAAGCggg                                                                                                                                                <  3:44588/146‑1 (MQ=255)
          cTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGt                                                                                                                                      <  4:133091/149‑1 (MQ=255)
           tGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGGGGGCCGGTGTCggt                                                                                                                                     >  1:339943/1‑147 (MQ=255)
             tCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGttt                                                                                                                                    <  3:146957/148‑1 (MQ=255)
             gcggcggCCAATGGCGGCAAGCCTAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGt                                                                                                                                      <  2:339943/145‑1 (MQ=255)
                        aggcggcAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGcaca                                                                                                                        >  5:30933/2‑149 (MQ=255)
                          gcggcAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAgc                                                                                                                      >  7:284905/1‑149 (MQ=255)
                           cggcAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAgcg                                                                                                                     >  8:33434/1‑149 (MQ=255)
                            ggcAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCAAGCGGGCCGTTGTc                                                                                                                                        >  4:219189/1‑129 (MQ=255)
                            ggcAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCAAGCGGGCCGTTGTc                                                                                                                                        <  3:219189/129‑1 (MQ=255)
                                    gAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAg                                                                                                            <  6:30933/149‑1 (MQ=255)
                                         aCAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTa                                                                                                                                 >  4:272860/1‑123 (MQ=255)
                                         aCAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTa                                                                                                                                 <  3:272860/123‑1 (MQ=255)
                                           aGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACaa                                                                                                     >  1:61879/1‑149 (MQ=255)
                                               acaTCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCATGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAAcgcc                                                                                                 <  8:284905/149‑1 (MQ=255)
                                                           tCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGctct                                                                                      <  2:61879/148‑1 (MQ=255)
                                                           tCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCta                                                                                     <  7:33434/149‑1 (MQ=255)
                                                                   gTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCgctgct                                                                                           >  5:156623/1‑135 (MQ=255)
                                                                   gTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCgctgct                                                                                           <  6:156623/135‑1 (MQ=255)
                                                                   gTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCg                                                                             >  2:212107/1‑149 (MQ=255)
                                                                      gAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTg                                                                                                                                   <  2:101929/92‑1 (MQ=255)
                                                                      gAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTg                                                                                                                                   >  1:101929/1‑92 (MQ=255)
                                                                      gAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGCCACAACGCCGCTGCTGCTCTATACTGCCGctc                                                                          >  6:6628/1‑149 (MQ=255)
                                                                       aaGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGctct                                                                         <  6:13420/149‑1 (MQ=255)
                                                                       aaGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGctct                                                                         <  6:151280/149‑1 (MQ=255)
                                                                                    cAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACACCGCCCCCTCAGGACACAACGCCCATGCTGCTCTATACTCCCGCTCTTTTTCATTCGCAc                                                            >  8:99610/1‑149 (MQ=255)
                                                                                       ggTCCTGGCAGCCCTGAGCTGTTCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCGCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTg                                                         <  7:99610/149‑1 (MQ=255)
                                                                                            tGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATAc                                                    >  1:236532/1‑149 (MQ=255)
                                                                                                       aGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATg                                          <  8:86459/148‑1 (MQ=255)
                                                                                                             tCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCAccc                                   <  5:221803/149‑1 (MQ=255)
                                                                                                             tCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCAccc                                   <  1:212107/149‑1 (MQ=255)
                                                                                                                          gCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGa                      <  7:34029/149‑1 (MQ=255)
                                                                                                                                 gACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGAGGACCACCCTGCTGGCCCTGGATACCGcc               >  3:234166/1‑149 (MQ=255)
                                                                                                                                  aCAACGGCCCGCAAGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCa              >  6:94512/1‑149 (MQ=255)
                                                                                                                                   cAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCAc             >  1:197336/1‑149 (MQ=255)
                                                                                                                                    aaCGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCAcc            >  1:295174/1‑149 (MQ=255)
                                                                                                                                    aaCGGCCCGCAAGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCAcc            >  3:312418/1‑149 (MQ=255)
                                                                                                                                              ttGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCACCGAAGCCTGtt  >  7:199392/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
CTACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCACCGAAGCCTGTT  >  NC_002947/1712809‑1713099

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: