Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,875,632 (T)7→8 intergenic (‑66/+287) glmU ← / ← atpC bifunctional N‑acetylglucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyltransferase/F0F1 ATP synthase subunit epsilon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,875,6251.T91.9% 219.1 / 10.5 62intergenic (‑59/+294)glmU/atpCbifunctional N‑acetylglucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyltransferase/F0F1 ATP synthase subunit epsilon
Reads supporting (aligned to +/- strand):  ref base . (2/3);  new base T (30/27);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAG  >  REL606/3875552‑3875691
                                                                          |                                                                        
aGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtt                                                                       <  1:634045/76‑1 (MQ=255)
 gCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttt                                                                     <  2:1898896/76‑1 (MQ=255)
 gCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttt                                                                     <  2:829947/76‑1 (MQ=255)
    ttatcCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtttttt                                                               >  1:1880082/1‑76 (MQ=255)
     tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttttt                                                              >  2:1539305/1‑76 (MQ=255)
       ttCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttttttc                                                            >  1:416733/1‑74 (MQ=255)
             aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTg                                                      <  2:846339/76‑1 (MQ=255)
                 gcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAAt                                                  <  2:839211/76‑1 (MQ=255)
                   gTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATgg                                                >  1:871474/1‑76 (MQ=255)
                     ccTGACTGTAATTTGAGAACGAATTTAAATCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGAc                                              <  2:678117/76‑1 (MQ=255)
                     ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCGCCTTGAAAAAACTACATTTATTTTCATCGTGAAGTgggc                                              <  2:1971321/76‑3 (MQ=255)
                     ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATGGTGAAATGGAc                                              <  1:68454/76‑1 (MQ=255)
                     ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGAc                                              <  2:1391776/76‑1 (MQ=255)
                      cTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACa                                             <  2:1038915/76‑1 (MQ=255)
                       tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACag                                            >  1:630523/1‑76 (MQ=255)
                       tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGGGAAATGGACag                                            <  1:50776/76‑1 (MQ=255)
                        gACTGTAATTAGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATTGTGAAATGGACaga                                           >  1:169440/1‑76 (MQ=255)
                            gTAATTTGAGAANGANTTTAAACCGCTTCACCTTNAAAAAAATACATTTTTTTTCATCGTGAAATGGACAGAGGAt                                       <  1:1665932/76‑1 (MQ=255)
                            gTAATTTGAGAACGAATTTCAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGAt                                       >  1:647291/1‑76 (MQ=255)
                            gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGAt                                       >  1:230586/1‑76 (MQ=255)
                            gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGAt                                       <  1:861498/76‑1 (MQ=255)
                             tAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGAtn                                      <  1:1543004/76‑2 (MQ=255)
                              aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaa                                     >  2:1077671/1‑76 (MQ=255)
                                tttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAt                                   <  1:1005006/76‑1 (MQ=255)
                                tttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAt                                   >  1:1026807/1‑76 (MQ=255)
                                 ttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACGTATTTTTACAACGTGAAATGGAGAGANGATAAAtt                                  >  2:341396/1‑76 (MQ=255)
                                 ttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAtt                                  <  1:497729/76‑1 (MQ=255)
                                  tGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTg                                 <  1:976477/76‑1 (MQ=255)
                                      aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCa                             <  1:877572/76‑1 (MQ=255)
                                      aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCa                             <  1:1531277/76‑1 (MQ=255)
                                        cGAATTTAAAGCGCTTCGCCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGTATCAATTGTTCAAt                           >  1:163009/1‑76 (MQ=255)
                                         gAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATa                          >  1:1053373/1‑76 (MQ=255)
                                         gAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATa                          <  1:1240991/76‑1 (MQ=255)
                                            tttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACg                       <  1:596199/76‑1 (MQ=255)
                                            tttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACg                       >  2:1417176/1‑76 (MQ=255)
                                             ttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACgg                      <  2:1772500/76‑1 (MQ=255)
                                             ttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCCATCGTGAAATGGACAGAAGATAAATTGTTCAATAACgg                      >  2:1632389/1‑76 (MQ=255)
                                              tAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa                     >  1:809015/1‑76 (MQ=255)
                                              tAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa                     >  1:1417586/1‑76 (MQ=255)
                                              tAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa                     >  1:1491446/1‑76 (MQ=255)
                                               aaaCCGCTTCACCTTGAAACAACTACATTATTTATCATAGTGAAGTGGGTAGAGGATAAATTGTTCACTAACAGAt                    >  1:1972448/1‑76 (MQ=255)
                                                aaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt                   <  1:1016583/76‑1 (MQ=255)
                                                aaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt                   >  2:1281194/1‑76 (MQ=255)
                                                aaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt                   <  2:209639/76‑1 (MQ=255)
                                                aaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt                   <  2:1901768/76‑1 (MQ=255)
                                                aaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt                   <  1:609603/76‑1 (MQ=255)
                                                aaCCGCTTCACCTTGAAAAAACTACATTTTATTTCATCGCGAAATGGACACAGGATAAATTGTTCAATAACGGAtt                   <  1:26951/76‑1 (MQ=255)
                                                 aCCGCTTCATCGTGAAAAAACTACATTTTTTTGCATCGTGAAATGGACAGAGGATAGCGTGTTCAATAACGGATTa                  <  1:23726/76‑1 (MQ=255)
                                                 aCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTa                  <  2:920401/76‑1 (MQ=255)
                                                  ccGCTTCACCTTGAAAAAACTACATTTTTTTTCATCATGAAATGGTCAGAGGATAAATTGTTCAATAACGGATTAt                 >  1:38458/1‑76 (MQ=255)
                                                   cGCTTAACATTGAANAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATc                >  1:136276/1‑76 (MQ=255)
                                                     cTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATccc              >  1:1384928/1‑76 (MQ=255)
                                                      ttCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcccc             >  1:703512/1‑76 (MQ=255)
                                                      ttCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcccc             <  2:511534/76‑1 (MQ=255)
                                                      gtCACCTTGAAAAAACTACATTTTTCTTTATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcccc             <  1:325553/75‑1 (MQ=255)
                                                       tCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg            >  1:1206011/1‑76 (MQ=255)
                                                       tCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg            >  2:1816671/1‑76 (MQ=255)
                                                         aCCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGGTCAATAACGGATTATCCCCGCt          <  2:802260/76‑1 (MQ=255)
                                                          ccTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTg         <  1:1429627/76‑1 (MQ=255)
                                                          ccTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTg         <  2:546952/76‑1 (MQ=255)
                                                           cTTGAAAAAACTACNTTTTTTTTCATCGTGAAATGGAAAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa        >  1:5954/1‑76 (MQ=255)
                                                           cTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa        >  1:145421/1‑76 (MQ=255)
                                                           cTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa        >  2:627906/1‑76 (MQ=255)
                                                             tGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGCTTATCCCCGCTGACg      >  1:1694618/1‑76 (MQ=255)
                                                              gAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGa     >  1:1046227/1‑76 (MQ=255)
                                                               aaaaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCACACACGGATTATCCCCGCTGACGaa    >  2:1421227/1‑76 (MQ=255)
                                                               aaaaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaa    >  2:1956427/1‑76 (MQ=255)
                                                               aaaaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaa    >  1:898438/1‑76 (MQ=255)
                                                                aaaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaaa   >  2:689269/1‑76 (MQ=255)
                                                                 aaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAg  >  2:1042375/1‑76 (MQ=255)
                                                                          |                                                                        
AGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAG  >  REL606/3875552‑3875691

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