Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 5,693,153 T→G 33.6% R228R (CGT→CGG pyrC' → Dihydroorotase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1435,693,1530TG33.6% 26.7 / 2.7 39R228R (CGT→CGGpyrC'Dihydroorotase‑like protein
Reads supporting (aligned to +/- strand):  ref base T (9/16);  new base G (4/10);  total (13/26)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCGACCTGGCCCAGGGCGGCCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCC  >  ACB143/5693018‑5693288
                                                                                                                                        |                                                                                                                                       
ccGCGACCTGGCCCAGGGCGGCCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGgcgtgcgt                                                                                                                                         <  6:175872/136‑1 (MQ=255)
     aCCTGGCCCAGGGCGGCCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGtgggggcggc                                                                                                                                    >  5:203227/1‑134 (MQ=255)
      ccTGGCCCAGGGCGGCCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTgcgccc                                                                                                                                   >  1:16931/1‑134 (MQ=255)
      ccTGGCCCAGGGCGGCCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCCTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGtgcgtgcg                                                                                                                                      >  5:253854/1‑133 (MQ=255)
          gCCCAGGGCGGCCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGATACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTc                                                                                                                                <  1:3480/135‑1 (MQ=255)
          gCCCAGGGCGGCCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTc                                                                                                                                <  2:65331/135‑1 (MQ=255)
            ccAGGGCGGCCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGCGCCCTTGAGc                                                                                                                             >  4:101217/1‑136 (MQ=255)
            ccAGGGCGGCCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGCGCCCGTCAGc                                                                                                                             >  4:201062/1‑136 (MQ=255)
                  cggcCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCCGCCAGAt                                                                                                                        >  4:232226/1‑135 (MQ=255)
                   ggcCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCCCGTCCGCCCGATCa                                                                                                                      >  7:202394/1‑136 (MQ=255)
                       tGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGCGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACGTCCGCCAGCTAACCAg                                                                                                                  >  5:57455/1‑136 (MQ=255)
                       tGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCAcccg                                                                                                                  >  6:214168/1‑134 (MQ=255)
                          cccACGATGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAgcgc                                                                                                               <  3:101028/136‑1 (MQ=255)
                                       cATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGCGCACTTCCGCCCGATCACCCGCCCCCGTGGCGCGCgg                                                                                                  >  2:5143/1‑136 (MQ=255)
                                             cAGCTTCCTGGGCCTGCCGGGCATCGCGGAGAGCGCCGAACCCGGGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGGGGGCGTGCGCACTTCAGCCAGACCACCAGCGCCCGTGGCGCGCGGCTGATc                                                                                            <  1:298654/136‑1 (MQ=255)
                                                 ttCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCCCTTCAGCCCGATCACCCGCGCCCGGGGCGCGCGGCGGCTCgccc                                                                                        >  3:37995/1‑136 (MQ=255)
                                                   ccTCGGCGTGCCGGCCCTCCGGGAGCGCCCCGAACCCGGGGCCCGGGCGCGCAACCTGGTGCTGGTGGAACTCTGCGGGGTGCTGGCGCACTTCAGCCAGACCCCCCGCGCCCGTGGCGCGCGGCTGATCgcccag                                                                                      <  4:229993/136‑1 (MQ=255)
                                                   ccTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGCGCCCATAAGCCCACTAACCACCGCCCCGGGCGCCCGGCAGCTAgcccac                                                                                      >  2:155245/1‑135 (MQ=255)
                                                     gcgccgtGCGGGGCATCCCGGAGAGCGCCTAAACCGTGGCCCGGGCGCGAAACCTGCTGCTGGGGCAACTGAGCGGCGTGCGTGCGCACTTCAGCCACACCCCCAGCGCCCGTGGCGCGCGGCTGATCgcccaggc                                                                                    <  2:56411/130‑1 (MQ=255)
                                                     gcGGCCTGCCGGGCATCCCGGAGAGCGCCGAACCGGGGGCCCGGGCGCGCAACCTGCTGCTGGGGGAACACCGCGGCGGGGGGGCGCACTACAGCCAGAACCCCCGCGCCCGGGGCGCGCGGCTGATCgcccaggc                                                                                    <  8:31371/135‑1 (MQ=255)
                                                       ggccgcccGGGCATCCCGGAGAGCGCCTAAACCGGGGCCCGGGCGCGCACCCGGCTGCGGGTGGAACCGCGCGGCGGGGGTGCGCAATCCCGCCAGAACACCCGCGCCCGGGGCGCGCGGCTGATCgcccaggc                                                                                    <  4:82272/128‑1 (MQ=255)
                                                       ggCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGGGGCCCTGGCGCGCAACGTGCGGCTGGTGGAACAGAGCGGGGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCgcccaggc                                                                                    <  5:195098/134‑1 (MQ=255)
                                                         ccggCCGGCCATCCCGGACCGCCCCGAACCGGGGGCCCTGGCGCGCAACCTGCTGCTGGGGTAACAGCGCGGGGGGCGGGCCCACTGCCGCCCGAGCCCCCGCGCCCGGGGCGCGCGGCTGATCgcccaggcccag                                                                                <  6:199519/133‑1 (MQ=255)
                                                               gggctgcccGGAGAGCGCCTAAACGGGGGCGGGGGCGCGCACCGTGGTGCTGGGGGAACCGCGGGGCGGGCGGGCGCAATCCCCCCCGACCACCCGGGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTg                                                                          <  2:107956/130‑1 (MQ=255)
                                                                       cGGTGAGCGCCTAAACCGGGGCCCTGGCGCGCACCCTGCGGCGGGTGGAACAGAGGGGCGTGCGGGCGCACTTCAGCCCGAGCCCCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGcc                                                                  <  5:69698/136‑1 (MQ=255)
                                                                        ggAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCg                                                                 <  7:255648/136‑1 (MQ=255)
                                                                             ggccgtaaaCGGGGGCCCTGGCGCGAAACTGCGTGGTGGTGCACCCGCGGGGCGGGCGGGGGCACTACAGCCAGACCACCCGCGCCGTGGGGGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAc                                                            <  2:266736/130‑1 (MQ=255)
                                                                             gcgcCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAc                                                            <  5:262354/136‑1 (MQ=255)
                                                                             gcgcCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCGGGTGGAACAGAGCGGCGGGCGGGCGCACTTCAGCCAGATCACCCGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAc                                                            <  7:52109/136‑1 (MQ=255)
                                                                             gcgcCGAAACCGGGGCCCTGGCGCGCAACCTGCTGCTGGGGGACCAGAGCGGCGTGCGTGCGCACTTCAGCCAGAGCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAc                                                            <  4:223893/136‑1 (MQ=255)
                                                                                       gggggCCGTGGCGCGCAACCTGCGGGTGGGGCAACCGCGGGGGGGGGGGGCGCACTGCAGCCAGATCCCCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                  <  3:201062/133‑1 (MQ=255)
                                                                                       ggTGGCCCTGGCGCGAAACGTGCTGCTGGTGTAACAGAGCGGCGTGCGTGCGAACTACCGCCATATCACCCGCGCCGGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                  <  5:251366/135‑1 (MQ=255)
                                                                                       cggggCCGTGGCGCGCAACCTGGTGCTGGTGGAACCGCGCGGCGGGCGGGCGCAATTCACCCAGACCCCCCGCGCGCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                  <  4:18293/133‑1 (MQ=255)
                                                                                       cggggCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                  <  6:274609/133‑1 (MQ=255)
                                                                                       cGTGGCCCTGGCGCGCACCCTGCTGCGGGGTGAACAGTGCGGCGGGGGTGCGAACTTCAGCCAGATCACCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                  <  2:90441/136‑1 (MQ=255)
                                                                                            cccTGGCGCGCAACCTGCTGCTGGGGTAACAGAGCGGCGGGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCt                                             <  6:62708/136‑1 (MQ=255)
                                                                                                   gcgcAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCCCTTCAGCCAGATCACCCGCGCCCGTGGCGCGCGGCTGCTCGCCCCGGCCCCGGCCGTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTAcccg                                      >  1:176015/1‑134 (MQ=255)
                                                                                                    cgcAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCCCATAAGCCCGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCCGGCCGTGGGCCTGCCGGTAACCGCCGCTGTAGCGCGGTTCCAGc                                     >  7:21306/1‑136 (MQ=255)
                                                                                                     gcAACCTGCGGCTGGGGGAACAGAGGGGCGTGGGTGCGCACTTCAGCCAGATCACCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCt                                    <  7:78405/136‑1 (MQ=255)
                                                                                                                            gagCGGCGTGCGTGCGCACTTCAGCCACACCACCCGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTg              <  1:37752/136‑1 (MQ=255)
                                                                                                                            gagCGGCGTGCGGGCGAACGGCACCCACATCACCCCCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCGGTACCAGCTGATCCTCACGGACGAGTCGCTg              <  6:84586/136‑1 (MQ=255)
                                                                                                                              gcggcgTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGCTCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACCGACGAGTCGCTgcg            >  7:78720/1‑136 (MQ=255)
                                                                                                                                   gtgcgtgcGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt       <  8:215020/136‑1 (MQ=255)
                                                                                                                                   gggggtgcgAACTTCAGCCATAGCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt       <  2:22748/132‑1 (MQ=255)
                                                                                                                                   gggcgtgcgAAGTTCGCCCAGATCCCCCGCGCCGGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt       <  4:294955/134‑1 (MQ=255)
                                                                                                                                   gggcgggcgCACTCCAGCCAGATCACCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt       <  2:70034/134‑1 (MQ=255)
                                                                                                                                   gggcgggccCACTTCCCCCCGACCACCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt       <  1:139922/134‑1 (MQ=255)
                                                                                                                                     gcgtgcgCACTTCAGCCAGATAACCAGCGCCCGTGGCGCGCGGCTGATAGCCCAGGCCCCGGCCTTGGGCCTGCCGGTAACCGCCGAGGTATCGCTGTACCAGCTGATCCTCACGGACTAGTCGCTGCGCGAGtta     >  2:81688/1‑135 (MQ=255)
                                                                                                                                     gcgtgcgCACTTCAGCCAGACCACCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTtc     <  6:203227/136‑1 (MQ=255)
                                                                                                                                     gcgtgcgCACATCCGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCGTGGGCCTTCCGGGAACCGCCGCTGTAGCGCTGTACCCGCTGATCCTCACGGACGCGTCGCTGCGCGAGTtc     >  7:150129/1‑136 (MQ=255)
                                                                                                                                        tgcgCACTTCAGCCAGATCACCAGCGCCCGTGGGGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTcc  <  8:25409/136‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                       
CCGCGACCTGGCCCAGGGCGGCCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCC  >  ACB143/5693018‑5693288

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