Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,438,375 G→A 8.0% L228L (TTG→TTA)  thiD → Hydroxymethylpyrimidine/phosphomethylpyrimidineki nase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,438,3750GA 8.0% 133.3 / 4.2 51L228L (TTG→TTA) thiDHydroxymethylpyrimidine/phosphomethylpyrimidineki nase
Reads supporting (aligned to +/- strand):  ref base G (20/26);  new base A (2/2);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.56e-01

CTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGCCCCTGGACTTCTGCTCCTGATACGAGGCCTTCA  >  ACB122/5438253‑5438502
                                                                                                                          |                                                                                                                               
cTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTgc                                                                                                                    <  8:386859/136‑1 (MQ=255)
cTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCGCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTgc                                                                                                                    <  7:394380/136‑1 (MQ=255)
cTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTA‑CATGGCGCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCGCGGCGAGCAGTTGGAAAGCGCCGTgcg                                                                                                                   <  7:54833/136‑1 (MQ=255)
 tttACACCCAGGACGGCCATCACTTCAATTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCTGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCTAGAAGTTTGAAAGCGCGGTgc                                                                                                                    <  2:287837/133‑1 (MQ=255)
 tgaacacCCCGGACGGCCACCACTACACTTGGACTGGCCCGCGCCGGCCGGGCAGCTACCATGGCGCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCGGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTgc                                                                                                                    <  5:51229/132‑1 (MQ=255)
 tGTACCCCCAGGCCGGCCAGCACTTCACTTGGACTGGCCAGCGCCTGCGGGGCAGCTACCATGGCTCGGGCTGCACCCTGCCCTGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAACGCGCCGTgc                                                                                                                    <  2:357714/135‑1 (MQ=255)
 tGTACACCCCGGACGCCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTGGGCGAGCAGTTGGAAAGCGCCGTgc                                                                                                                    <  2:320534/135‑1 (MQ=255)
 tGTACACCCAGTACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGGGGCACCCTGGCAAGCGCCCTGGCGGGCCGCATGGCGCTGGGCGAGAATTTGGAAAGCGCGGTgc                                                                                                                    <  5:295874/135‑1 (MQ=255)
 tGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTgc                                                                                                                    <  3:107697/135‑1 (MQ=255)
  gTACACCCAGGCCGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCCCACTGGCCCGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGACAGCGCCGTgcg                                                                                                                   >  5:377071/1‑135 (MQ=255)
      acCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCAcc                                                                                                               <  1:330485/135‑1 (MQ=255)
             aCGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGa                                                                                                       >  7:53207/1‑136 (MQ=255)
             aCGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGa                                                                                                       >  8:376009/1‑136 (MQ=255)
             aCGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTAGAAAGCGCCGTGCGCACCGCCCTgg                                                                                                        >  3:210196/1‑135 (MQ=255)
             aCGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTAGAAAGCGCCGTGCGCACCGCCCTgg                                                                                                        >  3:205778/1‑135 (MQ=255)
                 ccAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTa                                                                                                    <  2:22210/135‑1 (MQ=255)
                 ccAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTa                                                                                                    <  2:46982/135‑1 (MQ=255)
                    gCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTacac                                                                                                 >  2:70343/1‑135 (MQ=255)
                            cttGGACTAGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGgcgc                                                                                          <  2:248529/134‑1 (MQ=255)
                             ttGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGCCTACACCTGGCGCAc                                                                                        >  6:60465/1‑135 (MQ=255)
                             ttGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCAc                                                                                        >  6:65187/1‑135 (MQ=255)
                                  cTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGc                                                                                   >  5:288186/1‑135 (MQ=255)
                                     gCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCATTTGGAAAGCGCCGTGCGCACCGCCCAGGACTAGACCTGGCGCACCCTCCGTg                                                                                >  5:102915/1‑135 (MQ=255)
                                     gCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTgc                                                                               >  7:124780/1‑135 (MQ=255)
                                      cccgcgcCTGCCGGGCAGCTACCATGGCTCGGGCGGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTg                                                                                <  6:318605/131‑1 (MQ=255)
                                        aGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACg                                                                             >  3:195084/1‑135 (MQ=255)
                                        aGCGACTGCCGGGCAGCTACCATGGCTCGGTCTGCACCCTGTACAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGCGCAGTTGGTAAGCGCCGTGCGCACCGCCCTGGACTACACCTGTCGCGCCCAGCGTGACg                                                                             >  6:27288/1‑135 (MQ=255)
                                           gccggccgGCCAGCTACCTTGGCGCGGGCTGCACCCGGGCCAGCGCGCGGGCCGGCCGCCTGGCGCTCGGCGAGCTTTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCg                                                                          <  1:28648/135‑1 (MQ=255)
                                             cTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAg                                                                        <  3:181923/135‑1 (MQ=255)
                                             cTGCCGGGCAGCTACCATGGCGCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGc                                                                       <  8:306534/136‑1 (MQ=255)
                                              ggCCGGGCAGCTACCATGGGTCGGGCGGCACCCTGCCCAGCGCCGTGGCCGGCGGCCTGGCGCTCGGCGTGCAGTTGGAAAGCGCCGGGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGc                                                                       <  5:214587/134‑1 (MQ=255)
                                                  gggCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACt                                                                   >  1:251842/1‑135 (MQ=255)
                                                   ggCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCCCCCTGCGTGACGCCGAGCAACTg                                                                  >  3:140140/1‑135 (MQ=255)
                                                   ggCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTgg                                                                 <  7:274246/136‑1 (MQ=255)
                                                   ggCAGCTACCATGGCGCGGGCTGCACCCTGGCCAGCGCCCTGGCGGGCGGCCTGGCGCTCGGGGAGAAGTTGGAACGGGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTgg                                                                 <  7:376009/136‑1 (MQ=255)
                                                   ggCAGCGACCCTGGCGCGGGCTGCCCCCGGGCCCGCGCCCGGGCCGGCCGCCTGGCGCGCGGCGAGCTGTTGGAAAGCGCCGGGCGCACCGCCCTGGACTACACCGGGCGCACCCTGCGTGACGCCGAGCAACTg                                                                  <  6:377071/135‑1 (MQ=255)
                                                    gCAGCTACCATTGCTCGGGGTCCCCGGGGCCCAGCGCCCTGGCCGGGGGCCTGGCGGTCGGCGTGCTGTTGGAAAGCGCCGTGCGCACCGCGCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTgg                                                                 <  3:185988/135‑1 (MQ=255)
                                                    gCAGCTACCATGGCTCGGGCGGCACCCTGGCCAGCGCCCGGGCCGGCCGCCTGGCGCTCGGCTAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTgg                                                                 <  3:11211/135‑1 (MQ=255)
                                                           ccATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGc                                                         >  8:57704/1‑136 (MQ=255)
                                                             aTGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGcc                                                        <  6:288186/135‑1 (MQ=255)
                                                                     ggCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTg                                                >  4:348528/1‑135 (MQ=255)
                                                                             ccTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTgccgcgccg                                        >  5:94839/1‑135 (MQ=255)
                                                                                   cAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGcc                                  >  2:12761/1‑135 (MQ=255)
                                                                                         ccTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGccccgtga                            >  3:173728/1‑131 (MQ=255)
                                                                                                ggccgCCTGGCGCGCGGCGAGCAGTTGGAAAGCGCCGGGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGCCCCTGGACTTCTGc                     <  5:95141/135‑1 (MQ=255)
                                                                                                     ccTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGcc                                                                           <  2:354406/76‑1 (MQ=255)
                                                                                                     ccTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGcc                                                                           >  1:354406/1‑76 (MQ=255)
                                                                                                       tGGCGCTCGGCGAGCAGTTAGAACGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGAGCCAATACGTGCCGCGCCGCCTGCCCCTGGACTTCTGCTCCTGAt              <  4:205778/135‑1 (MQ=255)
                                                                                                       tGGCGCTCGGCGAGCAGTTAGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGAGCCAATACGTGCCGCGCCGCCTGCCCCTGGACTTCTGCTCCTGAt              <  4:210196/135‑1 (MQ=255)
                                                                                                                  gAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGCCCCTGGACTTCTGCTCCTGATACGAGGCCtt    >  5:204430/1‑134 (MQ=255)
                                                                                                                   aGCTGTGGGAACGCGCCGGGCGCACCGCCCTGGACAACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAAAACGTGCCGCGCCGCCTGCCCCTGGACTTCTGCTCCTGATACGAGGCCTTCa  <  3:228029/135‑1 (MQ=255)
                                                                                                                          |                                                                                                                               
CTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGCCCCTGGACTTCTGCTCCTGATACGAGGCCTTCA  >  ACB122/5438253‑5438502

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