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

TGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGC‑TGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAG  >  ACB143/5693041‑5693306
                                                                                                                                       |                                                                                                                                   
tGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGCGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGC‑TGCTGGTGGAACAGAGCGGCGTGCGTGCGCACGTCCGCCAGCTAACCAg                                                                                                                                    >  5:57455/1‑136 (MQ=255)
   cccACGATGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGC‑TGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAgcgc                                                                                                                                 <  3:101028/136‑1 (MQ=255)
                cATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGC‑TGCTGGTGGAACAGAGCGGCGTGCGGGCGCACTTCCGCCCGATCACCCGCCCCCGTGGCGCGCgg                                                                                                                    >  2:5143/1‑136 (MQ=255)
                      cAGCTTCCTGGGCCTGCCGGGCATCGCGGAGAGCGCCGAACCCGGGGCCCTGGCGCGCAACCTGC‑TGCTGGTGGAACAGAGCGGGGGGCGTGCGCACTTCAGCCAGACCACCAGCGCCCGTGGCGCGCGGCTGATc                                                                                                              <  1:298654/136‑1 (MQ=255)
                          ttCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGC‑TGCTGGTGGAACAGAGCGGCGTGCGTGCGCCCTTCAGCCCGATCACCCGCGCCCGGGGCGCGCGGCGGCTCgccc                                                                                                          >  3:37995/1‑136 (MQ=255)
                            ccTCGGCGTGCCGGCCCTCCGGGAGCGCCCCGAACCCGGGGCCCGGGCGCGCAACCTGG‑TGCTGGTGGAACTCTGCGGGGTGCTGGCGCACTTCAGCCAGACCCCCCGCGCCCGTGGCGCGCGGCTGATCgcccag                                                                                                        <  4:229993/136‑1 (MQ=255)
                            ccTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGC‑TGCTGGTGGAACAGAGCGGCGTGCGGGCGCCCATAAGCCCACTAACCACCGCCCCGGGCGCCCGGCAGCTAgcccac                                                                                                        >  2:155245/1‑135 (MQ=255)
                              gcgccgtGCGGGGCATCCCGGAGAGCGCCTAAACCGTGGCCCGGGCGCGAAACCTGC‑TGCTGGGGCAACTGAGCGGCGTGCGTGCGCACTTCAGCCACACCCCCAGCGCCCGTGGCGCGCGGCTGATCgcccaggc                                                                                                      <  2:56411/130‑1 (MQ=255)
                              gcGGCCTGCCGGGCATCCCGGAGAGCGCCGAACCGGGGGCCCGGGCGCGCAACCTGC‑TGCTGGGGGAACACCGCGGCGGGGGGGCGCACTACAGCCAGAACCCCCGCGCCCGGGGCGCGCGGCTGATCgcccaggc                                                                                                      <  8:31371/135‑1 (MQ=255)
                                ggccgcccGGGCATCCCGGAGAGCGCCTAAACCGGGGCCCGGGCGCGCACCCGGC‑TGCGGGTGGAACCGCGCGGCGGGGGTGCGCAATCCCGCCAGAACACCCGCGCCCGGGGCGCGCGGCTGATCgcccaggc                                                                                                      <  4:82272/128‑1 (MQ=255)
                                ggCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGGGGCCCTGGCGCGCAACGTGC‑GGCTGGTGGAACAGAGCGGGGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCgcccaggc                                                                                                      <  5:195098/134‑1 (MQ=255)
                                  ccggCCGGCCATCCCGGACCGCCCCGAACCGGGGGCCCTGGCGCGCAACCTGC‑TGCTGGGGTAACAGCGCGGGGGGCGGGCCCACTGCCGCCCGAGCCCCCGCGCCCGGGGCGCGCGGCTGATCgcccaggcccag                                                                                                  <  6:199519/133‑1 (MQ=255)
                                        gggctgcccGGAGAGCGCCTAAACGGGGGCGGGGGCGCGCACCGTGG‑TGCTGGGGGAACCGCGGGGCGGGCGGGCGCAATCCCCCCCGACCACCCGGGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTg                                                                                            <  2:107956/130‑1 (MQ=255)
                                                cGGTGAGCGCCTAAACCGGGGCCCTGGCGCGCACCCTGC‑GGCGGGTGGAACAGAGGGGCGTGCGGGCGCACTTCAGCCCGAGCCCCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGcc                                                                                    <  5:69698/136‑1 (MQ=255)
                                                 ggAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGC‑TGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCg                                                                                   <  7:255648/136‑1 (MQ=255)
                                                      ggccgtaaaCGGGGGCCCTGGCGCG‑AAACTGCGTGGTGGTGCACCCGCGGGGCGGGCGGGGGCACTACAGCCAGACCACCCGCGCCGTGGGGGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAc                                                                              <  2:266736/130‑1 (MQ=255)
                                                      gcgcCGAAACCGTGGCCCTGGCGCGCAACCTGC‑TGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAc                                                                              <  5:262354/136‑1 (MQ=255)
                                                      gcgcCGAAACCGTGGCCCTGGCGCGCAACCTGC‑TGCGGGTGGAACAGAGCGGCGGGCGGGCGCACTTCAGCCAGATCACCCGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAc                                                                              <  7:52109/136‑1 (MQ=255)
                                                      gcgcCGAAACCGGGGCCCTGGCGCGCAACCTGC‑TGCTGGGGGACCAGAGCGGCGTGCGTGCGCACTTCAGCCAGAGCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAc                                                                              <  4:223893/136‑1 (MQ=255)
                                                                gggggCCGTGGCGCGCAACCTGC‑GGGTGGGGCAACCGCGGGGGGGGGGGGCGCACTGCAGCCAGATCCCCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                    <  3:201062/133‑1 (MQ=255)
                                                                ggTGGCCCTGGCGCGAAACGTGC‑TGCTGGTGTAACAGAGCGGCGTGCGTGCGAACTACCGCCATATCACCCGCGCCGGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                    <  5:251366/135‑1 (MQ=255)
                                                                cggggCCGTGGCGCGCAACCTGG‑TGCTGGTGGAACCGCGCGGCGGGCGGGCGCAATTCACCCAGACCCCCCGCGCGCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                    <  4:18293/133‑1 (MQ=255)
                                                                cggggCCCTGGCGCGCAACCTGC‑TGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                    <  6:274609/133‑1 (MQ=255)
                                                                cGTGGCCCTGGCGCGCACCCTGC‑TGCGGGGTGAACAGTGCGGCGGGGGTGCGAACTTCAGCCAGATCACCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                    <  2:90441/136‑1 (MQ=255)
                                                                     cccTGGCGCGCAACCTGC‑TGCTGGGGTAACAGAGCGGCGGGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCt                                                               <  6:62708/136‑1 (MQ=255)
                                                                            gcgcAACCTGC‑TGCTGGTGGAACAGAGCGGCGTGCGTGCGCCCTTCAGCCAGATCACCCGCGCCCGTGGCGCGCGGCTGCTCGCCCCGGCCCCGGCCGTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTAcccg                                                        >  1:176015/1‑134 (MQ=255)
                                                                             cgcAACCTGC‑TGCTGGTGGAACAGAGCGGCGTGCGTGCGCCCATAAGCCCGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCCGGCCGTGGGCCTGCCGGTAACCGCCGCTGTAGCGCGGTTCCAGc                                                       >  7:21306/1‑136 (MQ=255)
                                                                              gcAACCTGC‑GGCTGGGGGAACAGAGGGGCGTGGGTGCGCACTTCAGCCAGATCACCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCt                                                      <  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)
                                                                                                                                       |                                                                                                                                   
TGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGC‑TGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAG  >  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: ‑