Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 656,905 A→C 25.6% T241P (ACG→CCG)  yfhM → Alpha‑2‑macroglobulin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT2440656,9050AC25.6% 48.0 / 3.1 49T241P (ACG→CCG) yfhMAlpha‑2‑macroglobulin
Reads supporting (aligned to +/- strand):  ref base A (25/10);  new base C (5/9);  total (30/19)
Fisher's exact test for biased strand distribution p-value = 2.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGACATGCTCTCGACCTTCCTGGCCAACTGGGGACGCAACAGCAGCCTGTACTACTACCAGTCCAAGGAAACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGTACCTACGATTATTCTCAGCCGGCAACCCTGTTCACCCTTAGCG  >  KT2440/656773‑657034
                                                                                                                                    |                                                                                                                                 
ccAGACATGCTCTCGACCTTCCTGGCCAACTGGGGACGCAACAGCAGCCTGTACTACTACCAGTCCAAGGAACCCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACAcg                                                                                                                                 >  7:12158/1‑135 (MQ=255)
 cAGACATGCTCTCGACCTTCCTGGCCAACTGGGGACGCAACAGCAGCCTGTACTACTACCAGTCCAAGGAAACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACAcgc                                                                                                                                >  8:192623/1‑135 (MQ=255)
            tcGACCTTCCTGGCCAACTGGGGACGCAACAGCAGCCTGTACTACTACCAGTCCAAGGAAACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCgtg                                                                                                                     >  5:334637/1‑135 (MQ=255)
            tcGACCTTCCTGGCCAACTGGGGACGCAACAGCAGCCTGTACTACTACCAGTCCAAGGAAACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCgtg                                                                                                                     >  4:156551/1‑135 (MQ=255)
                      tGGCCAACTGGGGACGCAACAGCAGCCTGTACTACTACCAGTCCAAGGAAACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCACCACCCCCGAGACCCTGATGCTGCCGa                                                                                                           >  8:131494/1‑135 (MQ=255)
                       ggCCAACTGGGGACGCAACAGCAGCCTGTACTACTACCAGTCCAAGGAAACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCCTGTAGCTGCCGAt                                                                                                          >  5:148987/1‑135 (MQ=255)
                         ccAACTGGGGACGCACCAGCACCCGGTACTACTACCACTCCACGGACACCCTCGACATGGCCGTGCTGGACTACAGCGGCCGCTTCGAGCGCAACCCGGCGCGCAACCCGCGCGAGACCGTGTTGCTGCCGAt                                                                                                          <  2:262582/133‑1 (MQ=255)
                              tGGGGACGCAACAGCAGCCTGTACTCCTCCCAGGCCAAGGAAACCCTGGACAGGGCCGAGCTGGGCTACAGCGGCCGCTTCGATCGCAACCCGGCGCCCAACCCGCGCGAGACCGTGTTGCTGCCGATTGCCGGc                                                                                                   <  4:309039/135‑1 (MQ=255)
                                gggACGCAACAGCAGCCTGTACTACTACCAGTCCAAGGAAACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACCCGCGCGAGAACCGGTTTCTGCCGATTTCCGGCa                                                                                                  >  3:212119/1‑134 (MQ=255)
                                              gcCTGTCCTACTACCAGCCCAACGAAACCTTGGACATGGCCGAGCGGGTGAACAGCGGCCGCTTCCCTGTAAACCCGGGGCGCAACCCGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAgg                                                                                   <  8:12158/135‑1 (MQ=255)
                                              gcCTGTAATACTACCAGACCAAGTAACCCCTGGACATGGCCTAGCTTGTCTACCGGGGCCGGTTCTATCCCACCCCGCCGCCCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAgg                                                                                   <  4:212119/135‑1 (MQ=255)
                                                     ctactccCAGCCCAAGGAAACCCTGGACATGGCCGAGCGGGGCTACAGCGGCCGCTTCGAGCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAgg                                                                            <  7:192623/135‑1 (MQ=255)
                                                      tactacCAGTCCAAGGAAACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGc                                                                           >  3:63620/1‑135 (MQ=255)
                                                            cAGTCCAAGGAAACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCCACACCCGCCCGACCGTGGTGCTTCCGATTTCCGGGATCAAGCCCTTTCCCGCGCCCGGCGTGTAc                                                                     >  8:288228/1‑135 (MQ=255)
                                                                ccAAGGAACCCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGCTCTCAACCCGGCGCGCAACCCGCCCGCGCCCGTGTTTCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCACGCGGGTTCCTgg                                                                 >  8:128141/1‑135 (MQ=255)
                                                                ccAAGGAAACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGGACCTgg                                                                 >  6:109322/1‑135 (MQ=255)
                                                                 cAAGTAACCCGTGCAAATGGCCCCGCTGGTCAACCGCGGCCGCTCCGATCCCACCCCGGCGCCAACCCCGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGc                                                                <  6:289030/135‑1 (MQ=255)
                                                                 cAAGAAACCCCTGCAATTGGCCGCGCGGGGCAACCGGGGCCGCTTCGATCACACCCGGGCGCGCACCCCGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGc                                                                <  6:283997/135‑1 (MQ=255)
                                                                  aaGGAAACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACCCGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCa                                                               >  2:266248/1‑135 (MQ=255)
                                                                  aaGGAAACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCa                                                               >  3:153343/1‑135 (MQ=255)
                                                                     gAAACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACCCGCGCGGGCCCGTGTTGCTGCCGATTGCCGGCATCACGCCGGTGCAGGCGCCCGGCGGGGACCTGGCAGTg                                                            >  2:241976/1‑135 (MQ=255)
                                                                      aaaCCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCCGGCGTGTACCTGGCAGTgc                                                           >  3:83647/1‑134 (MQ=255)
                                                                       accccTGGCCTGGGCCGCGCGCGTCAACTGGGGCGGCGTCGCTCGCCCCCGGGCGCCCACCAGCCCCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGAt                                                          <  2:201769/133‑1 (MQ=255)
                                                                            tGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGATgcgcg                                                     >  8:204114/1‑135 (MQ=255)
                                                                                 aTGGCCGAGCGGGTCAACCGCGGCCGCTTCGATCTCACCCCGGCGCGCAACCCGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCg                                                <  7:110832/135‑1 (MQ=255)
                                                                                      cGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGTAc                                           <  6:146506/135‑1 (MQ=255)
                                                                                      cGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGTAc                                           >  7:187670/1‑135 (MQ=255)
                                                                                       gAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCCGGCGGGGTGCTGGCAGTGCTGCGCGCCGCGGGTGcc                                          >  5:202620/1‑135 (MQ=255)
                                                                                           tGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGTACCTa                                        >  8:219473/1‑133 (MQ=255)
                                                                                                tACCGCGGCCGCTTCGACCTCACCCCGGCGCGCACCCCGCCCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGTACCTACGattat                                 <  5:145055/135‑1 (MQ=255)
                                                                                                tACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGTACCTACGattat                                 >  3:275200/1‑135 (MQ=255)
                                                                                                 aCAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCCGGCGTGTTCCTGGCAGTGATGCGCGCCTCGGGTACCTACGattatt                                >  8:206202/1‑135 (MQ=255)
                                                                                                   cgcgcCCGCTACGACCGCCCCCCGCCCCCCAACCCGCGCGAGACCGGGTTGCTGCCGATTGCCGGCATCAAGCCGTGGCAGGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGTACCTACGattatt                                <  5:98401/132‑1 (MQ=255)
                                                                                                     cGGCCGCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGTACCTACGATTATTCTCa                            >  1:61057/1‑135 (MQ=255)
                                                                                                       gCCGCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGTACCTACGATTATTCTCAGc                          <  6:334637/135‑1 (MQ=255)
                                                                                                        ccGCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCCGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGGACCTACGATTATTCTCAGcc                         >  4:137357/1‑135 (MQ=255)
                                                                                                          gCTTCGATCTCAACCCGGCGCGCAACCCGCGCGAGCCCGTGTTGCTGCCGATTGCCGGCATCAAGCCCTTGCAGGCGCCAGGCGTGTACCTGGCAGTGCTGCGCGCCTCGGGTACCTACGATTATTCTCAGCCgg                       >  3:75838/1‑135 (MQ=255)
                                                                                                          gCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAgg                                                                                   <  7:279489/75‑1 (MQ=255)
                                                                                                          gCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAgg                                                                                   >  8:279489/1‑75 (MQ=255)
                                                                                                          gCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCgt                                                                         >  6:77790/1‑85 (MQ=255)
                                                                                                          gCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCgt                                                                         <  5:77790/85‑1 (MQ=255)
                                                                                                          gCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGTACCTACGATTATTCTCAGCCg                        >  7:262403/1‑134 (MQ=255)
                                                                                                           cTTCGATCTCAACCCGGCGCGCAACCCGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGTACCTACGATTATTCTCAGCCGGc                      <  4:63620/135‑1 (MQ=255)
                                                                                                                      aCCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGTACCTACGATTATTCTCAGCCGGCAACCCt                <  4:275200/130‑1 (MQ=255)
                                                                                                                      aCCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGTACCTACGATTATTCTCAGCCGGCAACCCt                <  4:153343/130‑1 (MQ=255)
                                                                                                                             gcgcAACACGCGCGAGCCCGTGTTGCTGCCGATTGCCGGCATCAATCCCTTGCAGGCGCCAGGCGGGGACCTGGCAGTGATTCGCGCCTCGGGTTCCTACGATTATTCTGAGCCGGCAACCCTGGTCTCCCTTAg    >  1:78355/1‑135 (MQ=255)
                                                                                                                               gcAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAgg                                                                                   <  7:154012/54‑1 (MQ=255)
                                                                                                                               gcAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAgg                                                                                   >  8:154012/1‑54 (MQ=255)
                                                                                                                               gcAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGTACCTACGATTATTCTCAGCCGGCAACCCTGTTCACCCTTAGCg  >  6:151647/1‑135 (MQ=255)
                                                                                                                                    |                                                                                                                                 
CCAGACATGCTCTCGACCTTCCTGGCCAACTGGGGACGCAACAGCAGCCTGTACTACTACCAGTCCAAGGAAACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGTACCTACGATTATTCTCAGCCGGCAACCCTGTTCACCCTTAGCG  >  KT2440/656773‑657034

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