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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB122656,9050AC32.8% 37.4 / 6.3 45T241P (ACG→CCG) yfhMAlpha‑2‑macroglobulin
Reads supporting (aligned to +/- strand):  ref base A (23/6);  new base C (7/9);  total (30/15)
Fisher's exact test for biased strand distribution p-value = 2.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGACATGCTCTCGACCTTCCTGGCCAACTGGGGACGCAACAGCAGCCTGTACTACTACCAGTCCAAGGAAACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGTACCTACGATTATTCTCAGCCGGCAACCCTGTTCACCCTTAGCGACAT  >  ACB122/656774‑657038
                                                                                                                                   |                                                                                                                                     
cAGACATGCTCTCGACCTTCCTGGCCAACTGGGGACGCAACAGCAGCCTGTACTACTACCAGTCCAAGGAAACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCACCAcccg                                                                                                                                   >  3:302763/1‑136 (MQ=255)
                cTTCCTGGCCAACTGGGGACGCAACAGCAGCCTGTACTACTACCAGTCCAAGGAAACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACCCGCGCGAGACCGTGTtgct                                                                                                                    >  2:84955/1‑135 (MQ=255)
                     tGGCCAACTGGGGACGCAACAGCAGCCTGTACTACTACCAGTCCAAGGACACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCGGCGAGCTCAACCCGGCGCCCACCCCCCGCGAGCCCCGGTtgcgtccgc                                                                                                               >  6:261293/1‑129 (MQ=255)
                     tGGCCAACTGGGGACGCAACAGCAGCCTGTACTACTACCAGTCCAAGGAAACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACACGCCCGAGCCCGTGTTGCTGCCGAt                                                                                                              >  7:315331/1‑136 (MQ=255)
                      ggCCAACTGGGGACTCAACAGCAGCCTGTACTACTACCAGTCCAAGGAAACCCAGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACACGGCGCGCCACAACCGCCAGACCAGGTTTATGCCGaa                                                                                                              >  6:240318/1‑134 (MQ=255)
                      ggCCAACTGGGGACGCAACAGCAGCCTGTACTACTACCAGTCCAAGGAAACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCGTCGATCTCAACCCCGCCCGCACACCCCGCGAGAACCGGGTGCTGCCgca                                                                                                              >  5:328412/1‑133 (MQ=255)
                                    gCAACAGCAGCCTGTACTACTACCAGTCCAAGGAAACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCCGCATCAAGc                                                                                                >  2:308615/1‑135 (MQ=255)
                                            ccccTGTACACCTACCAGGCCACGCAACCCCTGGACATGGCCGAGCTGGTCTACCGCGGCCGCTGCGATCTCACCCCGGCGCGAAACCCGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAgg                                                                                       <  7:181916/134‑1 (MQ=255)
                                            agcCTGTACTACAACCAGTCCAAGGAACCCCTGGACAGGGCCGCGCGGGCAAACCGCGGCCGCGTCGAGCTCACCCCGGCGCGCAACCCGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAgg                                                                                       <  8:108113/136‑1 (MQ=255)
                                                tGTACAACAACCAGCCCACGGAACCCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTGCGATCTCAACCCGGCGCCCACCCCGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGgcgc                                                                                    <  6:344479/135‑1 (MQ=255)
                                                  tactactacCAGTCCACGCAAACCCGGCACATGGCCGCGCTGGGCAACAGCGGCCGCTGCGATCCCACCCCGGCGCGCAACCCGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGGTGCAGGCGCCa                                                                                  <  7:59569/135‑1 (MQ=255)
                                                     tactacCAGTCCAAGGAAACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCCACACGCGCGCGCCCGTGTTGCTGCCGAATGCCGGCATCACACCGTTGCAGGCGCCCGGc                                                                               >  5:26523/1‑135 (MQ=255)
                                                     tactacCAGTCCAAGGAAACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCCACACGCCCGAGCCCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGc                                                                               >  5:177047/1‑135 (MQ=255)
                                                     tactacCAGTCCAAGGAAACCCTGGACATCGCCGAGCTGGTCTACAGCGGCCGCATCGATCACAACCCGGCGCGCAACACGCGCGAGCCCGTGTTGCTGCCGCATGCCGGCCTCAAGCAGATGACGGCGCCAGCCg                                                                              >  8:364734/1‑136 (MQ=255)
                                                       atcccAGACCAACCAAACCCGGCAAAGGGCCGAGCTGGAATACAGCGCCCGCTCCGATCTCAACCCCGCGCGCACCACGCCCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCgt                                                                             <  2:174305/132‑1 (MQ=255)
                                                          ccAGTCCAAGGAAACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTAc                                                                         >  4:36778/1‑136 (MQ=255)
                                                           cAGTCCAAGGAAACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACCCGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCgtgtcc                                                                         >  6:98564/1‑133 (MQ=255)
                                                           cAGTCCAAGGAAACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACACGCGCGCGACCGTGTTGCTGCCGATTGCCGGCATAAAGCCGTTGCAGGCGCCCGGCGTGTAc                                                                         >  6:298934/1‑135 (MQ=255)
                                                            aGTCCAAGGAAACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCCACCCGCCCGACACCCTGTTGCTGACCATTGCCGGCCTCAACCCGTTGGAGGCGCCCAGCCTGTAcc                                                                        >  5:121216/1‑135 (MQ=255)
                                                               ccAAGGAAACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTg                                                                      >  1:27183/1‑134 (MQ=255)
                                                               ccAAGGAAACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACACCCGCGCGCACGTGTTGCTGCCGATTGCCGGCATCAAGACGTTGGACGCGCCCGGCGGGTACCTgg                                                                     >  2:150913/1‑135 (MQ=255)
                                                                 aaGGAAACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCCGGCGCCAGGCGTGGACCTGGCAg                                                                  >  7:28013/1‑136 (MQ=255)
                                                                           tGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACCCGCGCGCGACCGTGTTGCTGCCGATTGCCGGCCTCCAGCCGGTGCCGGCGCCAGGCGTGTACCTGGCAGTGATgcgcg                                                         >  6:36075/1‑135 (MQ=255)
                                                                           tGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCCTGTTTCTGCCGATTTCCGGCCTCAACCCCTTGCAGGCCCCCGGCGGGGACCGGGCAGTGCTgggcgc                                                        >  7:190752/1‑136 (MQ=255)
                                                                           tGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCCGCGCGCAACACGCGCGAGACCGTTTTGCTGCCGATTGCCGGCATCAAGCCGTTTCAGGCGCCAGGCGTGTACCTGGCAGTGATgcgcg                                                         >  1:31128/1‑135 (MQ=255)
                                                                           tGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTAAACCCGCCGCGCAACCCCCCCCAGACCGTGTTGCTGGCCAATTCCGGCAAGAACCCGTTGCATGCCCCAGGCCTTTAACAGGCACTGATtcgcg                                                         >  6:392842/1‑135 (MQ=255)
                                                                              aCAGGGCCGAGCGGGTCTACCGCGGCCGCTTCGATCTCAACCCGGCGCGCACCCCGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTc                                                     <  3:80832/136‑1 (MQ=255)
                                                                                  ggCCGAGCTGGTCTCCCGCGGCCGCTGCGCTCTCAACCCGGCGCGCACCCCGCGCGAGACCGTGTTGCTGCCGATTGCCGGCAGCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCgg                                                   <  4:280834/134‑1 (MQ=255)
                                                                                     cGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACACGCGCGCGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCCCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGTAcc                                              >  8:259995/1‑136 (MQ=255)
                                                                                        gCTGGTCTACCGCGGCCGCTTCTAGCTCACCCCGGCGCGCAACCCGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGTACCTa                                            <  5:69621/135‑1 (MQ=255)
                                                                                          tGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCCAAACGCGCGAGACCGTGTTGCTGCCGATTGCAGGCATCAACCCCTTGCAGGCGCCAGGCGTGTACCTGGCCGTCATGCGCGCCTCGTGTTCCTACGa                                         >  8:216114/1‑136 (MQ=255)
                                                                                          tGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACACGCGCGCGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTAGCTGGCAGTGCTGCGCGCCTCGGGTACCTACGa                                         >  3:76008/1‑136 (MQ=255)
                                                                                                accgCGGCCGCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGTACCTACGATTATtc                                   <  4:302763/133‑1 (MQ=255)
                                                                                                aCAGCGGCCGCTTTGATCGAACCCCGGCGCGCAACCCGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGTACCTACGATTATtc                                   <  8:189970/136‑1 (MQ=255)
                                                                                                aCAGCGGCCGCTTCGATCTCAACCCGGCGCGCCACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGGACCTGGCAGTGCTGCGCGCCTCGGGTACCTACGattatt                                    >  6:196595/1‑135 (MQ=255)
                                                                                                 ccgcgcCCGCTGCCACCCCACCCCGGCGCCCAACCCGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGTACCTACGATTATtc                                   <  5:398441/133‑1 (MQ=255)
                                                                                                       ccGCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGGTGCAGGCGCCCGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGTACCTACGATTATTCTCAGCCg                            >  4:158051/1‑136 (MQ=255)
                                                                                                         gCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGCCCGTGTTGCTGCCGATTGCCGGCATCAAGCCCTTGCACGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGTACCATCGATTATTCTCAGCCGGc                          >  3:344082/1‑136 (MQ=255)
                                                                                                         gCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGTACCTACGATTATTCTCAGCCgg                           >  5:153078/1‑135 (MQ=255)
                                                                                                         gCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGTACCTACGATTATTCTCAGCCGGc                          <  4:355500/136‑1 (MQ=255)
                                                                                                         gCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCCTGTTGCTGCCGATTGCCGGCAACACGCCCGTGCAGGCGCCAGGCGTGTACCTGGCACTGCTTCGCGCCTCGGCTACCTAAGATTATTCTCAGCCGGc                          >  4:159855/1‑136 (MQ=255)
                                                                                                         gCTACGATCTCACCCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGTACCTACGATTATTCTCAGCCGGc                          <  4:344082/136‑1 (MQ=255)
                                                                                                           ttggATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGTACCTACGATTATTCTCAGCCGGCa                         <  6:177047/132‑1 (MQ=255)
                                                                                                                        cGGCGCCCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGTACCTACGATTATTCTCAGCCGGCAACCCTGTTCACCCt           <  4:325328/136‑1 (MQ=255)
                                                                                                                            gcgcAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGTACCTACGATTATTCTCAGCCGGCAACCCTGTTCACCCTTAg        >  6:110775/1‑135 (MQ=255)
                                                                                                                                 acacGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCCGGCGTGTACCGGGCAGTGCTGCGCGCCGCGGGGACCTACGATTATTCTCAGCCGGCCACCCTGTTCCCCCTTAGCGACAt  >  7:273202/1‑136 (MQ=255)
                                                                                                                                   |                                                                                                                                     
CAGACATGCTCTCGACCTTCCTGGCCAACTGGGGACGCAACAGCAGCCTGTACTACTACCAGTCCAAGGAAACCCTGGACATGGCCGAGCTGGTCTACAGCGGCCGCTTCGATCTCAACCCGGCGCGCAACACGCGCGAGACCGTGTTGCTGCCGATTGCCGGCATCAAGCCGTTGCAGGCGCCAGGCGTGTACCTGGCAGTGATGCGCGCCTCGGGTACCTACGATTATTCTCAGCCGGCAACCCTGTTCACCCTTAGCGACAT  >  ACB122/656774‑657038

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: ‑