Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,455,055 T→C 9.8% K206K (AAA→AAG)  pyrB ← aspartate carbamoyltransferase catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,455,0550TC 9.8% 106.6 / 4.9 41K206K (AAA→AAG) pyrBaspartate carbamoyltransferase catalytic subunit
Reads supporting (aligned to +/- strand):  ref base T (18/18);  new base C (2/2);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01

GATCGCTGGCGCGAAGAACAAACTGCGCTTTCACGTTGGCGTACTCGGACGGGTCCAGACGCTCTTTTTGCACGCGGGTCATGTACAGGATGTCAACTTCCGCCATCACTTCTTCAATAGAGCTGTGCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACG  >  REL606/4454910‑4455202
                                                                                                                                                 |                                                                                                                                                   
gaTCGCTGGCGCGAAGAACAAACTGCGCTTTCACGTTGGCGTACTCGGACGGGTCCAGACGCTCTTTTTGCACGCGGGTCATGTACAGGATGTCAACTTCCGCCATCACTTCTTCAATAGAGCTGTGCAGACTCCAGGCGATCCCttttt                                                                                                                                                 >  2:196588/1‑150 (MQ=255)
                  ccacacgcgcTTTCACGTTTGCGCACTCGGACGGGTCCAGACGCTCTTTTTGCACGCGGGTCATGTACAGGATGTCAACTTCCGCCATCACTTCTTCAATAGAGCTGTGCAGTCTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGaa                                                                                                                               <  1:450682/144‑1 (MQ=255)
                   aaaCTGCGCTTTCACGTTGGCGTACTCGGACGGGTCCAGACGCTCTTTTTGCACGCGGGTCATGTACAGGATGTCAACTTCCGCCATCACTTCTTCAATAGAGCTGTGCAGACTCCAGGCGATCCC‑TTTTCATCAAGCATATCCAGAATg                                                                                                                             <  2:59011/150‑1 (MQ=255)
                       tGCGCTTTCACGTTGGCGTACTCGGACGGGTCCAGACGCTCTTTTTGCACGCGGGTCATGTACAGGATGTCAACTTCCGCCATCACTTCTTCAATAGAGCTGTGCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTAt                                                                                                                          <  2:636510/150‑1 (MQ=255)
                                                   ggTCCAGACGCTCTTTTTGCACGCGGGTCATGTACAGGATGTCAACTTCCGCCATCACTTCTTCAATAGAGCTGTGCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGa                                                                                              >  2:145867/1‑150 (MQ=255)
                                                   ggTCCAGACGCTCTTTTTGCACGCGGGTCATGTACAGGATGTCAACTTCCGCCATCACTTCTTCAATAGAGCTGTGCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGa                                                                                              >  1:583162/1‑150 (MQ=255)
                                                            gCTCTTTTTGCACGCGGGTCATGTACAGGATGTCAACTTCCGCCATCACTTCTTCAATAGAACTGTGCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTaaa                                                                                     >  2:428957/1‑150 (MQ=255)
                                                               ctTTTTGCACGCGGGTCATGTACAGGATGTCAACTTCCGCCATCACTTCTTCAATAGAGCTGTGCAGACTCCAGGCGCTCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAAc                                                                                  <  1:439146/150‑1 (MQ=255)
                                                               ctTTTTGCACGCGGGTCATGTACAGGATGTCAACTTCCGCCATCACTTCTTCAATAGAGCTGTGCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAAc                                                                                  >  2:174967/1‑150 (MQ=255)
                                                                 ttttGCACGCGGGTCATGTACAGGATGTCAACGTCCGCCATCACTTCTTCAATAGAGCTGTGCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGGCCGGCGCGATGAAGTAAAAACgg                                                                                <  1:196021/150‑1 (MQ=255)
                                                                            ggTCATGTACAGGATGTCAACTTCCGCCATCACTTCTTCAATAGAGCTGTGCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGaa                                                                     >  2:109779/1‑150 (MQ=255)
                                                                             gTCATGTACAGGATGTCAACTTCCGCCATCACTTCTTCAATAGAGCTGTGCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAAc                                                                    <  2:186510/150‑1 (MQ=255)
                                                                                aTGTACAGGATGTCAACTTCCGCCATCACTTCTTCAATAGAGCTGTGCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTCAAACCGGTTGCCGTCGAACTTc                                                                 >  1:234805/1‑150 (MQ=255)
                                                                                aTGTACAGGATGTCAACTTCCGCCATCACTTCTTCAATAGAGCTGTGCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTCCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTc                                                                 >  2:673886/1‑150 (MQ=255)
                                                                                 tGTACAGGATGTCAACTTCCGCCATCACTTCTTCAATAGAGCTGTGCAGACGCCAGGCGATCCCTTTTTCAGCAAGCATATCCAGAATGTATTGCGGCGTTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTAGAACATCg                                                                >  1:574717/1‑150 (MQ=255)
                                                                                   tACAGGATGTCAACTTCCGCCATCACTTCTTCAATAGAGCTGTGCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCt                                                              <  1:427742/150‑1 (MQ=255)
                                                                                       ggATGTCAACTTCCGCCATCACTTCTTCAATAGAGCTGTGCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACg                                                          >  2:435367/1‑150 (MQ=255)
                                                                                                   ccGCCATCACTTCTTCAATAGAGCTGTGCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAggg                                              <  2:131566/150‑1 (MQ=255)
                                                                                                     gCCGTCACTTCTTCACTCGAGCTGTGCAGACTCCAGGCGGTCCCTTTTTCAGCAAGCATATCCAGACTGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAg                                            <  1:10021/150‑1 (MQ=255)
                                                                                                     gCCATCACTTCTTCAATAGAGCTGTGCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAg                                            <  2:622062/150‑1 (MQ=255)
                                                                                                     gCCATCACTTCTTCAATAGAGCTGTGCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAg                                            <  2:508834/150‑1 (MQ=255)
                                                                                                     gCCATCACTTCTTCAATAGAGCTGTGCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAg                                            <  1:174468/150‑1 (MQ=255)
                                                                                                     gCCATCACTTCTTCAATAGAGCTGTGCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAg                                            <  2:631121/150‑1 (MQ=255)
                                                                                                     ccaatgccttctgCGATAGAGCTGTGCAGAGTCGATGCGATCCCTTTTTCATCCAGCCTGACCAGACTGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAg                                            <  1:297881/143‑1 (MQ=255)
                                                                                                          cACTTCTTCAATAGGGCTGTGCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAAc                                       <  1:176139/150‑1 (MQ=255)
                                                                                                          cACTTCTTCAATAGAGCTGTGCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAAc                                       <  2:431859/150‑1 (MQ=255)
                                                                                                            cttcttCAATAGAGCTGTGCAGACTCCAGGCGATCCCCTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACgg                                     >  2:413963/1‑150 (MQ=255)
                                                                                                            cttcttCAATAGAGCTGTGCAGACTCCAGGCGATCCCCTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACgg                                     >  2:656417/1‑150 (MQ=255)
                                                                                                              tcttcAATAGAGCTGTGCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTg                                   <  2:39180/150‑1 (MQ=255)
                                                                                                               cttcAATAGAGCTGTGCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGc                                  >  2:321817/1‑150 (MQ=255)
                                                                                                               cttcAATAGAGCTGTGCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGc                                  >  2:253096/1‑150 (MQ=255)
                                                                                                               cttcAATAGAGCTGTGCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGc                                  >  2:577742/1‑150 (MQ=255)
                                                                                                               cttcAATAGAGCTGTGCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCAAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGc                                  >  2:459330/1‑150 (MQ=255)
                                                                                                                  cAATAGCGCTGTGCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGc                               >  1:259132/1‑150 (MQ=255)
                                                                                                                  cAATAGAGCTGTGCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGc                               >  2:12522/1‑150 (MQ=255)
                                                                                                                  cAATAGAGCTGTGCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGc                               >  2:62446/1‑150 (MQ=255)
                                                                                                                            gtgCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCaa                     >  1:192053/1‑149 (MQ=255)
                                                                                                                             tgCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAgg                    <  2:449215/150‑1 (MQ=255)
                                                                                                                                   cTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCaccg              >  1:383452/1‑149 (MQ=255)
                                                                                                                                     ccAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCaccaac            <  1:320908/150‑1 (MQ=255)
                                                                                                                                      cAGGCGATCCCCTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCaccaacc           <  1:654487/150‑1 (MQ=255)
                                                                                                                                      cAGGCGATCCCCTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCaccaacc           <  1:412781/150‑1 (MQ=255)
                                                                                                                                               ccTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACg  <  2:384570/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GATCGCTGGCGCGAAGAACAAACTGCGCTTTCACGTTGGCGTACTCGGACGGGTCCAGACGCTCTTTTTGCACGCGGGTCATGTACAGGATGTCAACTTCCGCCATCACTTCTTCAATAGAGCTGTGCAGACTCCAGGCGATCCCTTTTTCATCAAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACG  >  REL606/4454910‑4455202

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