Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 5,693,167 A→C 19.7% Q233P (CAG→CCG)  pyrC' → Dihydroorotase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1435,693,1670AC19.7% 59.1 / 4.3 58Q233P (CAG→CCG) pyrC'Dihydroorotase‑like protein
Reads supporting (aligned to +/- strand):  ref base A (14/30);  new base C (9/4);  total (23/35)
Fisher's exact test for biased strand distribution p-value = 2.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

GGCCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGC‑TGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCC  >  ACB143/5693037‑5693288
                                                                                                                                   |                                                                                                                         
ggcCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGC‑TGGTGGAACAGAGCGGCGTGCGGGCGCCCGTCAGCCAGCTCa                                                                                                                      >  2:42936/1‑136 (MQ=255)
    tGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGC‑TGGTGGAACAGAGCGGCGTGGGTGCGGACTTAAGCCAGAtaccccg                                                                                                                  >  8:248880/1‑130 (MQ=255)
    tGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGC‑TGGTGGAACAGAGCGGCGTGCGGGCGCACATCCGCCACATACCCAg                                                                                                                  >  3:157338/1‑136 (MQ=255)
         cACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGC‑TGGTGGAACAGAGCGGCGTGCGTGGGCCCGTCAGCCCGATAACCAGCGccc                                                                                                             >  5:319702/1‑136 (MQ=255)
              aGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGC‑TGGTGGAACAGAGCGGCGTGCGTGCGCCCGTCAGCCAGATCACCAGCGCCCGTGgc                                                                                                        >  1:234859/1‑136 (MQ=255)
               ggTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGC‑TGGTGGAACAGAGCGGCGTGCGGGCGCCCTTCAGCCCGATCACCACCGCCCGTGgcg                                                                                                       >  7:71718/1‑136 (MQ=255)
                   ccATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCACCCTGCTGC‑TGGTGGAACAGAGCGGCGGGGGGGGGCCCTTCCGCCCCGTCACCCGCCCCCGGGGcgcgcg                                                                                                   >  4:85623/1‑136 (MQ=255)
                   ccATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGC‑TGGTGGAACAGAGCGGCGTGGGTGCGCAGATCAGCCAGATAAACCCCGCCCGGGGcgcgcg                                                                                                   >  6:95679/1‑136 (MQ=255)
                   ccATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGC‑TGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGcgcgcg                                                                                                   >  8:118491/1‑136 (MQ=255)
                    cATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGC‑TGGTGGAACAGAGCGGCGTGCGTGGGCACATAAACCCGAAAAACACCGCCCGTGgcgcgccg                                                                                                  >  4:178705/1‑134 (MQ=255)
                    cATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGC‑TGGTGGAACAGAGAGGCGTGCGGGCGCGCATGAGACAGATCAACCGCCCCCGTGgcgcccca                                                                                                  >  2:77925/1‑132 (MQ=255)
                          cAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCGGGCGCGCAACCTGCTGC‑TGGTGGACCAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGGGGCGCGCGGCTGATc                                                                                            <  8:128220/136‑1 (MQ=255)
                          cAGCTTCCTCGGCCGGCCGGGCATCCCGGAGAGCGCCTAAACCGTGGCCCTGGCGCGCAACCTGCTGC‑TGGTGGAACAGCGCGGCGTGCGGGCGCACTTCAGCCATATCACCAGCGCCCGTGGCGCGCGGCTGATc                                                                                            <  3:247480/136‑1 (MQ=255)
                           aGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGC‑TGGTGGATCAGAGCGGCGGGCGGGGGCCCGTCAGCCAGCTGCCCAGCCCCCGGGGCGCGCGGcggctag                                                                                           >  5:318228/1‑130 (MQ=255)
                                ccTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGC‑TGGTGGAACAGAGCGGCGTGCGGGGGGCCTTCAGCCCGATCACCCCCGCCCGTGGCGCGCGGCGGCTCgccccg                                                                                      >  7:58435/1‑134 (MQ=255)
                                 ctgggCGTGCCGGGCAGCCCGGAGCGCGCCTAAACCGGGGCCCTGGCGCGAAACCTGCTGC‑TGGGGTAACCGCGCGGCGGGCGGGCGCACTTCCGCCATACCCCCCGCGCCCGGGGCGCGCGGCTGATCgcccagg                                                                                     <  8:372233/133‑1 (MQ=255)
                                 cTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGC‑TGGTGGAACAGAGCGGCGTGCGTGCGAACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCgcccagg                                                                                     <  5:196037/136‑1 (MQ=255)
                                 cTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGC‑TGGTGGAACAGAGCGGCGGGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCgcccagg                                                                                     <  2:125293/136‑1 (MQ=255)
                                     gCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGC‑TGGTGGAACAGAGCGGCGTGCGTGGGCCCTTCAGCCCGCTCACCCGCGCCCGTGGCGCGCGGCGGATCgcccaggcc                                                                                   >  1:48431/1‑134 (MQ=255)
                                     gCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGC‑TGGTGGAACAGAGCGGCGTGCGTGCGCACATCAGCCAGATCACCCCCCCCCGTGGCGCGCGGCTGCTAgccccggcccc                                                                                 >  6:303693/1‑135 (MQ=255)
                                     gCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGC‑TGGTGGAACAGAGCGGCGTGCGGGCGCACATCAGGCACAGAACCCACCCCCGGGGGGCGCGGGGGCTCgcccacgcgca                                                                                 >  2:207625/1‑136 (MQ=255)
                                            gggCATCCCGGAGAGCGCCGAAACCGTGGCCCGGGCGCGCAACCTGCTGC‑GGGTGGACCAGAGCGGCGTGCGGGCGCACTTCAGCCAGCTCACCAGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTg                                                                          <  5:382688/136‑1 (MQ=255)
                                            gggCAGCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGC‑TGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGAGCACCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTg                                                                          <  1:339505/136‑1 (MQ=255)
                                                          gcgcCTAAACGGGGGCCCGGGCGCGCAACCTGCTGC‑TGGTGGAACCGCGGGGGGGGGGTGCGCAATTCACCCACATCCCCCCCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAc                                                            <  7:154438/136‑1 (MQ=255)
                                                          gcgcCTAAACCGTGGCCCGGGCGCGCAACCTGCTGG‑TGGGGCAACAGCGCGGGGGGCGGGCCCACTTCAGCCACACCCCCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAc                                                            <  6:211661/136‑1 (MQ=255)
                                                          gcgcCGAAACCGTGGCCCTGGCGCGCAACCTGCTGC‑TGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAc                                                            <  7:82271/136‑1 (MQ=255)
                                                           ggcCGAACCCGTGGCCCTGGCGCGCAACCTGCGGC‑TGGTGGAACAGAGCGGCGTGCGTGCGCATTTCAGCCAGAGCACCAGCGCCGGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAcc                                                           <  7:247526/135‑1 (MQ=255)
                                                               gaaccggTGGCCCGGGCGCGCACCCGGCTGG‑TGGTGGACCCGCGCGGGGGGCGTGCGCACTGCCCCCCGACCCCCCGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAccgccg                                                       <  1:206840/130‑1 (MQ=255)
                                                               gAAACGGTGGCCCGGGCGCGCAACCTGCTGC‑GGGTGGAACAGCGCGGCGTGCGGGCGCACTTCAGCCATATCACCAGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAccgccg                                                       <  7:206598/136‑1 (MQ=255)
                                                               gAAACCGTGGCCCTGGCGCGCAACCTGCTGC‑TGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAccgccg                                                       <  5:183900/136‑1 (MQ=255)
                                                               gAAACCGTGGCCCTGGCGCGCAACCTGCTGC‑TGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAccgccg                                                       <  7:118491/136‑1 (MQ=255)
                                                                 aaCCGTGGCCCTGGCGCGAAACCTGCTGC‑TGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGa                                                      <  2:334858/135‑1 (MQ=255)
                                                                   ccGTGGCCCTGGCGCGCAACCTGCTGC‑TGGTGGAACAGAGCGGCGTGCGGGGGGCCCTAAGCCCGCAAACCCACCCCCCGGGCGCGCCGCCGGTCGCCCCCGCCCCCGCCGTGGGCCGGCCGGGGAcccccgcggg                                                   >  2:167595/1‑132 (MQ=255)
                                                                    ggTGGCCCTGGCGCGCAACCTGCTGC‑TGGTTTACCATAGGGGCGGGGGTGCGCACTTCACCCATATACCCCCCCCCGGTGGGGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                  <  6:389901/135‑1 (MQ=255)
                                                                    cggggCCCTGGCGCGCAACCTGCTGC‑TGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCCGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                  <  8:317746/133‑1 (MQ=255)
                                                                    cgggccccGGGCGCGCAACCGGCGGC‑TGGGGCAACCTCCCGGGGGGCGGGCGCACTTCAGCCCGACCCCCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                  <  2:396533/131‑1 (MQ=255)
                                                                    cGTGGCCCTGGCGCGCAACGTGCTGC‑GGGTGGAACTGAGCGGCGGGCGTGCGCACTTCAGCCATATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                  <  4:304723/136‑1 (MQ=255)
                                                                    cGTGGCCCTGGCGCGCAACCTGCTGC‑TGGGGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                  <  7:65981/136‑1 (MQ=255)
                                                                    cGTGGCCCTGGCGCGAAACCTGCTGC‑TGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                  <  3:93616/136‑1 (MQ=255)
                                                                         cccTGGCGCGCAACCTGCTGC‑TGGTGGAACAGAGCGGCGTGCGTGCGCACGTCAGCCACATAACCACCCCCCGTTGCGCCCCGCGGCACGCCCCCGCCCACGCCGTGGGGCTGCCGGTAAACCCCGATTTAtcg                                               >  8:287091/1‑134 (MQ=255)
                                                                             ggCGCGAAACCGCCTGC‑TGGTGCAACCTCGCGGGGTGCTTGCGAACTCCTCCCCTATCACCCCCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTAc                                         <  8:324670/136‑1 (MQ=255)
                                                                              gcgcgcAACCTGCTGC‑TGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCCGGCCTTGGGCCTGCCGGGAACCGCCGATGTAGCGCTGTAcc                                        >  1:256237/1‑136 (MQ=255)
                                                                                 cgaacccTGCTGCGTGGT‑GAACATAGCGGCTTGCGGGCGCACTTCACCCAGACCACCCGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGc                                     <  4:303250/131‑1 (MQ=255)
                                                                                   cAACCTGCTGC‑TGGTGGAACAGAGCGGCGTGCGGGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTg                                   <  1:351070/136‑1 (MQ=255)
                                                                                     aCCTGGTGG‑GGGTGCAACAGCGGGGGGTGCGGGCGCATTCCAGCCAGAGCCCCCGGGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGAt                                 <  4:290884/136‑1 (MQ=255)
                                                                                                           gcggggGGCGGGCGCACTCCCGCCCCACCCCCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTgcg            <  6:2826/136‑1 (MQ=255)
                                                                                                            cggggTGCGTGCCCACTTCCCCCACATCCCCAGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTgcgc           <  5:14696/136‑1 (MQ=255)
                                                                                                             ggcgTGCGTGCGCACTTCAGCCAGATCACCCCCGCCCCTGGCGCGCGGCTGATCGCCCCGGCCCCCGCCGTGGGCCTGCCCGGAACCGCCGCTGTAGCGCTGTACCCCCTGATCCTCACGGGAGAGTCGCTGcgcg          >  8:60286/1‑136 (MQ=255)
                                                                                                             ggcgTGCGGGCGCACTTCCGCCCGCTCACCCGCGCCCGGGGCGCGCGGCTGATCGCCCCCGCCCCGGCCTTGGGCCTGCCGGTAACCGCCCCTGTAGCGCTGTACCAGCTGCTACTCACGGCCGCGTCgcggcggg          >  6:190152/1‑134 (MQ=255)
                                                                                                               cgtgcgtgCGCCCGTAAGCCAGATCACCCGCCCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCCGGAACCGCCGATGTAGCGCGGTACCAGCTGATACGCACGGGCGCGTCGCGGcgcgcg        >  4:261488/1‑134 (MQ=255)
                                                                                                                gtgcgtgcGCACTGCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt       <  8:83605/136‑1 (MQ=255)
                                                                                                                gtgcgggcGCACTTCAGCCAGATCACCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt       <  6:130176/136‑1 (MQ=255)
                                                                                                                gggggtgcgCAGTTCACCATAACCACCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt       <  2:357631/132‑1 (MQ=255)
                                                                                                                gggggtgcgAAGTTCCGCCAGATCACCCGCGCCGGGGGGGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt       <  6:319702/132‑1 (MQ=255)
                                                                                                                gggcgtgcgCACTTCCCCCATAACCCCCGGGCCGGTGGGGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt       <  3:92051/134‑1 (MQ=255)
                                                                                                                gggcgggcgCACTTCCGCCAGACCCCCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt       <  3:177453/134‑1 (MQ=255)
                                                                                                                  gcgtgcgCACTTAAGCCCGATCACCAGCCCCCGTGGCGCGCGGCTGATAGCCCAGGCCCAGGCCTTGGGCCTGCCGGGAACCGCCGATGTAGCGCGGTACCACCAGCTCCTCCCGGCCGAGTCGCTGCGCGAGGtc     >  6:297875/1‑136 (MQ=255)
                                                                                                                     tgcgCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTcc  <  7:110464/136‑1 (MQ=255)
                                                                                                                                   |                                                                                                                         
GGCCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGC‑TGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCC  >  ACB143/5693037‑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: ‑