Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,644,703 T→C 11.1% E13E (GAA→GAG)  lon‑I ← DNA‑binding ATP‑dependent protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,644,7030TC11.1% 89.1 / 3.2 36E13E (GAA→GAG) lon‑IDNA‑binding ATP‑dependent protease
Reads supporting (aligned to +/- strand):  ref base T (18/14);  new base C (2/2);  total (20/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

CATTGACGATCACCGGCAGTACCTGGGCCGGGAAGAACGGGCGGTTGTGGATCGGGATGACGTAGACCTTGTCCGGCAGCTGTTGACCGGGCAGGGCGAGGGCGTGGCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGTTCGTCAGGGTGTTCGGGGAAATCCTGCTGGTCGCTCATGGGGCACCTGCGTGAATGACTATGGTGCTTAGATGGGGCGGGGAGGGATGGGTTTCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAA  >  NC_002947/1644558‑1644850
                                                                                                                                                 |                                                                                                                                                   
cATTGACGATCACCGGCAGTACCTGGGCCGGGAAGAACGGGCGGTTGTGGATCGGGATGACGTAGACCTTGTCCGGCAGCTGTAGACCGGGCAGGGCGAGGGCGTGGCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTgtgt                                                                                                                                                     <  4:247795/146‑1 (MQ=255)
 aTTGACGATCACCGGCAGTACCTGGGCCGGGAAGAACGGGCGGTTGTGGATCGGGATGACGTAGACCTTGTCCGGCAGCTGTTGACCGGGCAGGGCGAGGGCGTGGCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGTTCGt                                                                                                                                                 >  6:316436/1‑149 (MQ=255)
 aTTGACGATCACCGGCAGTACCTGGGCCGGGAAGAACGGGCGGTTGTGGATCGGGATGACGTAGACCTTGTCCGGCAGCTGTTGACCGGGCAGGGCGAGGGCGTGGCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGTTCGt                                                                                                                                                 <  5:62350/149‑1 (MQ=255)
      cGATCACCGGCAGTACCTGGGCCGGGAAGAACGGGCGGTTGTGGATCGGGATGACGTAGACCTTGTCCGGCAGCTGTTGACCGGGCAGGGCGAGGGCGTGGCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGTTCGTCAgg                                                                                                                                             <  6:17449/148‑1 (MQ=255)
         tCACCGGCAGTACCTGGGCCGGGAAGAACGGGCGGTTGTGGATCGGGATGACGTAGACCTTGTCCGGCAGCTGTTGACCGGGCAGGGCGAGGGCGTGGCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGTTCGTCAGGgtgt                                                                                                                                         >  4:57416/1‑149 (MQ=255)
              ggCAGTACCTGGGCCGGGAAGAACGGGCGGTTGTGGATCGGGATGACGTAGACCTTGTCCGGCAGCTGTTGACCGGGCAGGGCGAGGGCGTGGCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGTTCGTCAGGGTGTTCggg                                                                                                                                    <  4:238887/149‑1 (MQ=255)
               gCAGTACCTGGGCCGGGAAGAACGGGCGGTTGTGGATCGGGATGACGTAGACCTTGTCCGGCAGCTGTTGACCGGGCAGGGCGAGGGCGTGGCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGTTCGTCAGGGTGTTCggg                                                                                                                                    >  4:213310/1‑148 (MQ=255)
                  gTACCTGGGCCGGGAAGAACGGGCGGTTGTGGATCGGGATGACGTAGACCTTGTCCGGCAGCTGTTGACCGGGCAGGGCGAGGGCGTGGCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGTTCGTCAGGGTGTTCGGGGaaa                                                                                                                                <  6:4929/149‑1 (MQ=255)
                     ccTGGGCCGGGAAGAACGGGCGGTTGTGGATCGGGATGACGTAGACCTTGTCCGGCAGCTGTTGACCGGGCAGGGCGAGGGCGTGGCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGTTCGTCAGGGTGTTCgg                                                                                                                                     <  4:273264/141‑1 (MQ=255)
                     ccTGGGCCGGGAAGAACGGGCGGTTGTGGATCGGGATGACGTAGACCTTGTCCGGCAGCTGTTGACCGGGCAGGGCGAGGGCGTGGCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGTTCGTCAGGGTGTTCgg                                                                                                                                     >  3:273264/1‑141 (MQ=255)
                      cTGGGCCGGGAAGAACGGGCGGTTGTGGATCGGGATGACGTAGACCTTGTCCGGCAGCTGTTGACCGGGCAGGGCGAGGGCGTGGCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGCTCGTCAGGGTGTTCGGGGAAATCct                                                                                                                            >  1:238524/1‑149 (MQ=255)
                      cTGGGCCGGGAAGAACGGGCGGTTGTGGATCGGGATGACGTAGACCTTGTCCGGCAGCTGTTGACCGGGCAGGGCGAGGGCGTGGCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGCTCGTCAGGGTGTTCGGGGAAATCct                                                                                                                            >  1:234081/1‑149 (MQ=255)
                         ggCGGGGAAGCACGGGCGGTTGTGGATCGGGATGACGTAGACCTTGTCCGGCAGCTGTTGACCGGGCAGGGCGAGGGCGTGGCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGTTCGTCAGGGTGTTCGGGGAAATCctgc                                                                                                                          <  2:76896/148‑1 (MQ=255)
                              ggAAGAACGGGCGGTTGTGGATCGGGATGACGTAGACCTTGTCCGGCAGCTGTTGACCGGGCAGGGCGAGGGCGTGGCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGTTCGTCAGGGTGTTCGGGGAAATCCTGCTGGTCg                                                                                                                    >  3:213463/1‑149 (MQ=255)
                                  gaaCGGGCGGTTGTGGATCGGGATGACGTAGACCTTGTCCGGCAGCTGTTGACCGGGCAGGGCGAGGGCGTGGCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGTTCGTCAGGGTGTTCGGGGAAATcc                                                                                                                             <  3:222949/136‑1 (MQ=255)
                                  gaaCGGGCGGTTGTGGATCGGGATGACGTAGACCTTGTCCGGCAGCTGTTGACCGGGCAGGGCGAGGGCGTGGCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGTTCGTCAGGGTGTTCGGGGAAATcc                                                                                                                             >  4:222949/1‑136 (MQ=255)
                                     cGGGCGGTTGTGGATCGGGATGACGTAGACCTTGTCCGGCAGCTGTTGACCGGGCAGGGCGAGGGCGTGGCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGCTCGTCAGGGTGTTCGGGGAAATCCTGCTGGTCGCTCATgg                                                                                                             <  2:238524/149‑1 (MQ=255)
                                     cGGGCGGTTGTGGATCGGGATGACGTAGACCTTGTCCGGCAGCTGTTGACCGGGCAGGGCGAGGGCGTGGCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGCTCGTCAGGGTGTTCGGGGAAATCCTGCTGGTCGCTCATgg                                                                                                             <  2:234081/149‑1 (MQ=255)
                                       ggCGGTTGTGGATCGGGATGACGTAGACCTTGTCCGGCAGCTGTTGACCGGGCAGGGCGAGGGCGTGGCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGTTCGTCAGGGTGTTCGGGGAAATCCTGCTGGTCGCTCATgggg                                                                                                           <  5:316436/149‑1 (MQ=255)
                                          ggTTGTGGATCGGGATGACGTAGACCTTGTCCGGCAGCTGTTGACCGGGCAGGGCGAGGGCGTGGCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGTTCGTCAGGGTGTTCGGGGAAATCCTGCTGGTCGCTCATGGGGCAc                                                                                                        <  1:181001/149‑1 (MQ=255)
                                              gtgGATCGGGATGACGTAGACCTTGTCCGGCAGCTGTTGACCGGGCAGGGCGAGGGCGTGGCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGTTCGTCAGGGTGTTCGGGGAAATCCTGCTGGTCGCTCATGGGGCACCTGc                                                                                                    >  4:20995/1‑149 (MQ=255)
                                                ggATCGGGATGACGTAGACCTTGTCCGGCAGCTGTTGACCGGGCAGGGCGAGGGCGTGGCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGTTCGTCAGGGTGTTCGGGGAAATCCTGCTGGTCGCTCATGGGGCACCTGCGt                                                                                                  >  7:27836/1‑149 (MQ=255)
                                                                   cTTGTCCGGCAGCTGTTGACCGGGCAGGGCGAGGGCGTGGCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGTTCGTCAGGGTGTTCGGGGAAATCCTGCTGGTCGCTCATGGGGCACCTGCGTGAATGACTATGGTGCt                                                                                  >  3:154744/1‑146 (MQ=255)
                                                                   cTTGTCCGGCAGCTGTTGACCGGGCAGGGCGAGGGCGTGGCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGTTCGTCAGGGTGTTCGGGGAAATCCTGCTGGTCGCTCATGGGGCACCTGCGTGAATGACTATGGTGCt                                                                                  >  2:59590/1‑146 (MQ=255)
                                                                            cAGCTGTTGACCGGGCAGGGCGAGGGCGTGGCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGTTCGTCAGGGTGTTCGGGGAAATCCTGCTGGTCGCTCATGGGGCACCTGCGTGAATGACTATGGTGCTTAGATGGGGCgg                                                                      <  3:57416/149‑1 (MQ=255)
                                                                             aGCTGTTGACCGGGCAGGGCGAGGGCGTGGCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGTTCGTCAGGGTGTTCGGGGAAATCCTGCTGGTCGCTCATGGGGCACCTGCGTGAATGACTATGGTGCTTAGATGGGGCggg                                                                     <  6:118764/149‑1 (MQ=255)
                                                                              gCTGTTGACCGGGCAGGGCGAGGGCGTGGCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGTTCGTCAGGGTGTTCGGGGAAATCCTGCTGGTCGCTCATGGGGCACCTGCGTGAATGACTATGGTGCTTAGATGGGGCgggg                                                                    <  3:20995/149‑1 (MQ=255)
                                                                                       gggcgagggcgagggcgTGGCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGTTCGTCAGGGTGTTCGGGGAAATCCTGCTGGTCGCTCATGGGGCACCTGCGTGAATGACTATGGTGCTTAGATGGGGCGGGGAGGGATggg                                                           >  3:283156/6‑149 (MQ=255)
                                                                                                     gcgTGGCTGACTGGGGCCGGGAGGGTGTGTTCGCCTTCGGTGTGTCCGGCAGGGGGTCCGGGCAAACCCTGCGGGGCGCTCATGGGCCACCTGCGTGAATGACTATGGGGCTTCGAGGGGGCGGGGAGGGATGGGTTTCAATGGTTCgg                                             <  4:213463/149‑1 (MQ=255)
                                                                                                         ggCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGTTCGTCAGGGTGTTCGGGGAAATCCTGATGGTCGCTCATGGGGCACCTGCGTGAATGACTATGGTGCTTAGATGGGGCGGGGAGGGATGGGTTTCAATGGTTCGGGagc                                         >  5:248744/1‑149 (MQ=255)
                                                                                                           cTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGTTCGTCAGGGTGTTCGGGGAAATCCTGCTGGTCGCTCATGGGGAACATGAGTGAATGACTATGGTGCTTAGATGGGGCTGGGAGGGATGGGTTTCAATGGTTCGGTagcag                                       >  8:45689/1‑149 (MQ=255)
                                                                                                            tGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGTTCGTCAGGGTGTTCGGGGAAATCCTGCTGGTCGCTCATGGGGCACCTGCGTGAATGACTATGGTGCTTAGATGGGGCGGGGAGGGATGGGTTTCAATGGTTCTGGagcagc                                      >  1:158859/1‑149 (MQ=255)
                                                                                                               ctgcGGCCGGGCGGTGGTGTCCTACTCGGCTGTGTCCGTCCGGGTGGCCGGTTCAAGCCGGCTGGGCCCTCATGGGGCACCTGCGTGAATTACTATGGGGCTTAGATGGGGCGGGGCGGGAGGGGTTTCAATGGTTCGGGAGCAGCCtt                                   <  8:232014/146‑1 (MQ=255)
                                                                                                                ttCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGTTCGTCAGGGTGTTCGGGGAAATCCTGCTGGTCGCTCATGGGGCACCTGCGTGAATGACTATGGTGCTTAGATGGGGCAGGGAGGGATGGGTTTCAATGGTTCGGGAGCAGCCTtg                                  >  1:7276/1‑149 (MQ=255)
                                                                                                                      cTGGAGGGTGTGTTCGACTTCGCTGTGTTCGTCAGGGTGTTCGGGGAAATCCTGCTGGTCGCTCATGGGTCACCTGCGTGAATGACTATGGTGCTTAGATGGGGCGGGGAGGGATGGGTTTCAATGGTTCGGGAGCAGCCTTGTGTc                              >  4:22484/1‑147 (MQ=255)
                                                                                                                                  ttCGACTTCGCTGTGTTCGTCAGGGTGTTCGGGGAAATCCTGCTGGTCGCTCATGGGGCACCTGCGTGAATGACTATGGTGCTTAGATGGGGCGGGGAGGGATGGGTTTCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTgcg                >  5:231108/1‑149 (MQ=255)
                                                                                                                                        ttCGCTGTGTTCGTCAGGGTGTTCGGGGAAATCCTGCTGGTCGCTCATGGGGCACCTGCGTGAATGACTATGGTGCTTAGATGGGGCGGGGAGGGATGGGTTTCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGcagcag          >  5:82946/1‑149 (MQ=255)
                                                                                                                                           gcggtgtTCGTCAGGGCGTTGGGGTAATTCCTGCTGGTCGGTCTTGGGGCACCTGCGTTAATGACTTTGGTGCTTCGATGGGGGGGGCAGGGAGGGGTTTCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGccc       <  3:22484/146‑1 (MQ=255)
                                                                                                                                             tgtgtgCGGCAGGTGGTTCGGGGAAATCCTGGTGGGCGCTCATGGGTCACCTGCGTGAATGACTAGGGTGCTTAGATGGGGCGGGTTGGGATGGGTTTCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCa     <  4:283156/149‑1 (MQ=38)
                                                                                                                                                gtTCGTCAGGGTGTTCGGGGAAATCCTGCTGGTCGCTCATGGGGCACCTGCGTGAATGACTATGGTGCTTAGATGGGGCGGGGAGGGATGGGTTTCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGaa  >  2:147326/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CATTGACGATCACCGGCAGTACCTGGGCCGGGAAGAACGGGCGGTTGTGGATCGGGATGACGTAGACCTTGTCCGGCAGCTGTTGACCGGGCAGGGCGAGGGCGTGGCTGACTTCGGCCTGGAGGGTGTGTTCGACTTCGCTGTGTTCGTCAGGGTGTTCGGGGAAATCCTGCTGGTCGCTCATGGGGCACCTGCGTGAATGACTATGGTGCTTAGATGGGGCGGGGAGGGATGGGTTTCAATGGTTCGGGAGCAGCCTTGTGTCGCGATTGGGCTGCGCAGCAGCCCCAGAA  >  NC_002947/1644558‑1644850

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