Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 5,693,175 A→C 54.2% S236R (AGC→CGC)  pyrC' → Dihydroorotase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1435,693,1750AC54.2% ‑4.2 / 18.8 44S236R (AGC→CGC) pyrC'Dihydroorotase‑like protein
Reads supporting (aligned to +/- strand):  ref base A (7/13);  new base C (3/21);  total (10/34)
Fisher's exact test for biased strand distribution p-value = 1.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.12e-01

TGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAG  >  ACB143/5693041‑5693306
                                                                                                                                       |                                                                                                                                   
tGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGCGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACGTCCGCCAGCTAACCAg                                                                                                                                    >  5:57455/1‑136 (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)
                                                                                                                gggcgcACTTCAGCAAGATCACCAGCGCCCTGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTtctc                     <  1:259860/134‑1 (MQ=255)
                                                                                                                 tgcgCACTTCAGCCAGATCACCAGCGCCCGTGGGGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTcc                    <  8:25409/136‑1 (MQ=255)
                                                                                                                 ggcgcCCGTCAGCCAGAACACCCGCGCCCGTGGCCCGCGGCTGATCGCCCAGGCCCAGTCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTcc                    <  6:129159/135‑1 (MQ=255)
                                                                                                                                gATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTACCCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTg     >  3:184818/1‑136 (MQ=255)
                                                                                                                                   cACCCGGGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAg  <  3:70107/136‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                   
TGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAG  >  ACB143/5693041‑5693306

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