Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 30,273 T→C 6.5% G593G (GGT→GGC carB → carbamoyl‑phosphate synthase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported30,2730TC 6.5% 169.7 / 3.9 61G593G (GGT→GGCcarBcarbamoyl‑phosphate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base T (27/30);  new base C (2/2);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

AGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACCGC  >  Exported/30215‑30339
                                                          |                                                                  
aGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCAt                                                         >  1:2246741/1‑70 (MQ=255)
aGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCAt                                                         >  1:4013542/1‑70 (MQ=255)
aGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCAt                                                         <  1:2972580/70‑1 (MQ=255)
aGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCAt                                                         >  1:16251/1‑70 (MQ=255)
 gggTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCAtt                                                        <  1:1827662/70‑1 (MQ=255)
 gggTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTa                                                       <  1:1037913/71‑1 (MQ=255)
  ggcaTCGAGTTCGACTACTGCTGCGTGCACGCGGCGTTGGCGATGCGCGAAGACGGCTACGAGACcatcat                                                      <  1:3723251/68‑1 (MQ=255)
        gAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGa                                                               <  1:2694759/56‑1 (MQ=255)
        gAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGa                                                               <  1:1131595/56‑1 (MQ=255)
        gAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTa                                                       >  1:1139172/1‑64 (MQ=255)
        gAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTaa                                                >  1:2304589/1‑71 (MQ=255)
         aaTTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAAc                                               >  1:1682105/1‑71 (MQ=255)
         aaTTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAAc                                               >  1:2570048/1‑71 (MQ=255)
         aaTTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAAc                                               >  1:3545945/1‑71 (MQ=255)
         aaTTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAAc                                               >  1:738083/1‑71 (MQ=255)
          gttCGACTACTGCTGCGTGCACGCCGCCCTCGCCATGCGCGAAGACGGCTACGAAACCATCATGGTCAACt                                              <  1:4201379/70‑1 (MQ=255)
           ttCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTAc                                                                 <  1:512885/51‑1 (MQ=255)
           ttCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTAc                                                                 <  1:2591/51‑1 (MQ=255)
           ttCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAAc                                                            <  1:501100/56‑1 (MQ=255)
                   ctGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAAcc                                        <  1:912572/68‑1 (MQ=255)
                   ctGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCg                                      <  1:4101332/70‑1 (MQ=255)
                    tGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCgg                                     >  1:1705805/1‑70 (MQ=255)
                    tGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCAgg                                     >  1:3129122/1‑70 (MQ=255)
                      ttGCGTACACGCCTCGCTGGCGCTGCGCGATGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGa                                    <  1:2294315/69‑1 (MQ=255)
                      ttGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAACCCATTATGGTTAACTGTAACCCGGaaa                                  >  1:422740/1‑71 (MQ=255)
                      ttGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGaa                                   <  1:1694610/70‑1 (MQ=255)
                      ttGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGa                                    <  1:713255/69‑1 (MQ=255)
                      ctGCGTGCACGCCGCCCTCGCCATGCGCGAAGACGGCTACGaaa                                                             >  1:1841063/2‑44 (MQ=12)
                      ctGCGTGCACGCCGCCCTCGCCATGCGCGAAGACGGCTACGAAACCATCATGGTCAACTGCAACCCCGaga                                  >  1:3010605/2‑69 (MQ=255)
                       tGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGaa                                   <  1:124875/69‑1 (MQ=255)
                        gCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAAGTGTAACCCGGAAAc                                 <  1:4293684/70‑1 (MQ=255)
                        gCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAAc                                 <  1:487222/70‑1 (MQ=255)
                         cGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAAcc                                <  1:2772537/70‑1 (MQ=255)
                          gTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGt                              <  1:188231/71‑1 (MQ=255)
                           tACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGtc                             <  1:1708802/71‑1 (MQ=255)
                             cacGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGtctc                           <  1:609562/71‑1 (MQ=255)
                             cacGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGtctc                           <  1:3288903/71‑1 (MQ=255)
                             cacGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGtctc                           <  1:2164209/71‑1 (MQ=255)
                                gCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGtt                                                  <  1:1760819/45‑1 (MQ=255)
                                   tCGCTGGCGCTGCGCGAAGACGGTTACGAAACCAt                                                         <  1:329136/35‑1 (MQ=255)
                                   tCGCTGGCGCTGCGCGAAGACGGTTACGAAACCAt                                                         <  1:1930321/35‑1 (MQ=255)
                                    cGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTaa                                                >  1:2369383/1‑43 (MQ=255)
                                    cGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTaa                                                >  1:4412662/1‑43 (MQ=255)
                                    cGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGag                    >  1:1328076/1‑70 (MQ=255)
                                     gCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTaa                                                >  1:4065620/1‑42 (MQ=255)
                                     gCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGAc                    >  1:2340269/1‑70 (MQ=255)
                                        ggCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGaaa                                  >  1:3301842/1‑53 (MQ=255)
                                        ggCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTAc                 >  1:4109237/1‑70 (MQ=255)
                                         gcgcTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGa               >  1:3133679/1‑71 (MQ=255)
                                         gcgcTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGa               >  1:2703623/1‑71 (MQ=255)
                                               cgcgAAGACGGTTACGAAACCATTATGGTTAACTGTaa                                          <  1:652843/38‑1 (MQ=255)
                                                 cgAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTcc        >  1:70327/1‑70 (MQ=255)
                                                 cgAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTcc        <  1:424676/70‑1 (MQ=255)
                                                 cgAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTcc        <  1:2103137/70‑1 (MQ=255)
                                                    aGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGAcc    >  1:2498252/1‑71 (MQ=255)
                                                    aGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGAc     >  1:1139407/1‑70 (MQ=255)
                                                     gACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACCg   >  1:1071186/1‑71 (MQ=255)
                                                     gACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACCg   >  1:3803574/1‑71 (MQ=255)
                                                     gACGGTTACGAAACCATTATGGTTAACTGTAACACGGAAACCGTCTCCACCGACTACGACACTTCCGACCg   >  1:3592400/1‑71 (MQ=255)
                                                      aCGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACCGc  <  1:4274364/71‑1 (MQ=255)
                                                       cGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTAc                 <  1:2066685/55‑1 (MQ=255)
                                                          |                                                                  
AGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACCGC  >  Exported/30215‑30339

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