Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,688,637 A→C 71.3% S236R (AGC→CGC)  pyrC' → Dihydroorotase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,688,6370AC71.3% 4.6 / 12.3 42S236R (AGC→CGC) pyrC'Dihydroorotase‑like protein
Reads supporting (aligned to +/- strand):  ref base A (7/5);  new base C (7/21);  total (14/28)
Fisher's exact test for biased strand distribution p-value = 7.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.55e-01

TGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAGCC  >  ACB122/5688503‑5688770
                                                                                                                                      |                                                                                                                                     
tGGCCCACGAAGGTGCCATGGCCCGCTTCCTCGGCCTGCCGGGCATCCCGGAGCGCGCCGAAACCGTGGCCCTGGCGCGCACCCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGCGCACGTCAGCCCGCTAACCa                                                                                                                                       >  6:145817/1‑135 (MQ=255)
      aCGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGGGCCCGTCAGCCCCCTCACCACCGccc                                                                                                                                 >  3:314060/1‑135 (MQ=255)
          aGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGCGCCCCTCCGCCCGATAACCAGCGCCCGTgg                                                                                                                             >  2:126442/1‑135 (MQ=255)
                  tGGCCAGCTTCCTCGGCCGGCCGGGCATCCCGGAGCGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGGGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGc                                                                                                                     <  3:382744/135‑1 (MQ=255)
                       aGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGCGCCCTTCAGCCAGATCACCCGCGCCCGGGGCGCGCGGCGGATCg                                                                                                               >  7:212660/1‑136 (MQ=255)
                            ccTCGGCCTGCCGGGCATCCAGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGT‑GAACAGAGCGGCGTGGGGGGGCACGTGCGCCATAACCCCCCAGCCCCTGGCGCCCCGCGGATCgccccg                                                                                                          >  6:70234/1‑133 (MQ=255)
                              tCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCGGCGGGTGGAACAGAGCGGCGGGCGTGCGCACTTCCGCCAGATCACCCGCGCCCGGGGCGCGCGGCTGATCgcccagg                                                                                                         <  3:186788/135‑1 (MQ=255)
                                ggCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGCGCCCGTCCCCCACATCACCCCCGCCCGGGGCGCCCGGGGGGTCgccccggcc                                                                                                       >  2:246654/1‑135 (MQ=255)
                                      ccGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGGGGGGGCCCTGGAGCCACACCCCCACCGCCCCGGGCGCGCGGGGGCTCgcccccgcccccgcc                                                                                                 >  3:253619/1‑130 (MQ=255)
                                                 gcagagCGCCTAAACGGGGGCCCGGGCGCGAAACCTGGTGCTGGTGGAACCGAGCGGGGTGCGTGCGCACTTCCGCCAGATCCCCCGCGCCGGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGcc                                                                                      <  3:123163/133‑1 (MQ=255)
                                                      gcgcCGAACCGGGGGCCGGGGGGCGCACCGTGGGGGTGGTGGAACCGCGGGGGGGGGGGGCGCACTTCCGCCACAGCCCCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAc                                                                                <  8:364602/136‑1 (MQ=255)
                                                      gcgcCGAACCGGGGGCCCTGGCGCGCAACCGGCTGCTGGTGTAACCGAGCGGCGTGGGTGCGCACTTCAGCCCTAACCCCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAc                                                                                <  7:315396/136‑1 (MQ=255)
                                                        gcCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCCCTTCAGCCAGATCACCCGCCCCCCTGGCGCGCGGCTGCTCGCCCCCGCCCCGGCCGTGGGGCTGCCgggtacc                                                                               >  1:32848/1‑130 (MQ=255)
                                                          cGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGGGGGGGGCCTGGAGCCACCTAACCCGCGCCCGGGGCGCGCGGCGGCTCGCCCCGGCCCCGGCCGTGGGGCGGCCGGGGAccgg                                                                             >  5:187535/1‑134 (MQ=255)
                                                          cGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGGGGGGGGCCGTCAGCCAGATCAACCCCCCCCCGGGCGCGCCGCGGGACGCCCCCGCCCCGGGCGTGGGCCTGCCCGGAAcccc                                                                             >  6:359857/1‑133 (MQ=255)
                                                          cGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGGGCACTTCAGCCCGAACACCCCCCCCCCTGGGGCGCGGCGGCTAGCCCCCGCCCACGCCGTGGGCCAGCCCGGAAccgc                                                                             >  6:112585/1‑135 (MQ=255)
                                                            aaaCCGGGGCCCTGGCGCGCAACCTGCTGCTGGGGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGAGCCCCCGCGCCGGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAccgccg                                                                           <  3:53954/135‑1 (MQ=255)
                                                             aaCGGTGGCCGTGGCGGGCAACCGGGTGCGGGTGCACCAGAGCGGCGGGCGGGCGCAGTTCAGCCAGACCACCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGa                                                                          <  4:215227/135‑1 (MQ=255)
                                                                gggggCCCGGGCGCGAAACCTGCGGGGGGGGCAACAGAGCGGCGGGCGGGCGCACTTCAGCCAGATCCCCAGCGCCGGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                      <  7:257385/133‑1 (MQ=255)
                                                                cggggCCCTGGCGCGCAACCTGCTGCGGGGGGAACAGAGCGGCGGGCGTGCGCACTTCCGCCACAGCCCCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                      <  8:370936/133‑1 (MQ=255)
                                                                 gTGGCCCTGGCGCGCAACCTGCTGCGGGGGGAACAGCGCGGCGTGCGTGCGCACTTCAGCCAGATCCCCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                      <  2:32848/135‑1 (MQ=255)
                                                                 gTGGCCCTGGCGCGCAACCGGGTGGGGGGGGAACTGAGGGGGGTGCGTGCGCACTGCTGCCAGAGCCCCGGCGCCCGTGGGGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                      <  4:314060/135‑1 (MQ=255)
                                                                              gcAACCTGGGGCTGGGGAACCAGCGCGGCGGGCGGGCGCACTTCCGCCAGATCACCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCt                                                        <  7:322633/136‑1 (MQ=255)
                                                                              gaaaCCTGCGGCTGGGGGAACCGGGCGGGGGGCGGGGGCACTTCTGCCAGAGCACCCGGGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGc                                                         <  1:246654/133‑1 (MQ=255)
                                                                                  ccTGCTGCTGGTGGAACAGAGCGGCGTGCGAGCGCCCTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGCTCGCCCCGGCCCCGGCCGTGGGCCTGCCCGTAACCGCCGATGTAGCGCTGTACCAGCTGAt                                                     >  1:57948/1‑135 (MQ=255)
                                                                                                  cAGAGCGGCGTGCGTGCGCACATCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCCGGCCGTGGGCCTGCCCGTAACCCCCCATGTAGCGCTGTACCAGCTGCTACTCACGGCCGAGTCg                                     >  8:319947/1‑135 (MQ=255)
                                                                                                     cgcggcgGGCGTGCGCACTTCCGCCAGAGCACCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTg                                  <  1:3372/134‑1 (MQ=255)
                                                                                                      gcggcgTGCGGGCGCACTCCCCCCAGACCACCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTgcg                                <  7:201757/136‑1 (MQ=255)
                                                                                                      gcggcgGGCGGGCGCACCCCCGCCACACCCCCGCCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTgcg                                <  7:346479/136‑1 (MQ=255)
                                                                                                       cggcgTGCGTGCGCACTTCAGCCCCAGCCCCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTgcg                                <  2:57948/135‑1 (MQ=255)
                                                                                                           gggcgggcgCACTTCCGCCATATCACCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                           <  7:319947/134‑1 (MQ=255)
                                                                                                           gggcgggcgCACTCCAGCCAGATCCCCAGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                           <  8:79315/134‑1 (MQ=255)
                                                                                                            tgggggcgcACTCCAGCCAGACCCCCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                           <  5:112585/130‑1 (MQ=255)
                                                                                                            tgcgtgcgAACTGCAGCCAGAGCACCCGGGCCGGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                           <  6:125761/135‑1 (MQ=255)
                                                                                                            ggcgggcgCACTCCCCCCAGAACACCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                           <  6:305021/134‑1 (MQ=255)
                                                                                                             gcgggcgCACTTCCGCCAGAGCACCCGCGCCCGTGCGGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGtt                          <  4:127569/135‑1 (MQ=255)
                                                                                                             gcgggcgCACTCCCGCCAGACCACCCGCGCCGGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTtc                         <  8:366354/136‑1 (MQ=255)
                                                                                                              cgtgccAACTGCAGCCAGAGCCCCCGCGCCCGGGCCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTATCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTtc                         <  6:190352/135‑1 (MQ=255)
                                                                                                              cgggcgcACTTCAGCCACTTCACCCGGGCCGGTGGGGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATGCTCACGGACGAGTCGCTGCGCGAGTtc                         <  6:304293/132‑1 (MQ=255)
                                                                                                                tgcgCAGTTCAGCCAGATCACCAGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTcc                      <  8:84320/136‑1 (MQ=255)
                                                                                                                                    ccAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTtc                         >  4:312524/1‑113 (MQ=255)
                                                                                                                                    ccAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTtc                         <  3:312524/113‑1 (MQ=255)
                                                                                                                                    ccAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAgcc  >  7:126775/1‑136 (MQ=255)
                                                                                                                                      |                                                                                                                                     
TGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAGCC  >  ACB122/5688503‑5688770

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