Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,482,473 A→G 35.5% L174P (CTC→CCC)  argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,482,4730AG35.5% 57.0 / 62.5 62L174P (CTC→CCC) argIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (15/25);  new base G (6/16);  total (21/41)
Fisher's exact test for biased strand distribution p-value = 5.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

CTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTT  >  W3110S.gb/4482408‑4482540
                                                                 |                                                                   
cTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATc                                                                <  1:3586380/71‑1 (MQ=255)
cTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATc                                                                <  1:3027827/71‑1 (MQ=255)
  tCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGGGCATCg                                                               >  1:1441111/1‑70 (MQ=255)
  tCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCg                                                               >  1:1777656/1‑70 (MQ=255)
   ccGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCTAGCATCg                                                               <  1:2008648/69‑1 (MQ=255)
   ccGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGGGCATCg                                                               <  1:359229/69‑1 (MQ=255)
   ccGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGa                                                              <  1:3221224/70‑1 (MQ=255)
   ccGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGa                                                              <  1:559358/70‑1 (MQ=255)
     ggCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAAtt                                                           <  1:3909205/71‑1 (MQ=255)
        catcACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGcc                                                        <  1:3315014/68‑1 (MQ=255)
        cAGCACGCTTGTGGCGCGACCAGCCGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAAt                                                            >  1:2170953/1‑67 (MQ=255)
        cAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGaa                                                             >  1:3113928/1‑66 (MQ=255)
        cAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGcc                                                        >  1:2397515/1‑71 (MQ=255)
        cAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGcc                                                        >  1:3133086/1‑71 (MQ=255)
         aGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGGGCATCGa                                                              <  1:1020091/64‑1 (MQ=255)
          gCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCg                                                               <  1:2357074/62‑1 (MQ=255)
          gCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCg                                                               <  1:533882/62‑1 (MQ=255)
            aCGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGGGCATCGAATTGCCCAtg                                                    <  1:3722393/71‑1 (MQ=255)
             cGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtgt                                                   <  1:805743/71‑1 (MQ=255)
               cTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtgttg                                                 >  1:3558747/1‑71 (MQ=255)
               cTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtgttg                                                 >  1:1358597/1‑71 (MQ=255)
               cTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtgttg                                                 >  1:1121607/1‑71 (MQ=255)
                  gtgGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGGGCATCGAATTGCCCAtgttgt                                                <  1:1419693/69‑1 (MQ=255)
                  gtgGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGGGCATCGAATTGCCCAtgttgt                                                <  1:1940976/69‑1 (MQ=255)
                  gtgGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGGGCATCGAATTGCCCAtgttgt                                                <  1:933068/69‑1 (MQ=255)
                  gtgGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtgttgt                                                <  1:2360019/69‑1 (MQ=255)
                  gtgGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtgttgt                                                <  1:2007480/69‑1 (MQ=255)
                  gtgGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtgttgt                                                <  1:1846635/69‑1 (MQ=255)
                   tgGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGGGCATCGAATTGCCCAtgttgt                                                <  1:4005819/68‑1 (MQ=255)
                         gACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGt                                       <  1:1275004/71‑1 (MQ=255)
                         gACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGt                                       <  1:3206813/71‑1 (MQ=255)
                            cAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCAc                                    >  1:3738018/1‑71 (MQ=255)
                            cAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCAc                                    >  1:1703809/1‑71 (MQ=255)
                            cAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCAc                                    >  1:3055417/1‑71 (MQ=255)
                             aGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCAcc                                   >  1:2443880/1‑71 (MQ=255)
                               aCGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGGGCATCGAATTGCCCATgttgtt                                               <  1:1875266/57‑1 (MQ=255)
                               aCGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGGGCATCGAATTGCCCATgttgtt                                               <  1:185672/57‑1 (MQ=255)
                               aCGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGGGCATCGAATTGCCCATgttgtt                                               <  1:1039576/57‑1 (MQ=255)
                               aCGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCg                                                               <  1:3959890/41‑1 (MQ=255)
                               aCGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTg                                 <  1:2800042/71‑1 (MQ=255)
                               aCGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTg                                 <  1:1045391/71‑1 (MQ=255)
                               aCGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTg                                 <  1:1868877/71‑1 (MQ=255)
                                cGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCAc                                    >  1:1021120/1‑67 (MQ=255)
                                 gCAAATCCAGACCGGTAAGCGCCGCAGCTTCGGGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                               <  1:3860745/71‑1 (MQ=255)
                                 gCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGc                                                         >  1:2854748/1‑45 (MQ=255)
                                 gCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                               <  1:606048/71‑1 (MQ=255)
                                  cAAATCCAGACCGGTAAGCGCCGCAGCTTCGGGCATCGa                                                              <  1:917782/39‑1 (MQ=255)
                                  cAAATCCAGACCGGTAAGCGCCGCAGCTTCGGGCATCGa                                                              <  1:222683/39‑1 (MQ=255)
                                  cAAATCCAGACCGGTAAGCGCCGCAGCTTCGGGCATCGa                                                              <  1:2098247/39‑1 (MQ=255)
                                        cAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTa                                              >  1:3985616/1‑49 (MQ=255)
                                             cGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTc                   <  1:1602788/71‑1 (MQ=255)
                                                  aGCGCCGCAGCTTCGGGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTc              <  1:1199393/71‑1 (MQ=255)
                                                  aGCGCCGCAGCTTCGGGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGGCATTTc              <  1:2798611/71‑1 (MQ=255)
                                                   gcgcCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGAc                          <  1:2046356/58‑1 (MQ=255)
                                                   gcgcCGCAGCTTCGAGCATCGAATTGCCCATGTTCTTACGCGCGTCACCTGCATAGAc                          <  1:846129/58‑1 (MQ=255)
                                                      ccgcAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTg          <  1:149647/71‑1 (MQ=255)
                                                        gcAGCTTCGGGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATa                             >  1:3192031/1‑50 (MQ=255)
                                                        gcAGCTTCGGGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATa                             >  1:3436776/1‑50 (MQ=255)
                                                        gcAGCTTCGGGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATa                             >  1:3160001/1‑50 (MQ=255)
                                                              tCGGGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCAttt               >  1:1604469/1‑58 (MQ=255)
                                                              tCGGGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCttt  >  1:720538/1‑71 (MQ=255)
                                                                gAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATa                             <  1:42649/42‑1 (MQ=255)
                                                                 |                                                                   
CTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTT  >  W3110S.gb/4482408‑4482540

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