Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,688,615 T→G 44.3% R228R (CGT→CGG pyrC' → Dihydroorotase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,688,6150TG44.3% 2.2 / 11.0 35R228R (CGT→CGGpyrC'Dihydroorotase‑like protein
Reads supporting (aligned to +/- strand):  ref base T (5/13);  new base G (9/7);  total (15/20)
Fisher's exact test for biased strand distribution p-value = 1.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CCCAGGGCGGCCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCC  >  ACB122/5688491‑5688750
                                                                                                                            |                                                                                                                                       
cccAGGGCGGCCTGGCCCACGACGGTGCGATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCTACCTGCTGCTGGGGGAACAGAGAGGCGTGCGTGCGCACTTc                                                                                                                                <  1:192999/134‑1 (MQ=255)
cccAGGGCGGCCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGGGGgcgcgcgtaa                                                                                                                               >  5:125761/1‑129 (MQ=255)
 ccAGGGCGGCCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAg                                                                                                                              <  3:340463/135‑1 (MQ=255)
        ggcCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGACCAGAGCGGCGTGCGTGGGCCCTTGAGACCGAtg                                                                                                                       >  4:340463/1‑134 (MQ=255)
        ggcCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATc                                                                                                                       <  4:124879/135‑1 (MQ=255)
            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)
                                        ccTCGGCCTGCCGGGCATCCAGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGAACAGAGCGGCGTGGGGGGGCACGTGCGCCATAACCCCCCAGCCCCTGGCGCCCCGCGGATCgccccg                                                                                      >  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)
                                                                                                                        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)
                                                                                                                            tgcgCAGTTCAGCCAGATCACCAGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTcc  <  8:84320/136‑1 (MQ=255)
                                                                                                                            |                                                                                                                                       
CCCAGGGCGGCCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCC  >  ACB122/5688491‑5688750

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: