Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,134,767 (A)5→4 intergenic (‑59/‑95) argE ← / → argC acetylornithine deacetylase/N‑acetyl‑gamma‑glutamyl‑phosphate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,134,7630A.91.0% 262.5 / 24.4 67intergenic (‑55/‑99)argE/argCacetylornithine deacetylase/N‑acetyl‑gamma‑glutamyl‑phosphate reductase
Reads supporting (aligned to +/- strand):  ref base A (5/1);  new base . (33/28);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

GGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCA  >  REL606/4134689‑4134833
                                                                          |                                                                      
ggCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATaa                                                                       <  2:1615256/76‑1 (MQ=255)
 gCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATaaa                                                                      >  2:827484/1‑76 (MQ=255)
     tAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATaaaatct                                                                  <  1:377920/76‑4 (MQ=255)
     tAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATa                                                                  >  2:1991665/1‑76 (MQ=255)
        tttGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaac                                                              >  1:1438628‑M1/1‑71 (MQ=255)
         ttGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacg                                                             <  1:1422292‑M1/76‑7 (MQ=255)
          tGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgt                                                            >  1:80172‑M1/1‑69 (MQ=255)
            tttttCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAAc                                                           >  1:1348900/1‑76 (MQ=255)
            tttttCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttg                                                          <  1:617602‑M1/76‑10 (MQ=255)
            tttttCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttg                                                          >  1:293631‑M1/1‑67 (MQ=255)
            tttttCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttg                                                          <  1:1639394‑M1/76‑10 (MQ=255)
             ttttCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttga                                                         <  1:1144604‑M1/76‑11 (MQ=255)
              tttCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttgag                                                        <  1:748722‑M1/76‑12 (MQ=255)
               ttCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagc                                                       >  2:1497434‑M1/1‑64 (MQ=255)
                 cATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagcgt                                                     >  1:477330‑M1/1‑62 (MQ=255)
                 cATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagcgt                                                     >  1:301569‑M1/1‑62 (MQ=255)
                 cATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagcgt                                                     <  1:1766907‑M1/76‑15 (MQ=255)
                 cATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagcgt                                                     >  1:275971‑M1/1‑62 (MQ=255)
                  aTTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagcgta                                                    >  1:1376599‑M1/1‑61 (MQ=255)
                   ttgttgACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaa                                                   <  1:1390‑M1/76‑17 (MQ=255)
                   ttgttgACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaa                                                   <  1:429465‑M1/76‑17 (MQ=255)
                   ttgttgACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaa                                                   <  1:335639‑M1/76‑17 (MQ=255)
                   ttgttgACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaa                                                   >  2:75044‑M1/1‑60 (MQ=255)
                     gttgACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaata                                                 <  1:1101254‑M1/76‑19 (MQ=255)
                       tgACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaa                                               >  1:907928‑M1/1‑56 (MQ=255)
                        gACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaa                                              >  1:643740‑M1/1‑55 (MQ=255)
                         acacacCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaac                                             >  2:815485‑M1/1‑54 (MQ=255)
                          cacacCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaacc                                            <  1:118760‑M1/76‑24 (MQ=255)
                           acacCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaaccc                                           >  1:1907101‑M1/1‑52 (MQ=255)
                           acacCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaaccc                                           >  1:1067767‑M1/1‑52 (MQ=255)
                               ctctGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaacccacca                                       <  1:1956965‑M1/76‑29 (MQ=255)
                                tctGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaacccaccag                                      <  2:623776‑M1/76‑30 (MQ=255)
                                 ctGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaacccaccagc                                     >  1:250065‑M1/1‑46 (MQ=255)
                                  tGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaacccaccagcc                                    <  1:1416440‑M1/76‑32 (MQ=255)
                                    gTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaacccaccagccgt                                  >  2:320827‑M1/1‑43 (MQ=255)
                                      cATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGcc                                 <  1:1198/76‑1 (MQ=255)
                                      cATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaacccaccagccgtaa                                >  1:388369‑M1/1‑41 (MQ=255)
                                      cATGATAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaacccaccagccgtaa                                <  2:645072‑M1/76‑36 (MQ=255)
                                        tGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGt                               >  1:556200/1‑76 (MQ=255)
                                          aTAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaacccaccagccgtaaggtg                            >  1:829617‑M1/1‑37 (MQ=255)
                                           tAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaacccaccagccgtaaggtga                           >  1:788739‑M1/1‑36 (MQ=255)
                                           tAGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaacccaccagccgtaaggtga                           >  1:788743‑M1/1‑36 (MQ=255)
                                            aGTATCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaacccaccagccgtaaggtgaa                          >  2:1199298‑M1/1‑35 (MQ=255)
                                               aTCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaacccaccagccgtaaggtgaatgt                       >  2:1099917‑M1/1‑32 (MQ=255)
                                               aTCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaacccaccagccgtaaggtgaatgt                       <  2:1302450‑M1/76‑45 (MQ=255)
                                               aTCAATATTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaacccaccagccgtaaggtgaatgt                       >  1:1637385‑M1/1‑32 (MQ=255)
                                                  aaTATTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaccccaccagccgtaaggtgaatgtttt                    <  2:1035717‑M1/76‑48 (MQ=255)
                                                  aaTATTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaacccaccagccgtaaggtgaatgtttt                    <  2:1273686‑M1/76‑48 (MQ=255)
                                                  aaTATTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaacccaccagccgtaaggtgaatgtttt                    <  1:824720‑M1/76‑48 (MQ=255)
                                                  aaTATTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaacccaccagccgtaaggtgaatgtttt                    >  2:1803159‑M1/1‑29 (MQ=255)
                                                   atatTCATTCAGTATTTATGAATAAAATctaacgttgagcgtaataaaccccaccagccgtaaggtgaatgtttta                   <  2:114622‑M1/76‑49 (MQ=255)
                                                   atatTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaaccccccagccgtaaggtgaatgtttta                   <  1:1355883‑M1/76‑49 (MQ=255)
                                                   atatTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaacccaccagccgtaaggtgaatgtttta                   <  2:1485962‑M1/76‑49 (MQ=255)
                                                    tatTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaacccaccagccgtaaggtgaatgttttac                  >  1:940424‑M1/1‑27 (MQ=255)
                                                    tatTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaacccaccagccgtaaggtgaatgttttac                  >  2:32124‑M1/1‑27 (MQ=255)
                                                    tatTCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaacccaccagccgtaaggtgaatgttttac                  >  1:16479‑M1/1‑27 (MQ=255)
                                                     atTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGttt                  >  1:1418342/1‑76 (MQ=255)
                                                       tCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaccccaccagccgtaaggtgaatgttttacgtt               <  1:1497434‑M1/76‑53 (MQ=255)
                                                       tCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaccccaccagccgtaaggtgaatgttttacgtt               <  2:434812‑M1/76‑53 (MQ=255)
                                                       tCATGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaacccaccagccgtaaggtgaatgttttacgtt               >  2:1684534‑M1/1‑24 (MQ=255)
                                                         aTGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaacccaccagccgtaaggtgaatgttttacgttta             >  1:1420459‑M1/1‑22 (MQ=255)
                                                         aTGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaacccaccagccgtaaggtgaatgttttacgttta             >  1:1420442‑M1/1‑22 (MQ=255)
                                                         aTGCAGTATTTATGAATAAAATctaacgttgagcgtaataaaacccaccagccgtaaggtgaatgttttacgctta             >  1:1420456‑M1/1‑22 (MQ=255)
                                                           gCAGTATTTATGAATAAAATctaacgttgagcgtaataaaacccaccagccgtaaggtgaatgttttacgtttaac           <  2:173629‑M1/76‑57 (MQ=255)
                                                             aGTATTTATGAATAAAATctaacgttgagcgtaataaacccccccagccgtaaggtgaatgttttacgtttaacct         <  2:491004‑M1/76‑59 (MQ=255)
                                                                 tttATGAATAAAATctaacgttgagcgtaataaaacccaccagccgtaaggtgaatgttttacgtttaacctggca     >  1:3642‑M1/1‑14 (MQ=255)
                                                                  ttATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTg     >  1:1748125/1‑76 (MQ=255)
                                                                  ttATGAATAAAATctaacgttgagcgtaataaaccccaccagccgtaaggtgaatgttttacgtttaacctggcaa    <  2:618444‑M1/76‑64 (MQ=255)
                                                                  ttATGAATAAAATctaacgttgagcgtaataaaacccaccagccgtaaggtgaatgttttacgtttaacctggcaa    >  1:1478063‑M1/1‑13 (MQ=255)
                                                                   tATGAATAAAATctaacgttgagcgtaataaaacccaccagccgtaaggtgaatgttttacgtttaacctggcaac   <  2:1818221‑M1/76‑65 (MQ=255)
                                                                    aTGAATAAAATctaacgttgagcgtaaaaaaccccaccagccgtaaggtgaatgttttacgtttaacctggcaacc  <  2:424763‑M1/76‑66 (MQ=255)
                                                                          |                                                                      
GGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCA  >  REL606/4134689‑4134833

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: